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An Improved High-Frequency 
Radiotelephone System Featuring 
Constant Net Loss Operation 


By J. M. FRASER, H. H. HAAS and M. G. SCHACHTMAN 
(Manuscript received December 9, 1966) 


A system ts described in which a special type of syllabic compandor 
and other features offer an improved solution to the problems of using HF 
radio facilities for overseas telephone trunks. A major disadvantage of 
existing systems is the use of voice-operated gates to prevent singing and 
echo effects. These gates interfere with the free flow of conversation. The 
new system allows stable full-duplex operation as in conventional telephone 
circutis because a nearly constant circuit loss is maintained between the 
trunk terminals. 

Significant improvements in circuit quality result from full-duplex 
operation and from a reduction in the effects of radio noise and interference 
by compandor action. Preliminary results of a field trial between New 
York and Buenos Aires have been highly favorable. 


I. INTRODUCTION 


The first transatlantic telephone cable was laid in 1956. For almost 
thirty years prior to that time, transoceanic telephone communication 
was provided almost exclusively via high-frequency radiotelephone 
circuits. Numerous submarine cables have since been laid, and HF 
radio is no longer used for telephone service over routes of heaviest 
traffic, e.g., New York-London, where large numbers of cable and 
satellite circuits exist. Many long transoceanic routes, however, are 
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not served by submarine cables or satellites; direct telephone service 
between the United States and forty foreign countries and areas (see 
Fig. 1) is still provided by means of HF radio facilities exclusively. 

Although a change in the composition of transoceanic telephone 
facilities is expected as satellites become more fully utilized, it is likely 
that the HF services will hold their own for economic reasons, partic- 
ularly for small circuit groups. 

For some time the HF radio spectrum, a rather limited resource to 
begin with, has been nearing the point of saturation. Thus, as new 
submarine cable and land* facilities make it possible to suspend exist- 
ing HF routes, the frequencies are reassigned to provide service to new 
areas, typically to the developing nations of the world. Thus, the total 
number of HF circuits is expected to remain about the same well into 
the communication satellite era. 

Since high frequencies propagate via the ionosphere, which is in a 
continual state of change, HF transmission is highly variable and 
requires special equipment and procedures to cope with the problems 
created by this variability. At times, HF signals (say, 4 to 27 MHz 
for overseas service) propagate easily over great distances; at other 
times, the medium will not propagate a given frequency at all. Any 
condition between these extremes is possible. 

The ionospheric medium is a difficult and erratic one, subject to 
many propagation anomalies. Moreover, unique interconnection prob- 
lems exist between point-to-point HF radiotelephone circuits and con- 
ventional land plant. Because of the unstable nature of the HF radio 
transmission medium, the combination of the 2-wire part of the land 
plant and a “4-wire” radio circuit presents problems that are not 
encountered, for example, in 4-wire cable or microwave radio circuits. 
Any interconnection of 4-wire and 2-wire circuits creates potential 
singing and echo problems. In the land plant, this problem is solved by 
operating the circuit at a net loss sufficient to prevent singing and to 
render echoes negligible.t Because such circuits are stable, any allow- 
ance for variations in loss is small, and the mode of operation may be 
characterized as “constant net loss.” If HF radiotelephone circuits 
were operated in this mode without special controls, the allowance for 
variations would, in gencral, have to be rather large. If the design net 
loss were sufficient for stability with minimum transmission loss in the 
radio portion of the circuit, even the average loss would result in low 


* HF radio is widely used for intracontinental communication in some less de- 
veloped parts of the world. 
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received volume. The high loss occurring during deep fading would, at 
times, make the received volume inaudible. This would be a very in- 
efficient mode of operation. 

Input speech volumes from connecting plant vary over a wide range, 
at least 1000:1 in power. Linear transmission of this full range is an 
inefficient use of radio transmitter load capacity. 

Because. signal-to-noise and signal-to-interference ratios are fre- 
quently marginal, fairly high power is required for reasonable relia- 
bility. For example, at least a 4-kilowatt output might be required on 
a given path for adequate reception of the weakest talker. Then a 
4-megawatt transmitter capacity would be required in a linear system 
for the strongest talker. The use of 4-megawatts perhaps one percent 
of the time would be uneconomical. In addition, the unnecessary radia- 
tion of excessive power would aggravate interference problems. As a 
counterexample, if 4 kilowatts were the maximum available to the 
strongest talker, the weakest talker would develop an output of only 
4 watts and might be lost in the noise. These examples make it obvious 
that in order to minimize transmitter power and still accommodate 
the large range of speech volumes, some form of volume regulation is 
necessary. Volume regulation, which varies the gain of the circuit, 
compounds the problem of singing stability and echo effects discussed 
above. 

The basic HF radiotelephone transmission problem may be sum- 
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marized at this point as one of continuously maximizing the signal- 
to-noise ratio under varying conditions while at the same time pre- 
venting a net gain around the circuit that would allow singing to 
occur. The conventional solution to this problem? comprises at least 
the following two elements: 


(t) VODAS, or Voice Operated Device Anti-Singing, which makes the 
circuit one-way at a time. With this arrangement, singing and talker 
echoes are avoided. 

(it) VOGAD, or Voice Operated Gain Adjusting Device, which loads 
the radio transmitter efficiently. 


In one form or another, these elements have been widely used through- 
out the world from the beginning of the service up to the present time. 
Fig. 2 is a simplified diagram of the radio terminal equipment used 
in overseas gateways of the Bell System. In addition to the transmit- 
ting VOGAD and VODAS, there is a receiving VOGAD and a “noise 
reducer.” Briefly, the terminal operates as follows. The receiving path 
is normally enabled and the transmitting path disabled. Outgoing 
speech, regulated to constant volume by the transmitting VOGAD, 
operates the transmitting speech detector, which enables the transmit 
path and disables the receiving path. When the speech train ends, the 
VODAS reverts to the receiving mode. The receiving VOGAD provides 
fading regulation of incoming speech, which has already been coarsely 
compensated by automatic gain control in the receiver. The noise re- 
ducer is a voice-activated expandor capable of moderate noise and 
interference reduction in the speech gaps if radio conditions are at 
least fair. The receiving detector, by inhibiting the transmitting detec- 
tor, prevents echoes of the received speech from switching the circuit 
to the transmitting mode. The receiving detector is provided with a 
sensitivity adjustment that may be set by a “technical operator.” 
Although its wide use for a long period attests to the soundness of 
the present terminal, a number of disadvantages have become ap- 
parent. Foremost is the fact that even under perfect transmission con- 
ditions, a one-way-at-a-time circuit inhibits the smooth flow of con- 
versation; in fact, when double-talking occurs, a significant amount 
of speech may be lost. This effect is sometimes called “lock-out”. Also, 
because of the difficulty in differentiating between speech and noise, 
some “clipping” by the speech detectors occurs under the best condi- 
tions, particularly of the initial portions of speech having low energy. 
VODAS designs employing a receiving detector for echo protection 
are particularly vulnerable to operation on high received noise or in- 
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terference. Such false operation is called “lock-up” and makes the 
circuit completely unusable as long as it persists. Frequent adjust- 
ment of the receiving detector sensitivity and other controls is re- 
quired if transmission quality is to be optimized with changing radio 
conditions. 


II. AN IMPROVED RADIOTELEPHONE SYSTEM 


The HF radiotelephone transmission problems can be solved by a 
method other than the conventional one described. Experimental sys- 
tems based on general principles known for some time? * recently have 
been investigated by the Bell Telephone Laboratories and others.» ® * 8 
The major disadvantages of the existing system are eliminated and 
significant improvements are realized in the new type of system. This 
paper describes the experimental system investigated by the Bell Tele- 
phone Laboratories, called the constant net loss (CNL) system. 

The principles of the CNL system may be explained with reference 
to Fig. 3, which illustrates a circuit equipped at both ends with a 
new type of terminal. Input speech is applied to a “complete compres- 
sor,’ which gives constant output volume over a wide range of inputs. 
The compressor smooths out the syllabic changes in speech loudness 
in order to fully load the radio transmitter and to optimize the re- 
ceived signal-to-noise ratio on a syllable-to-syllable basis. Thus, it is 
more effective than the slow-acting VOGAD, which regulates differ- 
ences between talkers but does little to the syllabic variations of a 
given talker. The control signal from the first stage of the compressor, 
a low-frequency (0-100 Hz) analog of the syllabic variations of the 
input speech, feeds the two stages to give complete compression as ex- 
plained in Appendix A. In addition, the control signal frequency- 
modulates a subcarrier, which is transmitted over the radio link along 
with the compressed speech, but in a separate narrow-band channel. 
Fig. 4 shows the frequency allocations of a single voice channel with 
its control channel. 

The control signal is used at the receiving end to control an ex- 
pandor, which restores the original variations in loudness. The FM 
control channel has sufficient margin to assure that even under severe 
noise, interference, and fading conditions, the expandor will properly 
“track” the compressor. The blocks labeled LOG and ANTILOG in the 
control channel comprise an instantaneous compandor® (nonlinear 
compressor and expandor) whose function is to make the control 
signal less susceptible to noise and interference in the control channel. 

Although restoration of the original speech loudness variations by 
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Fig. 4— One-channel operation. 


an expandor controlled from the sending end is of some benefit in 
itself, it has two other more essential purposes: (7) It maintains an 
overall constant net loss by making the expandor loss track the com- 
pressor gain; (2) it very effectively mutes radio noise and inter- 
ference during the gaps between syllables when the expandor loss is a 
maximum. A substantial subjective improvement accrues because the 
noise is reduced when it is most noticeable. Expandor tracking error 
caused by noise or interference in the control channel is minimized by 
the narrow bandwidth, FM carrier power, and other parameters of 
the control channel. The complementary compressor/expandor action 
via the control channel results in constant loss in the absence of fad- 
ing. A stable loss around the outgoing and return paths allows the 
circuit to be set up “full-duplex” without a VODAS. 

In the CNL terminal, fading variations that remain after the action 
of the receiver automatic gain control are absorbed by the fading 
regulator preceding the expandor. The fading regulator (which is 
also a syllabic complete compressor) operates independently of the 
control channel and somewhat more slowly than the transmitting com- 
pressor. The dynamic actions of the transmitting compressor, fading 
regulator, and expandor are such that, although gains and losses of 
these devices and of the transmission medium vary within the circuit, 
the net overall loss is approximately constant. 

Since the CNL system has no VODAS with its inherent problems 
of clipping, lock-out, and receiving detector lock-up, a smoother flow 
of conversation results. Its compandor action and fading regulation 
are more effective in optimizing S/N (signal-to-noise) and S/I 
(signal-to-interference) ratios than the uncoordinated actions of the 
VOGADs and noise reducer of present terminals. An additional ad- 
vantage is that no in-service adjustments of the terminal by technical 
operators are necessary, as is the case in the present terminals. A 
summary of the principal characteristics of experimental CNL systems 
is given in Table I. 
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TasBLE I—SummMary oF CNL HF RapioTeLEPHONE 


System Capacity 


Transmitting Com- 
pressor Operating 
Range 


Speech Bandwidth 


Control Signal 
Allocations and 
Power 


Time Constants 


Fading Regulator 


FM Control Signal 


Deviation /Loss 
Ratio 





TERMINAL CHARACTERISTICS 


Four speech channels plus four control channels can be 
accommodated on one radio carrier in conjunction with 
conventional 12-kHz radio and channelizing equipment. 


The range of input power over which compressor output is 
essentially constant is -+-8 dBm0 to —32 dBm0. Below 
—32 dBm0, the compressor acts as a linear amplifier; 
above +8 dBm0, some clipping of a sinewave occurs. 


The nominal speech bandwidth is 250 to 2750 Hz for 
channels next to the carrier, 250 to 2950 Hz for outboard 
channels. 


(¢) The FM control signal associated with an inboard 
channel occupies an 180-Hz band centered at 
2910 Hz. With no system input, the subcarrier 
frequency is 2850 Hz; maximum test tone input 
(+8 dBm0) deviates the subcarrier upward 120 Hz 
to 2970 Hz. 

(ii) The outboard channel, 180 Hz wide, is centered on 
3140 Hz, and the corresponding subcarrier fre- 
quency excursion is from 3200 Hz downward to 
3080 Hz. 

(iit) The FM subcarrier power in the composite signal is 
applied to the radio transmitter at a level 16 dB 
below its peak envelope power rating and 6 dB 
below a test tone in the speech channel, the test 
tone being within the constant output range of 
the transmitting compressor. 


The transmitting compressor has an attack time of 3 msec 
and a recovery time of 13.5 msec. The overall effective 
RC time constant of the compandor system is 20 msec. 


The fading regulator maintains essentially constant output 
over an input range 20 dB below and 10 dB above the 
nominal no-fade level. The attack and recovery times 
are 12 msec and 54 msec, respectively. These values are 
4 times the corresponding transmitting compressor 
time constants. 


The deviation/loss ratio is the transfer constant that 
relates the deviation of the FM subcarrier to the ex- 
pandor loss variations and is the best measure of the 
susceptability of the control channel to noise and fre- 
quency instability. The constant is 2 Hz per dB. 








The CNL terminal is compatible with existing channelizing equip- 
ment, privacy devices, radio transmitters, and receivers. The system 
uses no more bandwidth overall than present systems. Although the 
FM channel shares bandwith with the speech, the loss of speech band- 
width is slight. This is possible because of more efficient use of the 
available bandwidth by means of sharper filter cut-offs. In conjunction 
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with conventional 12-kHz radio and channelizing equipment, a full 
complement of four speech channels plus four control channels can be 
accommodated (Fig. 5). 

The control signal must share the available power of an existing 
transmitter with the compressed speech. Satisfactory performance is 
expected with control subcarrier magnitudes (see Table I) such that 
the total transmitter load is not significantly increased. 

There are several reasons for using FM to transmit the control in- 
formation. The foremost reason is that limiter action provides level 
compensation, making the control signal insensitive to fading (as 
regards purely level, rather than S/N variations). Also, bandwidth 
can be traded for S/N advantage. In addition, with FM, the de com- 
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Fig. 5 — Four-channel operation. 


ponent of the control signal may be preserved.* Finally, the narrow- 
band control signal is similar in many ways to narrow-band telegraph 
signals which have been transmitted successfully via HF radio for 
many years using FM or related methods. It is too early to predict 
the extent to which diversity would be useful or necessary in a CNL- 
type radiotelephone system, although it is widely used with HF radio 
telegraph systems. 

A successful trial of experimental CNL terminals has been con- 
ducted between New York and Buenos Aires by the American Tele- 
phone and Telegraph Company with the cooperation of the foreign 
correspondent. Conventional terminals on one circuit of a regular four- 
channel system were replaced with CNL terminals. A large number of 
test and demonstration calls were made including some comparisons 


* The control signal is derived by rectifying the speech output of the first com- 
pressor stage; itis therefore unipolar and has a de component. 
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between the CNL circuit and a conventional circuit operating in the 
same radio system. Radio conditions on this path generally are fair 
to good, neglecting propagation outages. The participants in these 
tests were largely persons familiar with both cable and conventional 
radio overseas circuit performance. In the judgment of most of the 
participants, the CNL circuit quality approached that of a submarine 
cable circuit because of the compandor action and lack of VODAS 
impairments. 

Following the above demonstrations in July, 1966, the CNL- 
equipped circuit began an extended period of commercial service. 
Traffic records indicate that the single CNL circuit carried one quarter 
of all calls on the New York-Buenos Aires route, which has a total of 
12 circuits. The CNL circuit handled more calls than the next two 
most active conventional circuits combined. In the opinion of operat- 
ing personnel, calls that would have encountered considerable im- 
pairments or operating difficulties with conventional terminals were 
handled without customer complaint. During periods of poor radio 
conditions, the CNL circuit was frequently “commercial” when some, 
if not all, of the conventional circuits were unusable. 

Similar experimental terminals have becn developed and tested by 
the British General Post Office (GPO) and French Postes Et Tele- 
communications (PTT). The British and American Administrations 
submitted reports of their work on the new type of terminal to the 
plenary meeting of the International Radio Consultative Committee 
(CCIR) in Oslo in the summer of 1966. This report included essential 
design parameters and operating conditions, and was submitted as an 
initial contribution toward adoption of international compatibility 
standards for the new system. 


III. DESCRIPTION AND ANALYSIS OF IMPROVED SYSTEM 


3.1 System Application—Overall Block Diagram 


The block diagram Fig. 6 illustrates an arrangement of the experi- 
mental CNL equipment that was installed on a working radio circuit 
between Buenos Aires and New York. The equipment shown as 
shaded boxes is the same as is used with the present terminals. Other 
arrangements of signaling, privacy, channel shifters, and filters are 
possible; the best arrangements will depend on the characteristics of 
the specific equipment to be used on a given circuit. 

One independent sideband (ISB) HF transmitter normally carries 
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Fig. 6 — System block diagram. 
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& maximum of four channels, two on each sideband. The particular 
arrangement of Fig. 6 was designed, among other reasons, to facilitate 
building up multichannel CNL systems. In such systems, the voice 
outputs of two separate terminals are combined in a channel shifter; 
a single band-elimination filter at the output of the shifter produces 
a slot to accommodate two control channels. The control subcarrier 
of each terminal is transmitted through a bandpass filter of the ap- 
propriate frequency; the subcarriers from two terminals are then 
combined with the output of the band-elimination filter. The com- 
bination of the two speech and two control channels is applied to one 
sideband of an ISB transmitter. The ISB transmitter has provision 
for two additional circuits on the other sideband, as shown in Fig. 5. 
Since only one circuit was equipped with CNL terminals in the trial 
installation, the second, third, and fourth terminals indicated on Fig. 
6 were conventional VODAS terminals. 

The filter and channel restorer arrangement of the receiving side is 
essentially a duplicate of the transmitting side, with a single band- 
elimination filter for one pair of speech channels and bandpass filters 
for each of two control channels. The attenuation required to control 
crosstalk between the speech and the control channels is provided by 
the combination of transmitting and receiving filters. It will be noted 
from Fig. 6 that the control channel was combined with the speech 
channel after the channel shifter and split before the channel restorer. 
This approach was used so that the characteristics of particular chan- 
nel shifters and restorers would not affect the control channels. 

Signaling for the trial installation was based on manual operation 
of the radio telephone circuit, with a ring-down signal of 1000 Hz mod- 
ulated at 20 Hz inserted on the radio side of the terminal after the 
privacy. A 1000/20 signaling receiver, bridged after the fading regu- 
lator and before the privacy and expandor on the receive side, operated 
the appropriate switchboard equipment. Suitable trunk, monitoring, 
and control circuits similar to those found in VODAS terminals were 
provided. 

The trial installation includes a flat delay unit to equalize the delay 
of the speech path with respect to the separate narrowband control 
path. (A discussion of delay equalization is given in Section 3.2.4.) 

Terminals of the CNL type are inherently full-duplex and all equip- 
ment exclusive of the 4-wire terminating set, or hybrid, is arranged on 
a 4-wire basis. Fig. 6 shows connection to a 2-wire switchboard; 4-wire 
switchboards are also used. Privacy equipment must be 4-wire to pre- 
serve the two-way nature of the CNL principle. Either the inversion 
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or split-band type of privacy equipment may be used. In the split- 
band privacy, the speech is separated into five 550-Hz bands and 
transposed in frequency in any arrangement that provides reasonable 
privacy. With the CNL-type terminal, it may be necessary to restrict 
the possible transpositions to insure that the control channel slot is 
always taken out of the speech around 3000 Hz in order to minimize 
the effect of removing this slot from the speech band. 

An echo suppressor is required on all long circuits because of the 
inherent propagation time. The receiving split-type echo suppressor 
used on the Buenos Aires-New York circuit has a speech gate only in 
the transmit leg; the gate blocks transmitted speech only when speech 
of greater magnitude is present in the receiving leg. This type of sup- 
pressor was preferred because of its break-in properties and also be- 
cause it tends to reduce reradiation of echoes. 


3.2 Analysis of Factors Affecting Performance 


3.2.1 General Considerations 


As mentioned previously, allocations of available bandwidth and 
transmitter power must be made in a CNL-type system, in which the 
compressed speech and the information to control the expandor are 
transmitted via separate channels. The objective in making this alloca- 
tion is to obtain adequate margin against noise and interference in the 
control channel without encroaching substantially on the speech band- 
width and power. While the speech bandwidth is not critical, it is not 
a simple matter to determine what bandwidth should be allocated to 
the control signal. The procedure used in selecting the parameters of 
the experimental CNL system was to make a somewhat arbitrary al- 
location of bandwidth and power, then to investigate what further 
trade-offs could be made within these constraints. 

The control channel bandwidth selected was 180 Hz. The net speech 
bandwidths obtained were as indicated on Fig. 5 and Table I. The im- 
pairment due to loss of about 200 Hz of high-frequency speech energy 
of the inboard channel is slight. There is no loss of speech bandwidth 
to the outboard channel because the control channel in this case lies 
in what was originally guard space (between 3 and 3.25 kHz). 

In the experimental CNL system, the available power per channel 
was divided in the ratio 1:4, i.e., the FM subcarrier power was 6 dB 
below a 0-dBm0 voice-frequency test tone. The total load, including 
speech and subcarrier, was then only about 1 dB greater than the 
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speech power alone. It would be undesirable to put so much power in 
the subcarrier as to necessitate a reduction in the speech drive to the 
transmitter and a consequent reduction in signal-to-noise ratio at the 
receiver. 

For further discussion of factors affecting performance of the im- 
proved system, attention is drawn to the fact that the transmitting 
compressor and receiving expandor together comprise a syllabic com- 
pandor, which is closely related to more conventional devices of this 
class. In these devices, as well as in the CNL system, the transmission 
gain (or loss) is varied at syllabic rates in order to gain advantage 
against noise and interference. The conventional syllabic compandor 
has a compression ratio of 2:1, and uses the residual speech envelope 
of the compressed signal to derive the information with which to con- 
trol the expandor. Thus, it does not require a separate control chan- 
nel, although there are circumstances where such an arrangement could 
be advantageous. For example, a control channel could improve track- 
ing of such a compandor in the presence of high noise at the expandor 
input. 

Speech can be transmitted through a conventional compandor with 
low distortion using no more bandwidth than that of the original 
speech as long as the rate of gain variation, or “speed,” does not ex- 
ceed syllabic rates. Hence, the bandwidth, within limits, is not a con- 
sideration in establishing the speed of a conventional compandor. 

The following general requirements govern the speed of a syllabic 
compandor. (7) Compressor speed: a fast compressor is more effective 
than a slow compressor in raising weak speech syllables with respect 
to the noise before transmission. If this action is too abrupt, however, 
significant distortion is created. (11) Expandor speed: the speed of the 
expandor must match that of the compressor. If it is too slow, the ex- 
pandor may mutilate the initial parts of a syllable, or its loss may not 
be fully inserted during the gaps so that the noise reduction effect of 
the expandor suffers. 

Within the range of syllabic rates, the speed of a conventional com- 
pandor may be varied over a fairly broad range and still satisfy the 
foregoing requirements. Furthermore, the compandor noise improve- 
ment, 1.e., how well the compressor picks up weak speech before trans- 
mission and the expandor mutes the noise between syllables, is not 
highly sensitive to speed. 

The use of a separate control channel to control the expandor in the 
CNL system introduces an additional factor governing speed of the 
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CNL compandor; otherwise, the criteria governing speed would be 
the same as a conventional compandor. The additional factor is the 
noise performance of the control channel, which is strongly dependent 
on compandor speed. The speed of the compandor in the CNL system 
may be regarded as an independent parameter upon which two main 
categories of system performance depend. These are: (2) the control 
channel noise performance and (ii) the compandor noise improvement. 
The first degrades with increasing speed and the second improves with 
increasing speed; therefore, it is necessary to select parameters such 
that a reasonable balance is achieved. 

The principal elements of the CNL system affecting this balance 
can be represented by the model of Fig. 7. The control signal is de- 
rived from the output of the first compressor variolosser by rectifying 
a portion of the speech voltage at this point. The unfiltered output 
undergoes an initial smoothing, or prefiltering, in the network RoCo. 
It then branches into two paths, one feeding back through the low pass 
amplifier (R,C,) to control both variolossers of the compressor and 
the other feeding forward via the control channel and the low pass 
amplifier (ReC.2) to control the variolossers of the expandor. RoCo band 
limits the spectrum of the modulating signal at the input to the con- 
trol channel. The primary role of RsCe is that of a post-detection noise 
filter. In addition to their filtering functions, the networks RoCo and 
R.Cz together establish the time response of the expandor. The FM 
channel itself can be made to have a negligible effect on the expandor 
response. A unique feature of this particular configuration is that, re- 
gardless of the overall speed of the compandor, the time response of 
the feedback path can be matched to the feed-forward path by simply 
making RC, = ReC2; dynamic tracking of the compressor and ex- 
pandor is thus assured.* 

RC, and R2Cez are single-pole feedback networks of identical opera- 
tional amplifiers. The bandwidth of these elements, their time re- 
sponses, and the equivalence between bandwidth and time response 
may be analyzed in straightforward manner. The 3-dB bandwidth fre- 
quency of an RC low-pass characteristic is 


1 
I= 3 RG: 


* The representative spectra shown in Fig. 7 are illustrative of the significant 
bandwidths present at several points in the CNL terminal. Note that the noise 
bandwidth at the expandor variolosser is fe while the effective bandwidth of the 
cone signal paths is less because of the influcnce of RoCo in tandem with RiCi 
or Nee. 
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from which a direct connection between speed of the compandor and 
noise performance of the control channel will now be apparent. That 
is, as the speed of the compandor is increased (by decreasing RiC; 
and R.Cz), fe widens to admit more noise to the control lead of the 
expandor. The fluctuation of expandor loss about its nominal value 
produced by random noise is, in effect, low-frequency amplitude mod- 
ulation of the speech. This imparts a unique fluttering or “gritty” 
quality to the speech, which is a source of impairment over and above 
the effect of additive noise falling in the speech band. Interfering sig- 
nals falling in the control channel also produce loss deviations, but of 
a less random nature. A limitation on the noise reaching the expandor 
via the control channel imposes a corresponding limitation on fo = 
1/(27R2Cez). ReCo strongly influences the resulting speed of the com- 
pandor, and therefore the compandor noise improvement. Both the ef- 
fect of noise in the control channel and the compandor noise improve- 
ment are for the most part ultimately subjective. The task of finding 
an objective basis for evaluating these effects and analyzing them in 
detail is taken up in the following sections. 

In Section 3.2.2 and related Appendix B, control channel noise per- 
formance is analyzed, a criterion of performance is established, and 
subjective limits are given. 

In Section 3.2.3 the speed of the CNL compandor is discussed on a 
quantitative basis; relationships between time response and bandwidths 
are given. The related problem of dynamic tracking and the choice of 
operating ranges are dealt with in Section 3.2.4. 

In Sections 3.2.5 and 3.2.6, respectively, the expandor noise improve- 
ment and the efficacy of the complete compressor in loading the trans- 
mitter are discussed quantitatively, both as a function of compandor 
speed. 


3.2.2 Noztse Performance of the Control Channel 


The performance of the control channel has two aspects: (7) “normal” 
operation when the subcarrier is above a certain noise threshold and 
(11) the breaking threshold below which the channel breaks down. 

The breaking threshold is reached in the FM channel of the CNL 
system when noise or interference at the discriminator unit is so strong 
relative to the subcarrier that the polarity of the signal is too often 
reversed. Under these conditions, impulsive noise currents are gener- 
ated at the discriminator output; they drive the expandor to its ex- 
treme loss or gain. As the threshold is approached, the impulsive noise 
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currents punch “holes” in the speech or produce annoying intervals of 
excessive loudness. 

The threshold occurs with random noise when the subcarrier-to- 
noise ratio is in the vicinity of 10 dB on an rms basis. This ratio is a 
function of the line, or predetection, bandwidth and the subcarrier 
power. In the case of the experimental CNL system, the line band- 
width is 180 Hz, and the subcarrier power is 6 dB below the mean com- 
pressed specch power. Therefore, with the assumption that the noise 
is flat across the band, breaking occurs when the ratio of compressed 
speech to noise in the 2500-Hz speech band is about 5 dB. 

The bandwidth occupied by the FM sidebands must not exceed the 
180-Hz bandwidth allocated. It is a complex function of the amplitude 
and spectral distributions of the modulating signal and the frequency 
swing of the FM subcarrier. To further complicate matters, the modu- 
lating voltage is not symmetrically distributed and its spectrum is 
modified by the nonlinear LOG circuit preceding the frequency mod- 
ulator. It was experimentally determined that the FM spectrum re- 
sulting from a control signal rolloff at 100 Hz (8 dB down) and a 
peak-to-peak swing of 120 Hz was in accord with the 180-Hz line 
bandwidth. The 100-Hz roll-off occurs in the control rectifier filter, 
represented symbolically by RoCo on Fig. 7. The 120-Hz swing divided 
by the expandor loss range (60 dB, as discussed in Section 3.2.4) es- 
tablishes the deviation/loss ratio. It relates frequency errors of the 
FM subcarrier to the system loss error. These errors are due to noise 
and other mechanisms. The ratio enters into the calculations of Ap- 
pendix B. : 

The magnitude of system loss fluctuations above the breaking thresh- 
old due to noise in the control channel has been calculated by means 
of an expression derived in Appendix B. The derivation assumes that 
the noise is random and flat across the band common to both the com- 
pressed speech and the subcarrier. Although a more usual situation in 
HF reception is a mixture of noise and tone-like interfering signals 
having a variety of characteristics together with selective fading, the 
simple model nevertheless gives useful results. 

Equation (40) from Appendix B is 


loss fluctuation magnitude _ 3(f.)? E 1 (™)} 


(dB peak-to-peak) 35 s/NL 4,7 °° \y 


Fig. 8 is a plot of the loss fluctuation as a function of S/N with fe as 
a parameter. S/N is the ratio (numeric) of the speech voltage to the 
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noise voltage in the 2500-Hz speech band and is an indirect measure 
of the noise falling in the control channel. Expressing results in this 
manner makes it possible to correlate quality judgments based on the 
speech-to-noise ratio of the circuit with the degree of impairment 
caused by fluctuating system loss. By subjective tests, the approximate 
magnitudes were determined at which the loss fluctuations were judged 
(t) to have just noticeable impairment and (iz) to be so severe as to 
render the circuit uncommercial. These magnitudes are indicated on 
Fig. 8 as cross-hatched bands. 

The abcissa of Fig. 8 refers to the S/N at the input to the receiving 
terminal. In a circuit equipped with CNL-type terminals, a S/N of 
about 10 dB would be judged uncommercial and 15 dB would be judged 
poor but usable. By the standards normally applied to HF radiotele- 
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Fig. 8— Net loss fluctuation as a function of speech-to-noise ratio with control 
channel bandwidth as a parameter. 
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phone circuits, 25 dB would be considered good. Fig. 8 shows that if 
fe becomes large, the effect of noise in the control channel would be 
detectable even under otherwise good conditions. On the other hand, 
an fe of 5 Hz or less provides more margin than is necessary. With an 
fe of 20 Hz, when the S/N drops to 10 dB, the loss fluctuations reach 
an uncommercial magnitude. This same S/N would be judged uncom- 
mercial with a noise-free control channel. Thus, with an fo of 20 Hz, 
the limit of commercial quality would be reached more or less simul- 
taneously due to noise mixed with the speech and noise falling in the 
control channel. The value of f. adopted for the experimental CNL 
system was 15 Hz. It tends to be conservative in respect to loss fluc- 
tuation at the expense of compandor speed. The other speed-dependent 
performance factors discussed in Sections 3.2.5 and 3.2.6 change slowly 
with speed. 

The CNL system employs nonlinear signal processing, or “instanta- 
neous companding,’”’* ® within the control channel. This is the function 
of the LOG and ANTILOG circuits shown in Fig. 7. Such processing 
reduces the effect of noise on weak control signals. With weak talkers 
or during idle circuit conditions the control channel can tolerate more 
noise for a given requirement on expandor loss fluctuations if instan- 
taneous companding is used. While instantaneous companding increases 
the bandwidth of the analog control signal, a net improvement in con- 
trol channel performance results in the CNL application. The response 
of the LOG circuit itself is plotted on normalized linear scales on Fig. 
9. The ANTILOG characteristic must be sufficiently complementary 
to the LOG characteristic to insure overall linearity of the control 
channel. 


3.2.3 Time Constants 

Before considering CNL compandor speeds, it is necessary to have 
a precise definition of time constants. Also, an accepted definition will 
permit comparison of the CNL with other types of compandors, in 
particular, results obtained by other organizations (British GPO, etc.) 
currently working on CNL-type systems. The CCITT* has developed 
a definition for message compandor time constants that has found 
wide acceptance. This definition has the virtue of being readily imple- 
mented with simple measuring techniques and is suitable for specifying 
time constants in the CNL compandors described herein. 

The definition of attack and recovery time constants of a compandor 


* International Telegraph and Telephone Consultative Committee. 
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Fig. 9 — Log circuit input-output characteristic. 


involves a measurement of the transient response of a 2:1 compressor. 
In the CNL system, the first stage is identical to the conventional 
compressor and the second variolosser is slaved to the first. This en- 
ables the CNL compressor time constants to be measured using the first 
compressor in accordance with the methods specified by the CCITT. 
The attack time is measured by increasing the level of the input tone 
by 12 dB from —16 dBm0 to —4 dBm0 and measuring the time inter- 
val between the occurrence of the step and the time when the output 
envelope reaches 1.5 times the final steady-state value. The recovery 
time is the time interval between the occurrence of a downward step 
of 12 dB (from —4 dBm0 to —16 dBm0) and the time when the out- 
put waveform reaches 0.75 times the final steady-state value. 

For message compandors used on wire and radio relay circuits of 
moderate length, the recommended CCITT values are:1%14 


Attack time = 38 + 2 msec 
Recovery time = 13.5 + 6.5 msec. 


In order to relate measured attack and recovery times of a com- 
pandor to circuit elements as an aid to design, R. O. Carter’? ana- 
lyzed a simple model of a compressor (Fig. 10) and developed formulas 
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for calculating attack and recovery times. In his model, the control 
rectifier consisted of a peak detector (diode) and a single RC smooth- 
ing network. He showed that the CCITT attack and recovery times 
were given in terms of the product RC (which has the dimensions of 
time) as follows: 


Attack time = 0.15 RC (1) 
Recovery time = 0.675 RC. (2) 


The value of RC, since it is common to both the attack and recovery 
time constants, provides a single measure of the speed of the compres- 
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Fig. 10— Simple model of 2:1 compressor. 








sor. The bandwidth of the control signal into the variolosser is also a 
function of RC: 


f = 1/2 RC). * (3) 


Using relations (1) and (2), the CCITT recommended time constants 
of 3-msec attack and 13.5-msec recovery correspond to a speed of 
RC = 20 msec. The objective of the following analysis is to examine 
the CNL system model (Fig. 7) and find an RCottective for the some- 
what more complicated CNL compandor. 

The transient response for the compressor configuration with RoCo 
and RC, networks was calculated approximately by divorcing the 


* This relationship permits a ready transformation between factors best con- 
sidered on a time (or speed) basis and factors best considered on a frequency 
(or bandwidth) basis. It is exact only when a simple RC network is used in the 
model. The bandwidth, f, corresponds to the frequency at which the control 
path gain is 3 dB down from the value at de. 
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time-constant. determining components from the nonlinear feedback 
circuit and solving their response to a step input. A similar step re- 
sponse was determined for the simple RC network of Carter’s model 
for the time t = 0.150 RC, the CCITT attack time. The two responses 
were made equal and the attack time for the CNL compressor, f4, was 
obtained by solution of 


( 0.150 RC 
LS 6xps, Se 


RC ) = 0.139 


1 1 
RCo exp (zh ts) — R,C, exp (ai ts) 
Rc, — RG, - 


The use of the step response of the networks alone to obtain an engi- 
neering estimate of the attack time was justified because the input to 
the speed-determining element under CCITT test conditions is essen- 
tially a step function in the time scales considered here. The speed of 
the CNL compressor (given by an RCerrective) is determined from this 
calculated attack time through use of (1). Similarly, the recovery time 
was obtained by solution of the following equation for ftp: 





=l]+ 


exp (—2875 RC) _ 9 509 
a iC; exp ec RC Ron tn) - RCo exp (-< RC G, ts) 
Ri, = Bac, (5) 


This calculated recovery time can also be used to obtain RCerrective US- 
ing (2). The validity of expressions (4) and (5) was confirmed by 
direct measurement of attack and recovery times. Expressions (4) and 
(5) reduce to the form of (1) and (2) when either RC is much greater 
than the other. When the two time constants are equal, 


t,(for RC, — RoC) — 0.652 R,Co (6) 
tr (for FiGa ad RoCo) = 1.649 RCo ’ (7) 


an attack time increase of over 4 times and a recovery time increase 
of almost 24 times over the values of a single RoCo, rather than a 
twofold increase as might be surmised at first glance. Plots of both the 
attack and recovery equations (4) and (5) for values of RoCo and 
R,C, in the range of interest are given in Figs. 11(a) and 11(b). These 
curves illustrate how the values of RoCo can be adjusted to yield a 
desired attack or recovery time when R,C, is specified. 
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Mig. 11— (a) Calculated attack time; (b) calculated recovery time. 
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In Section 3.2.2, the frequency determined by ReC2 was given as 15 
Hz; substituting this value in (3) gives an R,C,; = ReCe of 10.6 msec. 
The objective in selecting RoCo was to yield compressor attack and 
recovery times recommended by the CCITT. This resulted in a choice 
of 


R,C,(attack) = 1.6 msec 
R,Cy(recovery) = 5.0 msec.* 


From Figs. 11(a) and 11(b) these values give overall attack and re- 
covery times as follows: 


Attack time = 3 msec 
Recovery time = 13.5 msec. 


The corresponding RCoftective given by (1) and (2) was 20 msec. 
Fig. 12 shows the relationship between overall speed (RCcttective) 
and the characteristics of the lowpass networks (RiC; and RCo). 


3.2.4 Dynamic Tracking and Operating Ranges 

As mentioned in Section 3.2.1, the control paths from control recti- 
ficr to variolossers pass through essentially identical smoothing net- 
works. Therefore, the input control signals to the transmitting com- 
pressor variolossers and receiving expandor variolossers have the same 
transient characteristic. Matching the control signals into the vario- 
lossers insures that the expandor loss matches the compressor gain on 
a dynamic basis. However, an imperfect reconstruction of the speech 
waveform at the system output can occur if: 


(t) the transmission time delay of the control signal differs from 
that of the corresponding speech signal, causing a mismatch between 
the instantaneous expandor loss and the envelope of the speech signal, 
and/or 

(it) the fade regulator is too fast, thus interpreting the residual 
amplitude variations in the output of the transmitting complete-com- 
pressor as fades and introducing further compression (without corre- 
sponding expansion). 


These effects result in speech that sounds distorted and possible cir- 
cuit instability if round-trip propagation delay is short. 


* The two values of RoCo require that the simple model of Fig. 7 be modified 
to provide a control rectifier with different charge and discharge times. 
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Fig. 12— Overall speed versus low pass network time constant. 


It has been found that a fade regulator RC time constant approxi- 
mately 4 times RCetrective reduces the additional compression distortion 
to acceptable levels without affecting performance during fading. 

Flat voice-frequency delay equalization in the speech or pilot channel 
is used to correct the transmission time differences due to filter and 
privacy flat delays. For speech service, it is desirable to reduce the 
time difference to no more than 5 to 10 msec. For signaling and similar 
pulse transmission where pulses are amplitude modulated on a signal- 
ing tone, the distortion due to time differences is more severe; a l- or 
2-msec difference produces noticeable changes in pulse shape. Equali- 
zation requirements for actual signaling systems of this type would 
depend upon system sensitivity to pulse shape. 

The need for time delay equalization is reduced if the CNL com- 
pandor is slow. If R,C,; = ReCe is large compared with the time dif- 
ference, the rate of change of control current into the variolossers is 
slow in any interval corresponding to the transmission time difference. 
The expandor loss mistracking errors will therefore be small since the 
slope of the loss curve as a function of time is low. 
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The operating ranges of the variable gain clements of the CNL sys- 
tem are interdependent since the sum of the variable gains and losses, 
including medium variations, results in an approximately constant 
value under both dynamic and static conditions. This is essential for 
system stability and tracking. The transmitting compressor range is 
40 dB; the fading regulator downward range is 20 dB and the upward 
range is 10 dB; the expandor range is 60 dB, and is the sum of the 
transmitting compressor range and the fading regulator downward 
range (sce Fig. 13). The operating ranges chosen were a reasonable 
compromise based upon expected speech volume and fading ranges, as 
well as practical considerations. 


3.2.5 Syllabic Compandor Quieting Versus Speed 


The speed of the compressor affects the degree of compandor quiet- 
ing of system noise heard by the subscriber. A limited subjective com- 
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Fig. 13 — Static operating ranges. 
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Fig. 14 — Distribution of S/N differences for equal subject preference. 


parison of a fast and slow CNL compandor (RC = 5 msec; RC = 75 
msec) was made to determine the extent of the change in noise quiet- 
ing solely related to speed. To this end, the FM control channel was 
made perfect (replaced by a de connection), noise was added to the 
speech path, and recordings were made. The slow system was taken 
as the reference system, and an approximate 10-dB compressed speech- 
to-noise ratio maintained. The fast system was recorded with a vari- 
able speech-to-noise ratio over a range of 6 to 12 dB. Direct compari- 
sons were made between the two sets of recordings by 12 test subjects. 
Each subject determined the value of variable noise (in the fast sys-— 
tem) which in his judgment made the fast and slow systems equivalent 
in quality. The reference value of 10 dB was chosen because it is the 
approximate point at which the CNL system quieting would be most 
effective on a circuit that was noisy but still usable. Fig. 14 shows the 
resulting distribution of S/N differences between fast and slow CNL 
compandors where the subjects indicated equal preference. For the 
median subject, the fast system gave a 24-dB improvement. Since the 
speed-related difference in quieting is strongly affected by the amount 
of noise present, the 24-dB advantage will decrease when the S/N 
improves from the 10-dB test value. 
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Fig. 15— Chart record of compressor input and output for fast and slow time 
constants; recorder integrating time—approximately 100 milliseconds. 


The conclusions drawn from the tests were that changes in com- 
pandor speed in the range of RC = 5 to 75 msec do not produce a 
sufficient change in the overall noise improvement to warrant consid- 
eration as a major factor in system “trade offs.” Since the total noise 
improvement for the average talker is 20 to 30 dB, the change in per- 
formance with speed is of minor consequence. 


3.2.6 Transmitter Loading 


The complete compressor serves to reduce the variations in the speech 
signal into the radio transmitter. However, residual variations do re- 
main and these variations become greater as the speed of the com- 
pandor is reduced. The regulated speech output from the transmitting 
compressor was examined to determine the extent of the regulation im- 
provement as a result of increasing the speed. The method used in this 
examination was based upon recognition of current CCIR recommen- 
dations. In these recommendations the mean power of the speech 
signal with smoothly read text is used as a criterion of ISB trans- 
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mitter loading. The mean power is defined: “The power supplied to 
the antenna transmission line by a transmitter during normal opera- 
tion, averaged over a time sufficiently long compared with the period 
of the lowest frequency encountered in the modulation. A time of 0.1 
sec during which the mean power is greatest will normally be selected.” 
For purposes of comparative measurements, the mean power was in- 
terpreted as that value obtained by a running measurement of the 
speech signal power with about a 0.1-see integrating time and by se- 
lection of the value of signal power when the running measurement is 
greatest, 1.e., at the crest of the speech envelope. 

CNL output signals were examined using a power measuring chart 
recorder with about 100-msec averaging time to determine how the 
mean power and the distribution of signal peaks changed as the speed 
of the transmitting compressor was adjusted to a very fast time con- 
stant (RC = 5 msec) and then to a very slow time constant (RC = 
75 msec). The input speech and the regulated speech output from the 
compressor were recorded on a strip chart (see Fig. 15) for each time 
constant extreme. The strip chart shows the power on a dB scale, 
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Fig. 16— Input and output speech envelope peak distribution. 
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where an output of 0 dB corresponds to the steady state output of the 
compressor with an input 0 dBm0 1000-Hz tone. 

When speech is present, the regulated output signals with both fast 
and slow time constants are no more than 5 dB apart. More impor- 
tantly, the slow compressor (with more overshoots expected due to 
more sluggish regulation) has a mean power that is at most 2 dB 
higher than the fast compressor. Following accepted practice, the level 
into the transmitter for RC = 75 msec would then have to be dropped 
only 2 dB compared to the level for RC = 5 msec. If the peak values 
for each output are plotted to give a probability distribution (see Fig. 
16), the slight difference between the outputs as a function of com- 
pressor speed is evident. In Fig. 16 the spread of the speech envelope 
peak distribution of the input speech is 26 dB. This spread of 26 dB 
is reduced to 9 dB at the output of the slow compressor and to 6 dB 
at the output of the fast compressor. 

Based on the above results, it was concluded that any otherwise 
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Fig. 17 — Chart record of compressor input and output—RC = 20 milliseconds; 
recorder integrating time—approximately 100 milliseconds. 
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Fig. 18 — Chart record of compressor input and output—RC = 20 milliseconds; 
recorder integrating time—approximately 10 milliseconds. 


proper RC time constant in the range 5 to 75 msec would produce 
nearly equivalent transmitter loading. The value of RCcéfective = 20 
msec used in the trial equipment falls within the above limits and, in 
fact, produces a load that has almost as little variation as the fast 
RC = 5 msec compandor. Fig. 17 shows an example of the input and 
output waveforms for the RCertective = 20-msec system. 

To verify that the accepted 100-msec averaging time used in this 
investigation has not obscured any very short term effects, an addi- 
tional chart recording was made of system input and output with a 
10-msec integrating time. This record, a portion of which is reproduced 
in Fig. 18, shows that while some overshoots remain, their magnitude 
and duration are not expected to produce unwanted effects in the radio 
transmitter. 
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APPENDIX A 


Gain Relations in Single- and Two-Stage Compandors 


Application of conventional compandors™ to conventional wire, ca- 
ble, and radio relay circuits produces level characteristics such as those 
shown in Fig. 19. At the compressor, low-level inputs are amplified 
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Fig. 19 Example of 2:1 compandor action. 


more than high-level inputs; correspondingly, at the expandor, low- 
level outputs result from more attenuation than that given to high- 
level outputs. (Note: certain very high inputs are attenuated to reduce 
overload.) As shown in the example of Fig. 19, conventional syllabic 
compandors halve the volume range between compressor and expan- 
dor; thus, a 0 to —40-dB range becomes a 0 to —20-dB range in the 
transmission portion. These compandors are known as 2:1 compan- 
dors;* the 2:1 characteristic is easily obtained as will be shown below. 

The compressor portion of the compandor consists of only four sim- 
ple elements (Fig. 20). The variolosser is a circuit configuration where 


*In CCITT documents, this is stated as a compression ratio of 2. 
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the attenuation ratio in the speech path is inversely proportional to 
the unidirectional control current. The variolosser equation is 


C1 =h =, (8) 


where é€i, is the amplitude of the input envelope. The control rectifier 
acts as an envelope detector, producing control signals with a spectrum 
extending from 0 Hz to approximately 100 Hz.* The feedback path 
directs the control signal to control the variolosser attenuation. The 
2:1 compression effect results because the control signal developed 
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CONTROL 
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Fig. 20 — Elements of 2:1 compressor. 


from the output is fed back to control a gain adjusting variolosser, 
which in turn changes the output. The fixed gain amplifier is added to 
the speech path for level adjusting purposes. The amplitude of the 
output envelope is 


Cour = Keer . (9) 
The control current 1s related to the output envelope by 
12S heey (10) 
from which 
Cout = Feet in (11) 
Ve 
—_ kyketin 
Cout . ks€out (12) 
Cout = ks Cin . (13) 


* The upper cutoff frequency is not a clearly defined point, depending on the 
characteristics of the speech signal and the type of rectifier smoothing. 
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Taking logarithms 


20 log eu, = 10 loge;, + 10 log ae (14) 
3 

10 log éou, = (10 loge,,) + K (15) 

E,4.(dB) = 3#,,(dB) + K’. (16) 


The output volume changes 1 dB for every 2-dB change in input vol- 
ume. The foregoing analysis is based upon steady state conditions. The 
transient response (or time constant) of the compressor is determined 
by an RC network in the control rectifier. The transient response of 
compandors and the CNL system is discussed in the analysis section 
of the body of the paper. 

The expandor portion of the compandor consists of five simple ele- 
ments (Fig. 21). In this case, the variolosser is a circuit configuration 
where the attenuation ratio in the speech path is directly proportional 
to the unidirectional control current. The control rectifier acts as an 
envelope detector in the same way as in the compressor, but in this 
case, the control signal is applied in a forward acting manner rather 
than in a feedback manner. The variolosser equations 


€, = ketein (17) 


Cout = ksey = kskete€in (18) 


are then manipulated as in the case of the compressor and the follow- 
ing result is obtained: 


Cout = legkakskeein (19) 
E4(dB) = 2K,,(B) + K’’, (20) 
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Fig. 21— Elements of 2:1 expandor. 
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The output of the expandor changes 2 dB for every dB change in the 
input, thus undoing the effect of the compressor. 

Idealized input-output characteristics of the conventional compandor 
are shown in Fig. 22. When installed at each end of a circuit, the out- 
put will be a true replica of the input if the circuit is distortionless. 
Since the expandor control signal is derived from the output of the 
transmission path at the input to the expandor, examination of the 
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Fig. 22 — Input/output characteristic of 2:1 compandor. 


characteristics of Fig. 22 verifies that gain variations (with time) of 
the medium will be magnified in dB by a factor of two. For this rea- 
son, conventional compandors have not found application in HF radio 
service.* 

The CNL compandor arrangement (Fig. 23) solves the problem of 
net loss variations while providing essentially constant output from 
the transmitting compressor. Constant output is obtained by adding 
a second slave variolosser controlled by a replica of the regular con- 


* Compandors were used for the first transatlantic radiotelephone circuit® (long 
wave) where the net loss was relatively stable. 
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Fig. 23 — Example of CNL system gains and losses. 


trol signal present in the feedback path of the normal first stage com- 
pressor. Examination of Fig. 22 shows that for any input in the range 
of compression, a conventional compressor provides a variolosser con- 
trol signal so as to insert an amount of gain equal to one-half the dif- 
ference (in dB) between the “unaffected level” (or input level at which 
no gain or loss occurs in the compressor or expandor) and the input 
signal level. The slave variolosser gain is also controlled to be one-half 
the difference between the unaffected level and the input signal level, 
so the sum of the gains of two stages is sufficient to raise the level of 
the input signal to the unaffected level for any value of input signal 
in the range of compression. An example showing the CNL two-stage 
compressor gains and levels for two inputs is given in Fig. 23. 

At the receiving terminal, expandors are used to restore the speech 
distribution. Two normal expandor variolossers are used, since each 
can accomplish the inverse of the action of one of the compressor 
variolossers. However, the local control arrangements at the expandors 
can not obtain control information from the input signal (which is 
nominally constant for all speech in the absence of fading). Re- 
examining the compressor and expandor equations 


. out kike 
compressor gain = ““ = (21) 


Cin te 





expandor gain — et = kskete . (22) 
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Thus, if a compressor and expandor are in tandem, their net gain is 
constant if 

t, =%,. (23) 
A suitable separate control channel, which supplies each expandor with 
an 1, that is a replica of 2, , assures this identity. 

To maintain a constant net loss when there is fading in the trans- 
mission path, the CNL system uses a fading regulator that is similar 
to the transmitting compressor, except that the fading regulator has 
no corresponding expandor. The fading regulator action is based on the 
recognition of the requirement that the output of the radio receiver be 
constant since the input to the radio transmitter is controlled to a 
constant value by the transmitting compressor. If the output of the 
radio receiver is less than the nominally constant value, it is assumed 
that a fade has occurred. 

The fading regulator inserts sufficient gain to raise the output to 
the nominal nonfade value required by the expandor. Fig. 23 shows 
the gains and losses with a test tone input for an example without fad- 
ing and for an example with an assumed flat fade in the medium. In 
this simple model the overall effect is to provide a nominally constant 
overall net loss even with a time varying transmission medium. 


APPENDIX B 


Net Loss Fluctuation Caused by Noise in Control Channel 


The effect of noise in the control channel can be obtained by use 
of the model shown in Fig. 24 (not to scale). A flat band of random noise 
extends across the speech and subcarrier bands. The bandwidth of the 
noise affecting the speech directly, f, — f, , is 2500 Hz. The bandwidth 
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Fig. 24 — Signal and noise spectra. 
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B is the line, or predetection bandwidth, of the control channel, which 
is assumed to be defined by an ideal rectangular filter. It is assumed that 
the peak noise voltage in this band is always less than the peak sub- 
carrier voltage A, . 

In the following analysis, an expression is derived that relates the 
magnitude of fluctuations in the speech channel loss, to the noise level. 
The noise level enters this expression implicity in terms of S/N, the 
speech-to-noise voltage ratio in a 2500-Hz band. It is thus possible 
to relate judgments concerning S/N and loss fluctuations, respectively. 

The flat band of noise of Fig. 24 can be represented by an infinite 
number of equal-amplitude sinusoids having incommensurable but 
approximately uniformly spaced frequencies and incoherent phases. 
The noise in an incremental bandwidth df is represented by a single 
sinusoid whose mean-square voltage is equal to that of the noise. The 
mean-square noise voltage in df is just a;df, where a, is the noise voltage 
density constant in rms volts per unit square root bandwidth. Thus, 
if A, denotes the peak amplitude of the noise component, 

AS 2 
os = a, df. (24) 

The superposition of the subcarrier and a noise component at a 
frequency f, = w,/27 relative to the subcarrier can be represented by 
the phasor diagram (Fig. 25). The peaks of the subcarrier and noise 





Fig. 25 — Phasor diagram, carrier and noise. 
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components are denoted by A, and A,,, respectively. It is assumed 
that for A, «< A, the resultant angle modulation of the subcarrier 
produced by the noise will be equal to the summation of the modula- 
tions that would be produced by each noise component separately. 
Except for noise-produced modulation, it will be assumed that the 
subcarrier is otherwise unmodulated and centered in the band. Am- 
plitude variations may be ignored on the basis that the subcarrier 
frequency demodulator is ideal. 

Rotation of the phasor A, about A, causes the resultant FR to oscillate 
about A, with a peak phase deviation 








Ag = tan” An nw Ar 


A AS (25) 


It can be shown’® that if A, « A, as assumed, ¢(t) is essentially 
sinusoidal and equivalent to pure phase modulation of the form 
M(t) = A, sin @,¢ + Agsin pi + y), 


where the modulating signal is a sinusoid of arbitrary radian frequency p. 
By differentiation of arg M(t), the instantaneous frequency is obtained: 


oO; = © (wut + Agsin pt + y) 
(26) 
=, + p(Ag) cos pi. 


The instantaneous frequency deviation produced by the noise component 
at the frequency w, relative to the subcarrier can be deduced from (26) 
by using (25) and setting p = a, . 


Inst. freq. dev. = w, 4 COS wf, radians per sec. (27) 
The mean-square frequency deviation due to a single noise component 


is thus 


((freq. dev.)*),, = 1( : 4] ; (radians/sec)” 


(Se). ae 


To obtain the incremental mean-square frequency deviation (dé”),, due 
to noise in the band df at a frequency f, relative to the carrier, (24) 
is substituted into (28): 


(d82)4. = (sels) dj, Hy’. (29) 


I 


718 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1967 


If a discriminator sensitivity of 1 volt per Hz is assumed for convenience, 
the integral of (dé”),, over the line bandwidth may be equated to the 
positive-frequency discriminator output mean-square noise voltage 
spectrum integral: 


fet+B/2 B/2 
[am =f a at, (30) 
Se—B/2 0 
where 
_(!), vam, osssBp 
0, Lo B/2: 
The noise spectrum at the output of the single-pole post detection filter is 
FQ) =| YA |G),  V/Haz, (32) 
where 
| YO |= + G/p)T? (33) 


is the voltage transfer function of the filter and f, = 1/27R.C, is the 
frequency at which the response is down 3 dB. Substituting (31) and 
(33) into (82) and integrating over half the line bandwidth gives the 
total mean-square noise voltage output of the control channel that acts 
on the expandor variolosser: 


ie FAA) df = o( 2) [- Read dj, Vv’. (34) 


The effect of the ANTILOG circuit (see Fig. 7) on the noise output 
of the control channel involves a nonlinear transformation of a random 
process that is difficult to handle mathematically, hence the ANTILOG 
circuit has thus far been ignored. In spite of this omission, the approxi- 
mate results derived herein were in good agreement with measurements. 

It is to be noted, however, that because of the ANTILOG circuit 
the translation from a frequency deviation at the FM receiver input 
to a change in expandor variolosser loss in dB is linear with a slope 
(or deviation/loss ratio) of 2 Hz per dB. By referring the control channel 
noise output voltage, (34), back to the discriminator input, an equivalent 
frequency deviation can be obtained; this deviation, divided by the 
deviation/loss ratio, gives the fluctuation in expandor loss directly in 
dB. Accordingly, when this mean square deviation, denoted by 6”), , 
is equated to the right side of (84) (on the basis of the assumed dis- 
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criminator sensitivity of 1 volt per Hz) and the integral is evalu: 


we have: 
Oa = (2 a) Gay |2 OF. — tan” (2) |. Hz’. (‘ 


In the CNL system, B is 180 Hz. Substituting this value into (35) a: 
taking the square root to get rms, the result is 


“sf oy 
2 tan : Hz. 36 
aia | 20 — tan (36) 

This result is next recast in terms of S/N, the ratio of the rms com- 
pressed speech voltage to the rms noise voltage in the 2500-Hz speech 
band. First, note that 








Dias _ 


c 


V2 


since the subcarrier is transmitted 6 dB below the speech in the CNL 
system; note also that the rms noise voltage 





bo |~R 


(37) 


= rms subcarrier voltage = 





2500 h 
W-[f ee 
: (38) 
= 50a, . 
Then substitution of (37) and (88) into (36) yields 
Rue Ee a (22) } 
c= 2 S/N LF tan ‘1? Hz. (39) 


Dividing (39) by the deviation/loss ratio, 2 Hz/dB gives the expandor 
loss error in dB rms. The result can be expressed in dB peak-to-peak 
(by ignoring peaks in excess of 3 times rms) if (89) is multiplied by 6. 
When (39) is modified by these two factors, the result is 


Loss error = (4)(6) 6::s (40) 
= 3(f2)" E = (0 ae eak-to-peak. 
= 35SN 7, tan f, dB peak-to-pea 


Fig. 8 is a plot of the loss error versus 20 log (S/N) with f, as a parameter. 
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Experiments with an Underground 


Lens Waveguide 


By D. GLOGE 


(Manuscript received December 21, 1966) 


A laser beam was transmitted over a distance of 1/2 milein an underground 
tron pipe using glass lenses 400 feet apart. The beam deflection and the 
temperature gradient in the air-filled pipe was measured simultaneously 
over periods of several weeks. It was found that gradients up to 0.02°C/cm, 
changing with the season, accounted for beam displacements up to 2 cm on 
the lenses. Since all these deviations are slow and no other severe disturbance 
was noticed, the conclusion is that the transmission of optical beams under- 
ground is possible without evacuation of the conduit. Suggestions for the 
construction of such a beam waveguide are made which could reduce the 
mentioned temperature influence by a factor of 100. 


I. INTRODUCTION 


If a beam of coherent light is to be transmitted along the surface 
of the earth it will be necessary to redirect and focus it at intervals 
by means of lenses or mirrors to follow the terrain.’’? Furthermore, 
it must be shielded from atmospheric temperature fluctuations which 
result in variations of the index of refraction. A laboratory experiment 
which, in a 100-m metal pipe, folded a transmission path back upon 
itself gave the impression that the atmospheric effects could be overcome 
by choosing the proper beam enclosure.*’* Other experimenters, using 
a 1-km long pipe above ground, considered evacuation as the most 
reasonable means to avoid temperature effects.° In any case, the trans- 
mission path would most likely be installed underground where temper- 
ature variations are much smaller than in the open air. 

To gain information about the transmission characteristics in this 
case a 3-mile underground iron pipe was used to build a lens waveguide 
with focusing lenses 400 feet apart. The temperature field and the 
beam displacement in the tube were measured simultaneously and 
compared with theoretical estimates. This gave valuable indications 
for the construction of underground lens waveguides. 
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II. GROUND TEMPERATURE 


The temperature distribution in the ground depends very much on 
the type of soil, on the moisture content, and on the character of the 
surface. The common feature of any soil, however, is a high heat capacity 
c and comparatively low conductivity « which leads to very slow heat 
diffusion processes. To give a representative picture, a half space with 
constant ¢ and x is considered. Then heat diffuses only along the vertical 
direction y and the temperature 7’ can be found from the equation 


oT oT 
ot =a ay” ’ (1) 
where the ratio 
K ¢ 
a=: (2) 


is called the thermal diffusivity. 

Since the temperature changes at the surface are mainly periodic 
it is convenient to consider the Fourier amplitudes 7, of components 
with frequency w. The solution of (1) in complex form is then 


Tyo) = T., exp (twt — yV w/a). (3) 


The temperature is retarded and attenuated with increasing depth 
and increasing frequency. 

To find the diffusivity a involved in this experiment, thermocouples 
were installed 0, 1, and 3 feet below the surface. A recording over 2 
weeks is shown in Figs. 1(a) and (b). A comparison of the temperatures 
at the surface and 1 foot below show that a 1-foot layer of soil reduces 
the amplitude of the diurnal temperature cycle by a factor of 10 and 
introduces a phase lag of 10 hours. This corresponds to a diffusivity 
of a = 5.5 10°° em’/s. 

Data gathered at different locations in the United States yield an 
average diffusivity of 5.4 107° em’/s.* Handbook data for clay soil 
are shown in Table I. According to these data the diurnal temperature 
cycle at a 3-foot depth is reduced by three orders of magnitude and 
therefore not noticeable in Fig. 1(b). 

The recording of Fig. 1(c) shows the temperature gradient in the 
iron pipe 5 feet below the surface, which was used for the transmission 
experiment. This gradient was measured by a digital quartz thermometer 
whose sensors were installed at the top and the bottom of the tube. 
These sensors are connected with the equipment by 100 feet of cable 
so that the measurements could be taken far inside the tube and remote 
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Fig. 1— Temperature recordings at Bell Telephone Laboratories, Holmdel, 
N. J. from April 22 to May 4, 1966. (a) Temperature at the surface. (b) Tem- 
perature at a depth of 1 and 3 ft. (c) Thermal gradient in an iron pipe with a 
wall thickness of 3” and an i.d. of 3.5” at a depth of 5 ft. 


from any disturbance. No diurnal temperature changes were noticeable. 
From (8) one finds the temperature gradient in the soil to be 


oO, = sy = —Paaft= exp (twt — yV tw/a). (4) 


Gradients of very low and very high frequency disappear at a large 
enough depth. If the Fourier spectrum of the temperature on the surface 
is assumed to be flat then Fig. 2 gives the amplitude spectrum of the 
gradient at various depths. 

















TABLE I 
Still Clay Alumi- Fiber- 
air soil Steel num glass 
Conductivity « cal__ | 5.510-8] 2.810] 0.11 | 0.55 10-1 
em s°C 
Heat capacity c oo ON. 3107 | 0.5 1 0.65 — 
em? °C , 
: Les cm? : 
Diffusivity @ = 0.18 5.410-*/ 0.1 0.85 = 
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Fig. 2— Amplitude spectrum of the thermal gradient at various depths related 
to the amplitude at the surface of the ground. 


Actually, the temperature spectrum at the surface is not constant 
but has peaks with periods of 1 day and 1 year. About 1 year is also 
the period which causes maximum response of the temperature gradient 
at a depth of 5 feet. 

The mean temperature behavior valid for the area of New Jersey, 
as shown in Fig. 3, is based on data gathered in Ref. 6. The corresponding 
temperature cycle at a depth of 5 feet was calculated from (2). The 
dots show measured values. The presence of an air-filled iron pipe distorts 
the temperature field considerably. As is shown in the appendix, a 
function describing this field can be derived from an equivalent problem 
in the theory of electrostatic fields. Assuming that ?, , the gradient in 
the absence of a tube, is approximately constant over an area comparable 
to the tube crosssection, the field inside the tube is found to be homoge- 
neous having a thermal gradient 3; that is given by (18). 

For the tube used the inner radius r, = 1.75 inches, the outer radius 
ro = 2 inches and x, x, and x, are the thermal conductivities of air, 
iron, and soil, respectively (see Fig. 4). With the numbers given in 
Table I, one finds that the internal gradient 3; is reduced by a factor 
of 0.4 in comparison with #3, . Therefore, 


3, = 0.44 ee exp (iwt — yViw/a). (5) 
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Telephone Laboratories, Holmdel, N. J. in 1966. 





Fig. 4— Sketch of the pipe arrangement used to derive the internal tempera- 
ture field. 
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As long as the thermal gradient has no horizontal component and 
its value is below a certain threshold, convection cannot develop in 
a cylindrical container.’ For the air-filled tube with the given dimensions, 
this threshold is 0.02°C/cm. 


Wi. THE LENS WAVEGUIDE 


The experimental setup as shown in Fig. 5(a) makes use of a straight 
iron pipe of 4-mile length with an inner diameter of 3.5 inches at a 
mean depth of 5 feet. Roughly every 140 m there is access to the pipe 
through manholes. Here biconvex thin lenses with a focal length of 
70 m are installed. 

The fundamental mode of this confocal lens waveguide is a beam 
with Gaussian field distribution. Operating at the laser wavelength 
of 0.63 yu, the 1/e width of the field at the lenses is 21/140 m 0.63 p/m = 
10.4 mm. The lenses are supported by the pipe itself and mounted 
close to the wall of the manhole where the pipe comes out of the ground. 
No vibrations were noticed. The conduit is airtight over the whole 
transmission length. 

The lenses have a diameter of 60 mm and can be adjusted for about 
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Fig. 5— Sketch of the experimental setup: (a) The transmission path. (b) The 
launching system. 
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(a) (b) 


Fig. 6— Possible redirecting devices in a lens waveguide. (a) Mirror arrange- 
ment as proposed by Kompfner.’ (b) Rotatable split lens which contains a vari- 
able prism. 


+15 mm from the outside of the pipe. This is sufficient to guide the 
beam in the conduit which is straight over the total length of 4 mile. 
Of course, aligned sections cannot be expected in an actual field instal- 
lation. Redirecting devices must therefore be provided where the 
described model has only focusing lenses. 

In order to gain a feeling for the usefulness of the described model, 
some thought will be devoted to the features of redirecting devices. 
Fig. 6 shows two examples. In Fig. 6(a) a plain and a curved mirror 
are mounted opposite to each other.’ Fig. 6(b) shows two eccentric 
parts of semiconvex lenses which can be rotated independently. The 
space between the adjacent surfaces is filled with an oil film which 
has the same refractive index as the lenses. The complete device acts 
therefore as a lens plus a variable prism. 

The device in Fig. 6(a), as well as the one in Fig. 6(b), uses two 
(but not more than two) reflecting or refracting surfaces both to direct 
and to focus the beam. Those surfaces determine essentially the trans- 
mission loss in an optical waveguide.*’*’’ Therefore, redirection of the 
beam may be achieved without a large increase of the transmission 
loss, no matter if reflection or refraction is used. The experiment 
described here was started with refracting devices since they are easier 
to mount and to adjust, but an experiment with reflecting devices 
is planned as well. 

The lenses used in the experiment have a surface quality of \/10 
and are furnished with an antireflection coating which guarantees less 
than 0.2 percent reflection at 0.63 » wavelength. The bandwidth is 
about 500 A. The absorption loss of the lens material is in the order 
of 0.02 percent for a lens thickness of 5 mm. Assuming a total absorption 
loss of 0.1 percent in the glass and the layers of the two coatings, an 
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overall loss of 0.5 percent must be expected per lens. (This varies by 
about +0.05 percent since the reflections from both surfaces interfere 
with each other and the total reflection therefore depends on the lens 
thickness. ) 

A loss of the same order for mirrors was reported in Ref. 3. No 
significant variations of the loss are therefore to be expected in this 
case. For the mean lens spacing of 140 m a waveguide loss of 0.15 dB/km 
would occur. This, of course, holds only for perfectly clean surfaces. 
A considerable increase of the loss must be expected if dust accumulates 
on the lenses. No loss measurements have been performed up to the 
present time. 

The Gaussian field distribution is generated in the almost concentric 
resonator of a He-Ne laser and magnified by a factor of 8 in the lens 
arrangement shown in Fig. 5(b). These lenses are polished to a quality 
of 4/10 but are not corrected for aberrations. 

The launching system consisting of laser and magnifier as shown in 
Tig. 5(b) is mounted on a solid concrete table resting on a 1-foot thick 
concrete floor. Although the air temperature around the launching 
system is kept constant within 2°C the temperature changes in the 
environment are large enough to cause deviations of the proper launching 
angle of 1.5 107° rad or 2 mm at the first lens. To climinate this effect 
a servo-loop was built which holds the beam center at the first lens 
within an area of 30 X 30 uw. Of course no temperature effect can now 
be measured in the first pipe section, but by evacuating this section 
it is possible to attain a deviation angle of the forthgoing beam of no 
more than 30 4/128 m = 2.3 1077 rad. 

To measure the displacement of the beam the circuit in Fig. 7(a) 
was used. If the shown CdS8-cells have the resistances R, and R, the 
output voltage is 

V, R, ia R, 
Po Be Ry 6) 
Since these resistances, in the region used, are inversely proportional 
to the illuminating intensity J, (6) can also be written in the form 


Ve Veh. 
2 I, + J; 


This shows that the output voltage does not depend on the absolute 
laser intensity but only on the ratio of the intensities absorbed by the 
resistors, which is a function of the displacement d and the beam width. 
A typical sensitivity curve of the device is shown in Fig. 7(b). The 











(7) 
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Fig. 7— Device to measure the beam displacement with photoresistors. (a) 
The circuit. (b) The sensitivity curve. 


same circuits were used to produce the horizontal and the vertical 
feedback signals and to measure the beam movement at the end of 
the transmission path. 


IV. BEAM DISPLACEMENT 


Assume only small deviation of the beam from the ideal propagation 
axis z. Then the displacement d in a transverse direction y can be 
calculated from the equation 

a’ _ Lan, 
Oz” 


= a ay (8) 


Here n is the refractive index which for air at atmospheric pressure 
has the temperature dependence 


n= 1+ dn7 (9) 


with An = 3 10™*. This yields 


ad _ _An Lar 


azo T., oy ; (10) 


If 07 /dy is assumed constant over a transmission length L a simple 
integration of (10) gives the displacement 


_ An It af 


ae ie oy aD 
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at the end of the path Z, provided the beam is injected on axis at the 
beginning. To measure this effect the stabilized beam coming out of 
the evacuated pipe section was sent through the next section and its 
displacement recorded in the second manhole. At the same time the 
temperature gradient was measured in another pipe section. Inserting 
these values into (11) the expected beam displacement was calculated. 
The thermal gradient and the measured and calculated displacement 
are shown in Fig. 8. There are short-term movements superimposed 
on the main slope which are caused by diurnal temperature changes 
in the manholes. Here the pipe warms up over the day and cools down 
at night. This introduces convection currents along the pipe which 
cause a transverse temperature gradient and a deflection of the beam.’ 

Recordings of the beam position over the full transmission length 
again showed diurnal deviations (up to 4 mm) introduced by the 
temperature change in the manholes. 

So far only changes of the temperature gradient have been considered 
but nothing has been said about the absolute gradient. Fig. 3 shows 
that a maximum gradient of +0.02°C/cm must be expected. In this 
case, if the beam enters the waveguide on axis, the displacement on 
the second lens according to (11) is 2 cm. The third lens is passed at, 
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Fig. 8— Measured and expected beam displacement as calculated from the 
thermal gradient shown below for a laser beam transmitted over 128 m through 
an iron pipe at a depth of 5 feet in the ground. 
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Fig. 9— Beam path in the experimental lens waveguide at the time of maxi- 
mum temperature gradient. 


the center, but at the fourth lens again 2-cm displacement occurs. 
The beam undulates about the dotted line shown in Fig. 9. A proper 
launching of the beam would avoid the undulations and the beam 
would follow the dotted line. This behavior was indeed found quanti- 
tatively during measurements in June when a high temperature gradient 
occurred. 

Since the beam center has to pass a lens at least 1 cm from the lens 
edge to avoid noticeable diffraction, the described lens waveguide can 
just handle the occurring temperature effect. Of course, this waveguide 
runs below a uniform grass surface and the temperature gradient is 
fairly constant everywhere. An actual transmission path, however, 
will be subject to changes of the temperature gradient from section to 
section and that will cause undulations of the beam which increase 
with the transmission length. For this reason a thermal gradient as 
great as that measured would be a serious drawback for the waveguide 
and a considerable reduction of the temperature influence in the conduit 
seems necessary. 

This can be accomplished by employing a pipe material of higher 
conductivity. The reduction factor 3;/3, calculated from (18) is shown 
in Fig. 10 for iron and aluminum tubes with various diameter ratios. 
Using an aluminum pipe instead of the present iron pipe would reduce 
the inside thermal gradient by a factor of 5. 

The shielding effect can be improved even more by using several 
coaxial layers of strongly alternating conductivity. As a practical 
example for a triple-layer pipe, the combination aluminum-concrete- 
aluminum was calculated using a concrete volume three times larger 
than the aluminum volume. The thickness of the layers is chosen in 
such a way that the shielding effect is largest. As Fig. 10 shows, such 
a pipe can be more effective than a solid aluminum pipe if the layers 
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Fig. 10— The reduction factor versus the ratio of inside to outside diameter 
for a solid iron pipe, a solid aluminum pipe and aluminum pipes stacked with 
concrete and fiberglass using a volume ratio 1:3. 


are very thick. For the same amount of aluminum as in the solid tube, 
however, it is only a factor of } better. 

Using fiberglass instead of concrete improves the shielding effect 
considerably. A pipe with the same aluminum weight as before yields 
a reduction factor of 3 107°. Compared with the iron water pipe now 
in use, the shielding effect of the latter arrangement is two orders of 
magnitude better, and the maximum expected gradient in such 
pipe would therefore be 2 10°*°C/cm. 
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V. CONCLUSIONS 


The results of the described experiment show that the main dis- 
turbance of a coherent light beam transmitted in an air-filled under- 
ground conduit is a slowly varying displacement of the beam center. 
This displacement can be coordinated with thermal gradients in the 
atmosphere of the pipe and predicted theoretically. No other time 
dependent sources of disturbance were noticed in the experiment. The 
mounting, as well as the adjustment of the lenses, was not critical. 

To reduce the beam displacement to a tolerable amount some simple 
provisions must be made: (z) the conduit should be imbedded in the soil 
everywhere and exposure to temperature changes should be minimized; 
(it) the conduit should be airtight to avoid longitudinal air currents; 
(iat) a tube of high thermal conductivity or a multilayer tube with 
a good shielding effect should be chosen. With the proper design no 
larger gradients than 0.001°C/cem would occur in the pipe. If, for 
example, the gradient is not constant over the transmission path but 
varies with a rms value of 0.001°C/cm then, if the lens spacing is 
140 m and the beam enters on axis, beam undulations occur with a 
mean amplitude of 0.5 mm at the beginning of the waveguide, increasing 
to 5 mm after 100 lens sections. There the beam could be realigned by 
an electronic servo system.* 

Drift and vibrations of the ground were not investigated until now, 
but these might very well cause larger displacements than the ones 
just mentioned. This would necessitate beam aligning servo systems 
at even shorter intervals. 
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APPENDIX 


To find the temperature field in an around several coaxial cylinders 
of different material as shown in Fig. 4, the equivalent problem in 
the theory of electrostatic fields is considered. It is known that a circular 
dielectric rod with its axis perpendicular to a homogeneous electric 
field develops a homogeneous internal field determined by its polarization 
and an external field that consists of the superposition of the original 
field plus a field generated by the polarization of the rod. 
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In polar coordinates, a homogeneous temperature field with gradient 
3, has the form #,r cos ¢ (see Fig. 4). If now a circular rod is brought 
into the field, an additional term of the form (6/r’) cos ¢ arises originated 
by the “polarization” @ of the rod. In a more complicated configuration of 
several coaxial cylinders, the field in the various materials can still 
be described by terms of the above type. 

Labeling the n layers starting with » = 0 for the center rod, one 
has for the field 7’, in the rth layer 


T, = 2dr cose + “5 COs ¢ (12) 


for 
y=O0,1---n. 


7’, is the field in the surrounding and #, = #8, , the original thermal 
gradient. 6) = 0 since the field in the center bore is homogeneous. 
This leaves 2n constants 3, and 6,. They can be found by matching 
the fields at the boundaries which yields 2n linear algebraic equations 
of the form 





6, 6, 
i = =v, —- r (13) 
and 
Oye 6, 
cai(o,- + Bs) = ls, + 4), (14) 


To calculate the internal gradient 3; = ¢@, it is convenient to define 
the ratios 











+ y= + K, 
pi = Sate (15) 
pS (16) 
and a determinant 
+ rT) oe, n TY = 
hil oe Ef 7: 
Po ro Pi ( ) r, Pp 
r _ + n— To = 
—> Po Pi —(—1) ‘Ds 
D, = : ‘ - (17) 
nf es n-1 To + 
(=1) 7 Po (—1) Bi Dn 
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Then 


bo _ 1. 

o, OD, 

This is the ratio by which the gradient inside a multilayer tube is 

reduced in comparison with the gradient in the homogeneous soil 

(in absence of that pipe). For a solid tube as shown in Fig. 4, n = 2 
and one finds 

Lee a : (18) 

Be ko + ky ky + ky (*:) Ko Ki Ki — Ke 

2k, 


2k, 2k pe 
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CONS 


Axis-Crossings of the Phase of 
Sine Wave Plus Noise 


By A. J. RAINAL 
(Manuscript received November 16, 1966) 


This paper 1s concerned with the axis-crossings of the resultant phase 
—a S 6(t,a) S x of a sinusoidal signal of amplitude 2a and frequency fo 
plus Gaussian noise of unit variance having a narrow-band power spectral 
density which is symmetrical about fo. The discontinuous phase process 
O(t,a) is present at the output of the IF amplifier of a radio or radar receiver 
during the reception of a sinusoidal signal immersed in Gaussian noise. 
Also, the phase process 6(t,a) is basic in Rice’s recent analysis of noise 
in FM receivers. The following theoretical results are presented concerning 
the axis-crossings (level-crossings) of 0(t,a) at an arbitrary level 0: 


(2) The average number of upward (or downward) axis-crossings per 
second. 

(22) The conditional probability that an upward axis-crossing occurs 
between t + 7 and t + 7 + dr given a downward axis-crossing at t. 

(272) The conditional probability that a downward axis-crossing occurs 
between t + r andt + 7 + dr given an upward axis-crossing at t. 

(wv) The conditional probability that an upward axis-crossing occurs 
between t + r and t + 7 + dr given an upward axis-crossing at t. 

(v) The variance of the number of axis-crossings observed in a time r. 


The theoretical probability functions are presented in graphs as a continuous 
function of 7 for various values of 6 and ‘‘a’’ for the case when the Gaussian 
noise has a Gaussian power spectral density. 


I. INTRODUCTION 


Consider the stationary random process /(é,a) consisting of a sinusoidal 
signal of amplitude ~/2a and frequency fy plus Gaussian noise ,y(é), 
of unit variance, having a narrow-band power spectral density W,(f — fo) 
which is symmetric about f, . Rice’s’ graphical representation of J(t,a) 
is illustrated in Fig. 1 in order to define the Rayleigh envelope process 
R(t,a) and the resultant phase —7 S O(t,a) S wr. The purpose of this 
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Fig. 1 — Graphical representation of I(t,a) = V2a cos 2xfot + Iy(t) = R(t,a) 
cos [2rfot + @(t,a)]. The point P wanders around, as time goes by, in the plane of 
the figure and generates the phase process 6(t,a). 


paper is to present some theoretical results concerning the axis-crossing 
points of the stationary, discontinuous phase process 6(t,a). In the 
literature, these same points are also called level-crossings. The axis- 
crossing points and the axis-crossing intervals of 6(t,a) are defined 
in I"ig. 2. The axis-crossing points and the axis-crossing intervals of 
R(i,a) are defined in a similar manner and were discussed by Rice’ 
and Rainal.*** The Rayleigh process R(t,a) and the phase process 
6(t,a) are present at the output of the II amplifier of a typical radio 
or radar receiver during the reception of a sinusoidal signal immersed 
in Gaussian noise. Also, the phase process @(t,a) is basic in Rice’s” 
recent analysis of noise in I'M receivers. 

Using a notation consistent with Refs. 3 and 6, we shall present 
the following theoretical results, in terms of well-known tabulated 
functions, concerning the axis-crossings of @(¢,a) at an arbitrary level 6 
and arbitrary signal-to-noise power ratio ‘‘a’’: 


(2) Ne, the average number of upward (or downward) axis-crossings 
per second. 

(22) QU(7,0,a) dr, the conditional probability that an upward axis- 
crossing occurs between ¢ + 7 and ¢ + + + dr given a downward axis- 
crossing at t. 

(t7t) Q4(7,6,a) dr, the conditional probability that a downward 
axis-crossing occurs between ¢t + 7 and ¢ + 7 + dr given an upward 
axis-crossing at t. 

(iv) [U,(7,86,a) — Q,(7,6,a)] dr, the conditional probability that 
an upward axis-crossing occurs between ¢ + 7 and ¢ + 7 + dr given 
an upward axis-crossing at ft. 


Il. AVERAGE NUMBER OF AXIS-CROSSINGS PER SECOND 


N,, the average number of upward axis-crossings per second of the 
level 6 by the phase process @(t,a), follows directly from some results 
due to Rice. Rice’ showed that 
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Ng / dR / do’ 6'P(R, 6,6’), (1) 
0 0 
where 
ee co R? _ (@'RY ¢] 
PiR,0,0") = == ox | 8 — “98 + QR cos 6 — 9 


= V2a 
an? [ "Wilf — fof — fo)? af 


I 


B 
—7 S087. 
W.(f — fo) = one-sided narrow-band power spectral density of Jy(é). 
Performing the integrations we find that 


va 


No = exp [— asin’ 6]®(Q cos 6), (2), 
where 
1 [ —y?/2 
og) = aes f) 
(x) A/D an y. 
Iquation (2) was also derived by Tikhonov’. 

Since 6 = 0 is a level of symmetry we have that Ny» = N_,. Also, 
the average number of downward axis-crossings per second is given 
by the right-hand side of (1) with the upper limit of integration of 6’ 
set to —o. Thus, the average number of downward axis-crossings 
per second is also equal to Ny. 
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Fig. 2 — The level 6 defines the axis-crossing points and the axis-crossing in- 
tervals of the discontinuous phase process @(t,a). 
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When the level 6 = +-7 and ‘‘a” is large, 2N 47 represents the average 
number of clicks observed in a time 7 at the output of an ideal FM 
receiver’’® during the reception of a unmodulated carrier in the presence 
of receiver noise. The variance of the number of clicks observed in a 
time 7 is discussed in Section IV. 

III. CONDITIONAL PROBABILITY FUNCTIONS 


The reader should refer to Rice” for the definition of all notation 
which is not defined in this paper. For the phase process @(t,a), the 
conditional probability Q5(7,6,a) dr, the conditional probability that an 
upward axis-crossing of the level 6 occurs between ¢ + 7 andt + 7+ dr 
given a downward axis-crossing of the level 6 at t, is given by an equation 
analogous to Rice’s” equation (86): 


GG) Dar aN, i. dR! i: dR! [ dR, - dR, 
ee ee ) 
ff. dé! a d0; of6tp(R, , Ri, 0, 0 , Rx , RE, 6, 0), 
where 
p(R, , RE, 0, 6, Rs, RE, 6, 0) 

_RER? aL : ; 
= 10 On)" exp {— aM Bar + R2 — 2Q(R, + R.) cos 6 + 2Q"] 
+ 2M,.,[R,R{ — RR, — Q(R{ — Rb cos 6+ Q(R, 6, — R,63)sin 6] 
+ 2M,,[R,R3 — RR; — Q(RS — Ri) cos 6 + Q(R.63 — R,61) sin 4] 
+2M,,[R,R. — Q(Rz + R,) cos 6 + Q*] 


+ M[Ri + Ro + Rie? + R363"] + 2M,,[R{Ri + R aatoy |}. 
The M’s are given in Rice’s’ Appendix I with . 
m() = [Wal = fr) cos 2x(f = fae af. (4) 
By performing the integrations with respect to R{ and R we find that 


Qi(r,0,a) = —Nz' | aR, | aR, 
: : (5) 


0 ee) 
.f ae, | d9; 610!p(R, , 0, 8 , Rs , 8, 0), 
—0 0 
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where 
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p(k, ) 4, 6 , fk, ) 6, 63) = 





+ 2Q sin 6M, — My )[R,6! — r08 |} exp (—G,/2M) 
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+ 20(Q — R, cos 6 — R, cos 6)(M1,, + M,;) 
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-((R, — Q cos 6)” + (Rz — Q cos 6)"] 


j= Mj.M os =a Mi,Mos “fr 2M 2M 3M 22) 
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-(2(R, — Q cos 6)(R, — Q cos oy}. 


By introducing the variables x, y, in place of @/ , 65 with the following 
transformation 
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p= of ate [1D ain 
Finally, by introducing the new variables 2, yo in place of R,, Rz 
with the following transformation 
R, = x + Q cos 8 (9) 
Re = Yo + Q cos 8, (10) 


we find, after some simplifications using Jacobi’s® theorem, that 


Qi (7, 8,a) 


n= 





(11) 
|x, ho), 


—2asin® 6 


= -1 fxd 24-3 : 
= [2rN,] [1 m]*MoeJ(r , hi) exp | 1+ m 


where 


1 ice] eo re 
K(m, ho) onV1 — ne i. dx i dYo € 


_ 2 + Yo — 2mrYo 
2(1 — m’) 


= m’ 1— ay : 
h, = -4 ; sie -|| Mes sin @. 


The conditional probability Q{(7,0,a) dr, the conditional probability 
that a downward axis-crossing of the level @ occurs between ¢ + 7 
and ¢ -++ 7 + dr given an upward axis-crossing at t, is obtained from (8) 
by changing the signs of the ’s in the limits of integration of 0/ and 63. 
We find that Q{(7,4,a) is equal to the right-hand side of (11) with h, 
replaced by —h,. This latter result also follows from the symmetry 
relation Q7{(7,0,a) = Q](7,—9,a). 

The conditional probability [U,(7,6,a) — Q,(7,9,a)] dr, the conditional 
probability that an upward axis-crossing occurs between ¢ + 7 and 
t + 7 + dr given an upward axis-crossing at t, is obtained from (38) 
by changing the lower limit of integration of 6{ to +-«. We find that 
U,(7,0,a) — Q,(7,0,a) is equal to the right-hand side of (11) with the 


20 = 
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function J(r, , h,) replaced by the function J,(r, , i), where 


1 =o foo) . 
he wal, dx a we=WG@=he G2 


The conditional probability that a downward axis-crossing occurs 
between ¢ + 7 and ¢ + 7 + dr given a downward axis-crossing at ¢ 
is obtained from (3) by changing the upper limit of integration of 
6, to —o. The result is that this conditional probability function is 
equal to the conditional probability function U,(7,6,a) — Q,(7,9,a) as 
one would expect from symmetry. 

The functions J(r,, hi), K(m,ho), and Ji(r:, hi) are expressed in 
terms of Karl Pearson’s *°’'**” tabulated function (d/N) in Refs. 2 and 3. 
Thus, the conditional probability functions Q7(7,6,a), Q{(7,8,a), and 
U,(7,0,a) — Q,(7,0,a) are expressed in terms of well-known tabulated 
functions. 

Since 6 = 0 isa level of symmetry, we need only discuss the conditional 
probabilities when @ is restricted to the interval 0 < @ < wz. The 
corresponding results when @ is in the remaining interval —7 S @ < 0 
can be deduced from the following symmetry conditions: 


Qi(7,,a) = Qi(7,— 6,a) (13) 
U,(7,6,a) — Q,(7,0,a) = U,(7,—86,a) — Q,(7,—9,a). (14) 


IV. VARIANCE OF THE NUMBER OF AXIS-CROSSINGS IN A TIME 7 


éc ) 


For an arbitrary level @ and arbitrary signal-to-noise ratio ‘‘a, 
let N(r,6,a) denote the number of axis-crossings observed in a time r. 
Then, we have that 


EN(r,0,a) = 2Nor (15) 
and 
Var N(r,0,a) = EN’(r,0,a) — [2Nor]’, (16) 
where 
iE = Expectation 


Var = Variance. 


Using McFadden’s” general result, also see Rice’s derivation in 
Bendat,'* we have that 


EN?(r,0,a) = 2Nyt + 4Ny i G2OUGIHis 09 
0 
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In this latter equation, U,(r,6,a) dr denotes the conditional probability 
that an axis-crossing occurs between t + 7 and ¢ + 7 + dr given an 
axis-crossing at time ¢. Since the joint probability that an axis-crossing 
occurs between t + 7 and ¢ + 7 + dr and an axis-crossing occurs between 
¢ and ¢ + dt can be expressed as 


2N,U,(7,0,a) dt dr = 2N,[U,(7,0,a) — Q,(7,0,a)] dt dr 
+ N,Q%(7,6,a) dt dr “|= N.Qi(7,0,a) dt dr, 


(18) 


we have that 


U(r, 9,a) a [U.(7,6,a) i Q,(7,6,a)] + $Q1(7,6,a) + 3Qi(7,6,a). (19) 
Thus, Var N(7,9,a) can be computed by using (16), (2), (17), and (19) 


Var N(r,0,a) = 2Nor + 4N, [ (r — x)U,(x,6,a) dx — [2Nyr]’ (20) 
0 


I 


2N eo +2 [ E as * live. ~ 2Nj] aw}. (21) 
For large 7, (21) becomes 
Var N(7,6,a) = Neo + 2 | [U,(2,0,a) — 2No] awh. (22) 


When twice the value of the integral in (22) is small compared with 
unity we have that 


Var N(r,0,a) = 2Nor. (23) 


This is the relation one would expect if the axis-crossing points represent 
a poisson point process for which U,(7,0,a) = 2N, for all +. 

Rice’ assumed a poisson point process for the case 6 = m and “a” 
large in order to use (23) in his analysis of noise in FM receivers. Indeed, 
for the case of a Gaussian autocorrelation function (22) serves to justify 
Rice’s use of (23) for large 7, 6 = z, and a = 4. For this case, with 
a = 4, numerical integration showed that the value of the integral 
in (22) is approximately 0.05. 

Notice that (22) not only applies to the point process defined by 
6(t,a) but also applies to more gencral stationary point processes. 

Incidentally, the probability function U,(7,6,a) can also be used to 
compute, approximately, the probability density p(7,0,a) of the axis- 
crossing intervals z; by using the following basic integral equation of 
renewal theory: 


po(t,0,a) = Ui(r,0,a) — po(r,6,a) * U,(7,6,a). (24) 
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The symbol * denotes the convolution operation, that is, 
jeg= | sar - dae. (25) 


Equation (24) is based on the assumption that a given axis-crossing 
interval is statistically independent of the sum of the previous (m + 1) 
axis-crossing intervals for all non-negative integral m. A theorem in 
Paragraph 5.2 shows that the assumption is false when m = 0. Thus, 
(24) can only yield approximate results. 

The exact probability density of the axis-crossing intervals x; is at 
present unknown. However, the first moment of this probability density 
is equal to [2N4]~* 


V. SOME SPECIAL CASES AND A THEOREM 


In this section we shall state some special cases of the conditional 
probability functions. We shall also present a theorem concerning the 
dependence of two successive axis-crossing intervals. 

5.1 Large r and Fixed 6, a 

As 7 becomes large we find that Qj (7,6,a), Q}(7,0,a), and U,(r,0,a) — 
Q,(7,9,a) approach the value NV, as one would expect. 

5.2 Small r and Fixed 0, a 

By expanding m(r) as 











3 b, : 5 ? 
mr) =1-8 4" fety bry beled (26) 
we find that as r — 0 from the right with b; ¥ 0 
Qilr,8,a) > Qi(r,8,0) > 28 tbe | BV Pe a ar | (27) 
127 
(Gia = Qi(r,8,0) > 2 [eve | (28) 
127 


Equation (28) suggests that wiggles having infinite rapidity and infini- 
tesimal amplitude are associated with the phase process 6(t,a) when 
b; ¥ Oor Wi(f — fo.) = O(f *) asf @. 

We also find that for small 7 with b; = 0: 


“Ge J(1,hi)7 exp | -£ (rQ sin a | (29) 





Qi(7,9,a) = 


746 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1967 


Q*(r,6,0) = oe Jie |-é (7Q sin ay | (30) 


U.(7,0,a) — Qi(r,0,a) = mee Ji(r; , hy)r exp | (7Q sin a |, (31) 


where 


ee BQ sin @ 
V4 ae B° 


It is interesting to compare the above results with the corresponding 
results at the level J of a Gaussian process J(t) having the normalized 
autocorrelation function m(r). That is, compare the above results 
with Rice’s’ equation (63) or Rainal’s® equations (44), (52), (53), and 
(54) when J = Q sin 6. The results are identical. 

Thus, a theorem’ concerning the dependence of two successive axis- 
crossing intervals of the Gaussian process J(¢) also applies to the phase 
process 6(t,a). That is, if @(t,a@) is a phase process, defined in paragraph 
one, having a finite expected number of axis-crossing points per unit 
time at any level 6, then two successive axis-crossing intervals at that 
level 6 are statistically dependent. 

The theorem implies that successive axis-crossing points do not 
form a Markov or Poisson point process. 


5.3 Q*(7,0,a) for small r, bs = 0, and large Q sin 0 
For small 7 and large Q sin 6 with b; = 0 or W,(f — fo) # O(f-*) 
as f — ©, we find from (80) that 
Qi(r,0,a) = FQ sin 6)” exp | -2 (7Q sin ay |. (32) 


Thus, Q{(7,0,a) is approximated by a Rayleigh probability density 
identical to Rice’s’ equation (65) when I = Q sin 8. 
5.4 a = 0 and arbitrary 6, 7 

When a = 0 we find that 

Qi(7,0,0) = Qi(7,9,0) 


a pee iV a5 
= 26°7(1 — m]° on? [rr 


+ V1 — rl — cos m] 





— cos ' 1) (33) 
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U,(r,0,0) — Qi(r,0,0) = 267 [1 — m2]? oa eee 
+ V1 ~r le — cos” ml, 
where 
0S cos 1X7 
0S cos' mr. 
Thus, when a = O the conditional probabilities are independent of 


the level @ as one would expect. 

5.5 Large a, 6 = 0, and arbitrary 7 
When ‘‘a’” is large and 6 = 0 we find that 
Qi(7,0,a) = Q1(7,0,a) 








(35) 
- art 2,-3 Mos -1 Pe) 
= 8 [1 — m |]? on [71(7 — cos nyt Vi —2] 
U,(7,0,a) oo Q,(7,0,a) (36) 
= B41 — my} ve [—r, cos r, + V1 = ri]. 


Thus, when ‘‘a” is large and 6 = 0, the conditional probabilities are 
independent of er 

Again, it is interestine to compare the above results with the corre- 
sponding results at the level J = 0 of a Gaussian process J(t) having 
the normalized autocorrelation function m(r). That is compare the 
above results with Rice’s’ equations (62) and (85a). The results are 
identical. One would expect identical results from Rice’s’ equation (3.6). 


5.6 6 = mw and arbitrary a, r 


When @ = a we find that 
Qi (7,7,a) Q3(7,7,a) 
= [6'(1 — m’)?o(—Q)|" 


I 


oe [r.(r — cos’ r,) 


+ Vil-nlK(m,Q) (7) 
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U,(7,7,a) — Q,(7,7,a) 


Fi u Moe 


= [6(1 — m)'a(—@Q)]| * [—r, cos rt V1 —r]K(m,Q). (38) 


VI. RESULTS FOR A GAUSSIAN AUTOCORRELATION FUNCTION 


For purposes of computation we shall take W,(f — fo) and m(r) as 


follows: 
Wi ~ ) = b= exp |G | (39) 


I 


and 


m(r) = exp | = eee)". (40) 


This particular selection was also made by Rice’ and Rainal* in their 
study of the duration of fades associated with the Rayleigh process 
Rita). 

From (40) we see that it is convenient to define normalized time as 
Uy = 2,7. All our results are plotted with respect to normalized time 
ts. The units of NV, are now ‘‘crossings per unit of normalized time.” 

Figs. 3 through 11 present the resulting conditional probability 
functions for various values of the level @ and for various values of 
signal-to-noise power ratio ‘‘a’’. For large values of wu, all of the 
conditional probability functions approach the value of N, in accordance 
with Paragraph 5.1. 

Figs. 9 and 11 compare Q*(7,0,a) for 6 = 7/2 and a = 4, 10 with a 
corresponding Rayleigh density in accordance with (382). Thus, we 
conclude that the Rayleigh probability density is a 2 Gone approximation 
when 7 is small and Q sin @ = +/2a sin 6 = 2 

Fig. 7 compares well with Figs. 2 and 3 of Ref. : Thus, we conclude 
that (35) and (36) are good approximations when a 2 4. 


VII. CONCLUSIONS 


The theoretical probability functions Q7{(7,6,a), Q%(7,8,a), and 
U,(7,0,a) — Q,(7,0,a) are expressible in terms of well-known tabulated 
functions. These results can be used to compute Var N(z,6,a), the 
variance of the number of axis-crossing points observed in a time r. 
These results can also be used to compute, approximately, the prob- 
ability density of axis-crossing intervals x; via renewal theory. The 
exact probability density is at present unknown. 
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Fig. 3 —Plots of the probability functions Qi7(w, 6, a), Qit(us, 9, a) and 
Ui(us, 8, a) — Qi(us, 6, a) associated with the axis-crossing points defined by the 
level 6 of 6(t,a). 6(t,a) denotes the resultant phase of a sinusoidal signal plus narrow- 
band Gaussian noise having autocorrelation function m(r). ‘“‘a’’ denotes the signal-to- 
noise power ratio. 
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Fig. 4 — Plots of the probability function Q,~(us, 6, a) associated with the axis- 
crossing points defined by the level @ of 6(t,a). 6(t,a) denotes the resultant phase of 
a sinusoidal signal plus narrow-band Gaussian noise having autocorrelation function 
m(r). “a” denotes the signal-to-noise power ratio. 
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Fig. 5 — Plots of the probability function Qit(w, 6, a) associated with the axis- 
crossing points defined by the level 6 of 6(t,a). 0(t,a) denotes the resultant phase of 
a sinusoidal signal plus narrow-band Gaussian noise having autocorrelation function 
m(r). “a” denotes the signal-to-noise power ratio. 
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Fig. 6 — Plots of the probability function U:(u,, 6,a) — Qi(us, 0, a) associated 
with the axis-crossing points defined by the level @ of 6(t,a). @(t,a) denotes the 
resultant phase of a sinusoidal signal plus narrow-band Gaussian noise having auto- 
correlation function m(7r). “a” denotes the signal-to-noise power ratio. 
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Fig. 7 —Plots of the probability functions Q:-(us, 6,a), Qit(us, 6, a) and 
Ui(us, 0,4) — Qi(us, 6, a) associated with the axis-crossing points defined by the 
level 6 of 6(t,a). 0(t,a) denotes the resultant phase of a sinusoidal signal plus narrow- 
band Gaussian noise having autocorrelation function m(r). ‘a’ denotes the signal- 
to-noise power ratio. 
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Fig. 8 —Plots of the probability function Q:-(us, 6, a) associated with the axis- 
crossing points defined by the level 6 of 6(#,a). 0(t,a) denotes the resultant phase of 
a sinusoidal signal plus narrow-band Gaussian noise having autocorrelation function 
m(r). “a” denotes the signal-to-noise power ratio. 
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Fig. 9 —Plots of the probability function Q:+(us, 6, a) associated with the axis- 
crossing points defined by the level @ of 6(t,a). 0(t,a) denotes the resultant phase of 
a sinusoidal signal plus narrow-band Gaussian noise having autocorrelation function 
m(r). “a”? denotes the signal-to-noise power ratio. 
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Fig. 10 — Plots of the probability function Ui(uz, 0, a) — Qi(us, 6, a) associated 
with the axis-crossing points defined by the level @ of 6(¢,a). 6(t,a) denotes the 
resultant phase of a sinusoidal signal plus narrow-band Gaussian noise having 
autocorrelation function m(r). “a’ denotes the signal-to-noise power ratio. 


752 


SINE WAVE PLUS NOISE 753 





t 
a=i0 





u 
m(7T) =EXPI|- 3S 





Qy (up, 6, a) 



































ie) 0.5 1.0 1.5 2.0 25 3.0 He) 4.0 4.5 5.0 5.5 6.0 
Up=27 Opt 


Fig. 11 —Plots of the probability function Q:+(us, 8, a) associated with the axis- 
crossing points defined by the level @ of 6(t,a). 6(t,a) denotes the resultant phase of 
a sinusoidal signal plus narrow-band Gaussian noise having autocorrelation function 
m(r). “a? denotes the signal-to-noise power ratio. 


Because the level @ = 0 is a level of symmetry, results forO < 6S 7 
imply results for —r S 6 < 0. 

When W,(f — fo) = O(f *) asf > ©, wiggles having infinite rapidity 
and infinitesimal amplitude are associated with the phase process 6(é,a). 

When @ = 0 with the signal-to-noise power ratio a = 4, the con- 
ditional probability functions associated with the phase process 6(t,a) 
are equal, approximately, to the corresponding results for a certain 
Gaussian process. 

When W,(f — fo) # O(f-*) asf > ~, and Q sin @ is large, Q}(7,4,a) 
for small 7 is approximated by a Rayleigh probability density. 

When @ = + anda = 4, Var N(r,6,a) for large 7 is equal, approxi- 
mately, to 2Nor, the variance resulting from a poisson point process. 

When WN, is finite, two successive axis-crossing intervals of @(é,a) 
are statistically dependent. Thus, the axis-crossing points do not 
represent exactly a Markov point process or a poisson point process. 
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Operating Characteristics for a Linear 
Detector of CW Signals in 


Narrow-Band Gaussian Noise 


By G. H. ROBERTSON 
(Manuscript received July 14, 1966) 


This article presents a set of curves that can be used to evaluate in detaal 
over a wide range of operating conditions the performance of systems 
using envelope detectors to search for CW signals in narrow-band Gaussian 
noise. The charts discussed relate the probability, Pp , of detecting such 
signals to (2) the signal-power /noise-power ratio, S/N, (ti) the proportion of 
false detections, Pp, , and (iit) the number, M, of independent samples 
of the envelope of the combined signal and noise that are averaged in making 
one attempt at detection. The curves and scales were calculated and drawn 
entirely by computer. The computation program was designed so that 
nearly linear curves could be produced, thereby increasing the accuracy 
and ease of interpolation. 


I, INTRODUCTION 


The curves* given in this article relate statistical properties of the 
output of an envelope detector to the signal-power/noise-power ratio, 
S/N, of a CW signal in narrow-band Gaussian noise at the input. 

The signal will be represented here by a finite section of a sine wave 
which may be divided further into segments, for convenience in proc- 
essing, so that most of the signal energy will be contained in a rela- 
tively narrow band, comparable in width to the reciprocal of the seg- 
ment length. A filter will be used to select the part of the spectrum in 
which the signal may be found. The output of this filter will be ap- 
plied to an envelope detector, and the amplitude of the output of the 
detector will be measured. Only if the magnitude measured at this 


* This work was supported by the U. 8. Navy, Bureau of Ships under contract 
no. N600(63133)64940. 
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point exceeds a chosen threshold value will a signal be presumed 
present. 

For other detection tests two or more samples of the output of the 
detector, taken far enough apart in time that the noise values are in- 
dependent, will be averaged. Only if the value of such an average 
exceeds a chosen threshold will a signal be presumed present. The 
curves given here for the averages of multiple samples depend on the 
assumption that a detected signal remains steady while the samples 
are taken, and that the noise background is weakly stationary. 

This form of detection criterion has been discussed by several au- 
thors,” **»# who give major emphasis to the use of a square-law de- 
tector because it is easier to analyze. A linear detector, however, is 
much more practical and is extensively used, so results that specifically 
apply to it are desirable. 

The curves available in the literature are difficult to use for checking 
system performance at more than a few scattered points. Those given 
here, on the other hand, are believed to be suitable for defining system 
performance in detail over a considerable range of threshold values 
and, furthermore, are easy to use. 


II. DISCUSSION OF THE MATHEMATICAL MODEL 


We assume that the spectrum of noise in the vicinity of the signal 
is smooth enough so that the output of the filter can be considered a 
narrowband Gaussian process with the possible addition of a sinusoidal 
signal. S. O. Rice®° derived an expression for the probability density 
of the envelope of such a waveform. Rice’s expression, paraphrased, is 


P(p) dp = 2h*p exp [—h'°(° + p’)]Io(2h’rp) dp (1) 
which gives the probability that the envelope lies within an interval, 
dp, of p. P(p) is thus the probability density. Io(z) is the modified 
Bessel function of zero order defined by 


= (Ba) 


1@) = DL Gas 2) 





Note that in (1) the term h’r’, which represents the signal/noise power 
ratio, and 1/~+/2h, which represents the rms noise level, can be identified. 
In this article S/N always means the ratio of the signal power to the 
noise power accompanying it in the specified narrow band. 
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When there is no signal (1) becomes 
P(p) dp = 2h*pexp (—h’p’) dp (3) 
which represents a Rayleigh distribution having 
20° = 1/h’ . (4) 


for the second moment about zero. 

The distribution of the envelope of the filter output is thus given by 
(1) when a signal is present, and by (8) when there is no signal. From 
these expressions distribution curves can be drawn for cases where 
only noise is present, and where a known signal occurs with the noise. 
To generate the performance charts given here, one has to be able to 
calculate the area lying above some threshold for distribution curves 
corresponding to no signal, and for those corresponding to all the values 
of S/N of interest. 

Only in the case of noise alone, for single-sample detection, is it 
possible to get an explicit formula for the area under portions of such 
a curve. Moments of the distribution given in (1) can be calculated 
with the help of Ref. 7. We give here the nth moment about zero: 


Vn 


[ Poe ap 


2 exp (—h’r’) i h’p"** exp (—A’ p’)Io(2h7rp) dp (5) 
0 





me 2 
(V2eyr(" a ).r,( md: -17'). 
In this h?r? is the signal/noise power and 


20° = 1/h? 


as stated earlier. 
T(z) 1s the Gamma function and ,;F, (a;b;z) is the Confluent Hy- 
pergeometric function defined by 


a(a + 1) 2 
b(b 1) 2-1 
aa+1\a+2) 2 


TS + 16 42) o21 °° 





iF (a;b3z) = 1 + Pi A 
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When the signal/noise power becomes zero we have 


= (Vaoyr(2 +2), (7) 


This gives the nth moment about zero for the distribution due to noise 
alone. We therefore have 


ie {2 P13 - i") (8) 


which gives the nth zero moment of the distribution when S/N = h?r’. 
We can thus derive a type-A Gram-Charlier® (G-C) series approxima- 
tion to the distribution curve for any S/N, and, since the series can be 
integrated, the area under any desired part of the curve can be found. 
Knowing the moments of the primary distribution, we can calculate 
the moments of a secondary distribution formed by averaging M sam- 
ples of the primary distribution. We thus form a G-C series which can 
be integrated to get the area under any desired part of the secondary 
distribution curve. 

Charts are included that cover cases in which, when only noise is 
present, the probability is as low as 10 that the threshold will be ex- 
ceeded, even with the use of single samples in detection attempts. 

Producing a G-C series that would give results within 1 percent of 
the true value when the threshold is great enough to give such a low 
probability required evaluating many more moments than is usual in 
G-C series. In fact, 31 moments were calculated, which resulted in a 
series that oscillates about the true value of the distribution as the 
threshold is varied. A method that sometimes reduces these oscillations 
is to use special groupings of the terms, forming a so-called Edgeworth 
series.* It occurred to the author, however, that a better fit might be 
obtained by averaging the approximations given by several G-C series 
carried out to different numbers of terms since their oscillations as 
the threshold varies are not in phase. This involves very little more 
work than deriving the series that has the most terms, since it merely 
requires that the terms be weighted to represent the number of times 
they occur in forming the average. Fig. 1 shows the results obtained by 
forming G-C series that use up to 31 moments and averaging over the 
highest 23 and 138 terms, respectively. Both show less oscillation over 
restricted—though appreciable—ranges of threshold than does the 
Edgeworth series that uses up to 30 moments. The curve showing 
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CURVE ABOVE IS LOG Q 
(READ RIGHT SCALE) 
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(READ LEFT SCALE) 
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THRESHOLD (RAYLEIGH DISTRIBUTION) 


AVERAGING SERIES FOR 9 THROUGH 
31 MOMENTS 


B AVERAGING SERIES FOR 19 THROUGH 
31 MOMENTS 


EDGEWORTH SERIES TO MOMENT 30 





THRESHOLD —> 


Fig. 1 — Comparison of methods for evaluating G-C series. 


the logarithm of Q, the function approximated by the series, allows 
the useful range to be assessed. It makes possible the calibration of the 
horizontal scale, which is the same for all four curves, in terms of 
the probability that the threshold will be exceeded. The logarithm to 
the base ten of this probability is given by the right scale of Fig. 1. 


III. DESCRIPTION OF THE CHARTS 


We will now describe briefly how the charts (Figs. 2 through 15) 
were produced and point out their main features. Some of these are 
believed novel and show the advantage of being able to use a digital 
computer to generate and draw charts of this type. 

For each value of S/N at which a curve was desired, the detection 
threshold, 7’, was expressed in two ways: (2) measured with respect 
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to the mean of the distribution for noise only, and standardized by 
dividing by the square root of the variance, and (22) measured with 
respect to the mean of the distribution appropriate to the S/N, and 
standardized with respect to this distribution. 

The form in (t) is appropriate for calculating the probability of 
false alarm (i.e., the probability that noise alone will cause the thresh- 
old to be exceeded), and the form in (7) is appropriate for calculating 
the probability of detection (the probability that a signal at the chosen 
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Fig. 2— ROC curves. 
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Fig. 3— ROC curves. 


S/N will cause the threshold to be exceeded). The relationship between 
these two forms is linear, so that if the threshold were varied a straight 
line would give the relationship between a scale for measuring the 
variation based on (1) and another based on (iz). If these scales were 
calibrated in terms of probability of false alarm, Pr4 , and probability 
of detection, Pp , useful charts could be obtained that use these prob- 
ability parameters as horizontal and vertical axes, respectively. 
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But, for any Pr, scale, the corresponding Py scale depends on S/N. 
Fortunately, the Pp scale does not change rapidly with S/N for the 
range we are interested in, and, to allow several S/N values to be in- 
cluded on each chart, a constant normal probability scale was chosen 
for Py. The penalty incurred by doing this is the slight curvature of 
some of the S/N lines plotted in Figs. 2 through 15. 

It can be seen that the horizontal scale, P;4 , changes as the num- 
ber of samples averaged changes. This is because the probability dis- 
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Fig. 4— ROC curves. 
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Fig. 5— ROC curves. 


tribution for the average of a number of samples depends on the num- 
ber averaged as well as on the basic distribution. Since an explicit 
formula for the cumulative probability is known in the case of the 
basic (Rayleigh) distribution, the Rayleigh distribution was used to 
get the Pra scale for the chart based on single-sample decisions. All 
the other scales were calculated from G-C series approximations to the 
appropriate distributions. 
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To make interpolation easier, a column of heavy dots was put near 


the right edge of each chart. These dots mark 0.2 dB increments in 
S/N between the lines. Since any two adjacent S/N lines are very 
nearly straight and parallel it is easy to interpolate to within 0.1 dB 
between the lines using a parallel ruler guided by the heavy dots. 


The horizontal and vertical grids were plotted as rows and columns 


of dots in order to carry the accuracy of the scales into all parts of 
the chart without also producing a confusing mesh of lines. 
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Fig. 6— ROC curves. 
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Fig. 7— ROC curves. 


Each decade of the horizontal, Pr4 , scale, reading from left to right, 
is given by 


10.0 (1.0) 4.0 (0.5) 2.0 (0.2) 1.0, 


where the increments are in parentheses. 
The vertical, Pp , scale, reading from top to bottom, is 
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0.05 (0.01) 0.10 (0.02) 0.90 (0.01) 0.950 (0.005) 0.980 
0.980 (0.0025) 0.990 (0.001) 0.998 (0.0005) 0.999, 


the increments again being in parentheses. 


th 


On Fig. 16 are curves for two pairs of Pp , Pra values that show how 
e S/N for each pair varies with the logarithm to the base 2 of the 
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Fig. 8— ROC curves. 
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Fig. 9— ROC curves. 


number of samples averaged. Over any restricted range in which the 
number of samples is doubled, either curve could be approximated 
quite satisfactorily by a straight line. 

The slopes of these two curves are almost the same at the same S/N, 
although the Pp, Pr4 values are quite different. This means that, 
within any small range of Pp, Pr4 values, curves like those shown in 
Fig. 16 can be assumed parallel with little sacrifice in accuracy. Con- 
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sequently, the change in S/N corresponding to a small change in log,.M 
indicated by such a curve will hold quite accurately for a small range 
of Pp, Pra values around the pair for which the curve is actually 
drawn. 


The S/N corresponding to any pair of Pp , Pra values can therefore 


be calculated for a different number of samples averaged by assuming 
a linear relationship between S/N and the logarithm to the base 2 of 
the number of samples. 


256 INDEPENDENT SAMPLES AVERAGED 
S/N FROM -8.0 TO -2.0 IN 1dB INCREMENTS 


0.050°::°:: 28 8 were + ae 


0.100-:: tinea dates Gaara nebo bine 


BAO SE ween 


O2300 88 Sas Stk eas aety 


Ocd008 Sees Se OU oes thane tated Mule Moe ha ci beaten ae. 
Os BOOS ees = Soe Melee ZZ, atatvken eee eG Sar ceabsine hecvecncs aise hts teaehec tates 


0.600-°-: ----y 


0.700---: , 


O.900s He cr ttt frre cee teen een a sens teens reens sean se farrnee tres sete sceee sees mecsores sees fle mers coos 


0.950°+-: +: 


0.975: - 


0O.990°-**+: "7°: 


0.995°-:> 





FA 


STANDARDIZED THRESHOLD (UNIT IS dy) 


Fig. 10 — ROC curves. 
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Fig. 11— ROC curves. 


IV. CONCLUSION 


4.1 Lfficiency of Post-Detector Integration 

When more than one sample can be averaged to improve detection 
sensitivity the process is usually called post-detector integration. The 
curves of Fig. 16 show how efficient this process is at different values 
of S/N and Pp, Pras operating points. The gradients of these curves 
show the rate at which detectability is improved, measured in dB/ 
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double-the-number-of-samples, ds, for two operating points. This ratio 
is convenient to use in interpolating between adjacent charts, since the 
number of samples doubles from one chart to the next. 

It can be seen that the gradients are almost the same where the 
S/N is the same. The rate of improvement is better than 2 dB/ds for 
S/N above 10 dB, presumably reaching a maximum of 3 dB/ds for 
high S/N. The rate apparently becomes constant at 1.5 dB/ds for S/N 
values less than —10 dB. 
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Fig. 12— ROC curves. 
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Fig. 13 — ROC curves. 


4.2 Extending the Scope of the Charts 


The detector input is restricted to a relatively narrow band of fre- 
quencies to ensure that the envelope can be measured accurately after 
detection. By bandshifting high enough, this result can be obtained 
even for bands of quite appreciable width. The charts can therefore 
be used in situations where the significant signals may be quite short 
bursts of nearly sinusoidal waveform. 
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The same standardized threshold scale applies to all the charts. It 
is given in units of 


On ss o/ VAM, (9) 


where o is the rms noise into the detector and M independent samples 
are averaged. 

Although the curves given here refer strictly to envelope detectors 
they can be applied when square-law-detectors are used if a small 
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Fig. 14— ROC curves. 
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Fig. 15 —— ROC curves. 


error (always less than 0.2 dB) in S/N is taken into account, and the 
standardized threshold scale is changed. A calibration curve for this 
error 1s given in Ref. 1. Reasonable agreement with this curve was 
found by generating ROC curves for a square-law detector after suit- 
ably modifying the computer program. It was possible thus to get 
ROC curves for 1 sample and for averaging 128 or more independent 
samples. 
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Fig. 16 — Sensitivity related to number of samples averaged. 
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Phase and Amplitude Modulation in High- 
Efficiency Varactor Frequency Multipliers — 


General Scattering Properties 


By C. DRAGONE 
(Manuscript received April 20, 1966) 


The presence of phase and amplitude modulation in the signals of a 
varactor frequency multiplier is analyzed, and some general multiplier 
properties are derived. The following summarizes one of the most important 
results of this paper. 

Consider a frequency multiplier which has the following characteristics: 
(2) a varactor which has a square-law characteristic, (it) the order of multz- 
plication is N = 2" = 2, 4, etc., (tit) the minimum number of idlers, and 
(iv) zt 1s lossless and tuned. It is shown that for this multiplier there ts no 
conversion between small-index, low-frequency amplitude and phase modula- 
tion. Therefore, since narrow-band noise from external sources will be 
present at the input and output of the multiplier, the noise components 
corresponding to phase modulation of the carriers can be treated inde- 
pendently of the noise components corresponding to amplitude modulation. 

Consider now the input and output noise sidebands corresponding to 
phase modulation (PM). It is shown that the multiplier behaves with 
respect to these sidebands as an amplifier with the following properties: (7) 
a forward voltage transmission equal to N, (it) unity reverse transmission, 
(272) an output reflection of 1, —1, 3, respectively, for N = 2, 4, 8, and 
(iv) no input reflections. As a consequence of these properties the multiplier 
ts “potentially” unstable with respect to PM. 

The utility of the multiplier properties derived in this paper will be 
illustrated by the discussion in a companion paper in this tssue which 
shows how, in practical cases, instability arises and how it can be avoided. 


I. SUMMARY OF RESULTS 


The frequency multipliers to be considered are harmonic generators 
which use varactor diodes as nonlinear elements. Noise produces un- 
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wanted amplitude and phase modulation in the signals of a frequency 
multiplier. In this paper, the presence of these modulations is analyzed 
and some general multiplier properties are derived. The following 
summarizes one of the most important results of this paper. Consider 
a frequency multiplier which has the following characteristics: 


(<) a varactor which has a square-law characteristic, 
(iz) the order of multiplication is N = 2", 
(77) the minimum number of idlers, and 
(iv) it is lossless and tuned. 


It is shown that for this multiplier there is no conversion between 
small-index, low-frequency amplitude, and phase modulation. Therefore, 
since narrow-band noise from external sources will be present at the 
input and output of the multiplier, the noise components corresponding 
to phase modulation of the carriers can be treated independently of the 
noise components corresponding to amplitude modulation. 

Consider now the input and output noise sidebands corresponding 
to phase modulation (PM). It is shown that the multiplier behaves 
with respect to these sidebands as an amplifier with the following 
properties: 


(zt) a forward voltage transmission equal to JN, 

(ii) unity reverse transmission, 
(zit) an output reflection of 1, —1, 3, respectively, for N = 2, 4, 8, and 
(iv) it has no input reflections. 


As a consequence of these properties the multiplier is ‘‘potentially”’ 
unstable with respect to PM. This summarizes one of the most important 
results of this paper. Now let us consider the results of this paper in 
more detail. 

This paper is concerned with the presence of amplitude and phase 
modulation in the multiplier signals. Suppose, for the moment, that 
these two types of modulation are independent of each other. That is, 
suppose that the multiplier does not produce AM-to-PM conversion, 
and vice-versa. Suppose, furthermore, that only PM is present. Then 
each signal will consist of a carrier and of a pair of sidebands in quad- 
rature with respect to the carrier. Since either sideband can be obtained 
from a knowledge of the other, then one may consider only one of the 
two sidebands and ignore the other one. For example, let the upper 
sideband be chosen as the variable, and let it be described in terms of 
propagating waves. Then the input variables of the multiplier are the 
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two waves* vy; )U5,4, Which constitute the upper sideband of the input 
carrier wy . Similarly, the two output waves vy.,45) Uvo.+p represent the 
output variables. wy is the input ‘‘carrier’’ frequency of the multiplier 
and p is the frequency of the fluctuations. At this point, the scattering 
formalism furnishes a convenient way of describing the properties of the 
multiplier. More precisely, one may define the PM scattering param- 
eters of the multiplier as the reflection and transmission coefficients 
which relate the “scattered” waves v7,,,, Vvo,+> to the ‘‘incident”’ 
ONES Vp 4>) Yvo,+p- Ln this way one obtains the PM scattering matrix 
S, defined by 


ve v, UN 
SotP =? 2 WotD ye DN wotp 
+7; 


ian ae View Vu 





(1) 
UNwortp = T, Voorn + po UN worn 


—= 146 0 ’ 
0 
Vue Vins Vue. 


where V.,, is the Fourier coefficient of the input carrier w) , and Vy., is 
the Fourier coefficient of the output carrier Nw). Note that the side- 
bands have been normalized with respect to the relative carriers. Equa- 
tions (1) are equivalent to (8). 

In a completely similar way one defines the AM scattering matrix S, . 

It has been found that a multiplier with characteristics (2) through 
(zv) above does not produce AM = PM conversion and that it has 
the following scattering matrices: 


a pp Ts 0 (1) 


So = —— (2a) 
ei = 1 — (—1)"N 
rox) |w Lo Guw 
of Tt Noel eT 
S, = _ 3N N (2b) 
/ ae 1 0 


if the modulation frequency p is small enough. 

Notice that (2a) gives the properties stated at the beginning of this 
summary. As a consequence of (2a), the multiplier is potentially un- 
stable with respect to PM. AM instabilities cannot occur. This follows 
from (2b). 


*y,, v4 designate the Fourier coefficients of the voltage components of fre- 
quency w. They propagate in the directions indicated by the arrows. 

{+ Note that v{,4,/V.o,, etc. correspond to the modulation indexes j0@,’, etc. 
defined in the next section. 
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Equations (2a), (2b) also have some other interesting consequences. 
For instance, if one injects a single tone Nw) + p into the output port 
of a multiplier, then the multiplier will reflect a pure PM wave. That 
is, the multiplier behaves as an ideal reflection-type limiter in the vi- 
cinity of the output carrier. This property is common to all types of 
multipliers [see (19)] and may be useful if one wants to modulate the 
output phase of a multiplier without generating AM. Note that 
|p; | > 1 for N > 4. This means that, if N > 4, then the output port 
of the multiplier reflects amplifying the PM components of the output 
sidebands. 

It is also shown that some of the properties expressed by (2a), (2b) 
are common to all lossless multipliers. A consequence of this is that 
most of the results obtained for a multiplier with the characteristics 
listed above can be qualitatively extended to all efficient multipliers. 


Il. INTRODUCTORY REMARKS 


This paper is concerned with the presence of small amplitude and 
phase fluctuations in the electrical variables of the networks to be 
considered. By using the familiar terminology and concepts of modula- 
tion theory, the problem can be defined as follows. 

The electrical signals will consist of sinusoidal carriers which are 
both amplitude- and phase-modulated. The amplitudes of the modulat- 
ing waves* will be supposed to be small enough to guarantee superposi- 
tion to hold, and to allow the modulated waves to be approximated 
by the sum of the carriers and their first-order sidebands. Furthermore, 
since superposition holds, consideration will be limited to sinusoidal 
modulating waves. Then, if p is the frequency of the modulating waves 
(i.e., of the fluctuations), each carrier frequency w will be surrounded 
by two small side-frequencies w + p, w — p representing the sidebands 
of w. One can say that the object of this analysis is to study the presence 
of these side-frequencies around the carriers. A pair of sidebands is 
completely specified either by its spectrum or by the spectrum of the 
modulating waves associated with it. More precisely, consider a small 
voltage v(t) consisting of two side-frequenciesf w + p, w — p 

* The “modulating waves” are the “fluctuations” and will be clarified later. 


{ Through this analysis a real function g(t) consisting of sinusoids will be repre- 
sented by the complex Fourier series 


g(t) = ¥ Clas) exp [je] = 2(Re) ¥ C(wi) exp [atl 


The complex coefficients C will be called “the Fourier coefficients’’ of g(t). Further- 
more, only positive frequencies will be considered, as illustrated by the second 
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v(t) = 2(Re)Wa+» exp [jw + p)t] + v.-» exp [j@ — p)E]) (3) 


and let V,(¢) be the carrier associated with v(t). Let w, V. be the carrier 
frequency and the Fourier coefficient of V,(é), respectively. Then the 
Fourier coefficients of the modulating waves a(t), @(t) associated with 
v(t) are given by 





° (4) 
= nS (Pass = vt) 
27\ VY. Vi 
In fact, if one uses (4) to express v(t) in terms of a, 6, one obtains* 
u(t) = 2(Re)[(a + 70)V. exp [7 + p)é] 
+ (a — j0)*V. exp [j@ — p)é]] (5) 


I 


a VQ + 2 vy, 


Since it is supposed | v,4, |< | V. |, one has|a|]<1,|6|<1. There- 
fore, the last equality of (5) gives 


1) + V0 = 0+ av e+ 2], (6) 


Equation (6) shows that a(é), 0(¢) are the modulating waves which 
produce the sidebands v(é).t a(t) is the amplitude modulating wave, 
and @(t) is the phase modulating wave. 

a, 6, the Fourier coefficients of a(t), 9(t), will be called, respectively, 
the AM index and the PM index of v(t). Extensive use will be made of 
the AM, PM indexes to represent a pair of sidebands.{ 


summation. This convention will not be followed in the Appendix where (and only 
there) both positive and negative frequencies will have to be considered for con- 
venience. Note that C(w:) = C*(—w,) because g(t) is real. The asterisk ( )* will 
always indicate the complex conjugate of a complex quantity. (Re) means “real 
part of’. 

* a(t) = 2(Re) (a exp jpt), 0(t) = 2(Re) (@ exp jpt), V(t) = 2(Re) (V. exp jut) 
according to the preceding footnote. 

t When a(t), @(¢) are applied to the carrier V(t). 

t As already noted, a, 6, represent the spectra of the modulating waves relative 
to v(t). It should also be noted that the representation of v(t) by means of a, @ cor- 
responds to the familiar representation of a pair of sidebands by means of the so-called 
‘fn-phase” and “quadrature” components about a given carrier. This is shown by 
(5) in which aV,., a*V, are the “in-phase”? components and j@V., j@*V. are the 
“quadrature” components. Therefore (a) can also be regarded as the normalized 
(with respect to 7V-) amplitude of the “quadrature” upper-sideband. 
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2.1 The Variables of the Interconnection Between Two Stages 

With reference to Fig. 1 consider now the interconnection between 
two general stages. Let V(t) be the interstage voltage and assume that 
V(d) consists of a sinusoidal carrier V,(f) and a pair of sidebands v(Z). 
Now suppose that a short piece of line of characteristic impedance 
Z, is inserted between the two stages. Its electrical length is assumed 
to be so short that it may be inserted without altering the electrical 
properties of the cascade connection of the two stages. This should be 
regarded as an artifice made to allow the signals of the interconnection 
to be described in terms of propagating waves.* The sidebands v(t) 


Volt) +v (t) 
v(t) 


V(t) =V,(t) +v(t) =V_(t) v(t) v(t) 


Fig. 1— Description of the interstage voltage in terms of propagating waves. 


now consist of two waves uv’ (é), v(t) propagating in opposite directions, 
as indicated by the arrows —, <, and 


ot) =v (4) +v Md). (7) 


Finally, consider the respresentation of v”(¢), v(t) by means of their 
PM, AM indexes with respect to the carrier V,. By using (4) one 
obtains} two pairs of modulation indexes 


a’, 6’, the indexes of v(/) 
a’, 6, the indexes of v*(Z) 


(a”, etc.) provide a complete description of the fluctuations present at 
the generical interconnection, as illustrated in Fig. 2. 


2.2 Definition of the PM, AM Scattering Matrices of a Stage 
In most cases to be analyzed{ each stage will have two basic prop- 
erties: 


* Z, will always be chosen such that only one propagating wave exists at the 
carrier frequency. Therefore, the voltage of the interconnection consists of a modu- 
lated wave propagating in one direction, and of a pair of sidebands without carriers 
propagating in the other direction, as it is illustrated in Fig. 1. The carrier V(t) is 
used as a reference for deriving the modulation indexes of both sets of sidebands. 


t By replacing v(t), v44p,04—p With v(t), Vo +p » Ve—p ONC Obtains a”, 6. Ina similar 


way one obtains a’, 0°. 
t In Ref. 1, AM = PM conversion is considered. 
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(z) it will be linear with respect to the modulating signals (since they 
are supposed to be small) 
(zz) it will not generate any AM-to-PM (or vice versa) conversion. 


As a result of property (iz) the PM indexes of a stage are independent 
from the AM indexes, and vice versa. Therefore, the two cases AM, 
PM can be treated separately. 

Consider for instance the PM case. With reference to the stage of 
Fig. 3, let 6; , 6% be the input (¢ = 1) and output (¢ = 2) PM indexes. 
Then, because of (7), the scattered indexes 67, 62 can be related to 
the incident ones 67 , 6% through a set of linear equations of the type 


6, = poi + T002 
6, = T560; + po82. 


(8) 


> 


a,6 


Fig. 2— Description of the interstage sidebands in terms of modulation in- 
dexes. 


Notice that (8) are equivalent to (1). The scattering matrix S, 
po Ts 
Ts po 


will be called the PM scattering matriz* of the stage. 7"; will be called 
the PM forward transmission coefficient. T’, is the PM reverse trans- 
mission coefficient. p, is the PM input reflection coefficient. p is the 
PM output reflection coefficient. 

In a completely similar way the AM scattering matrix S, of a stage 


is defined tf 


S, == (9) 








8, = (10) 





paw “Ly: | 
Dee: py. 

* Note that the complex coefficients p~, 77, T~, p~ correspond to the familiar co- 
efficients Si, Sie, So, S2 which are usually employed to represent scattering 
parameters. (See Ref. 2). As can be directly verified from (8), if a unit index is 
incident on either port, 7’ gives the index transmitted out of the other port, and p 
gives the index reflected at the same port. The arrows —, < are applied to 7’, p to 
indicate the directions of the incident waves to which 1’, p refer. 


+ The convention of using the letters a, @ for designating the AM, PM eases will 
always be followed through the analysis. 
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Fig. 3— Description of the PM properties of a stage by means of the PM 
scattering coefficients. 


In conclusion in this analysis: the variables are the PM, AM indexes; 
the frequency of the variables is p; the parameters specifying the behavior 
of the various networks are the PM, AM scattering coefficients. 


III. GENERAL PROPERTIES OF A LOSSLESS MULTIPLIER WHICH DOES NOT 
PRODUCE AM = PM CONVERSION 


It is important to emphasize that this analysis applies both to a 
multiplier consisting of a single stage and to a chain of multipliers. 

Furthermore, in this paper, the analysis will be limited to the case 
of a lossless multiplier which does not produce AM @ PM conversion. 
This is the most important case for the following two reasons. The first 
reason is that, in general, a multiplier which is tuned does not produce 
AM = PM conversion if the sidebands are close enough to the carriers 
(see Appendix). Note that it is assumed that the bias network is properly 
designed so that effects such as those described in Ref. 3 are absent. 
Furthermore, results obtained for the case of no AM = PM conversion 
can be extended to the general case of AM — PM, as it is shown 
in the second paper. The second reason is that the question of stability 
generally arises only when the losses of the multiplier are small. Further- 
more, if the losses are small, the results which are obtained by neglect- 
ing them can be readily extended to include them, as is pointed out in 
the second paper. 

In this section it is demonstrated that a lossless multiplier of order 
N has the following properties: 


Tol .* = N(1 — pop.* 
<7 il e- 
TTo* = 77 (Ll — popa* (11) 
Pa* po Pa* , Po 
fe 4 yl 9, NP 4% 9, 
TT; TT; 


The equations of (11) have three important properties. 
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One is that they depend on only one multiplier parameter: the order 
of multiplication NV. They do not depend in any way on the actual 
circuit configuration or on the type of varactor, etc. 

Another property is that they are frequency independent. That is, 
they do not contain p explicitly even though 7'> etc. will in general 
be frequency dependent. 

The last important property of (11) is that they represent 4 in- 
dependent relations among the 8 unknowns 7, etc. For instance, 
if the AM coefficients are known, then the PM coefficients can be 
calculated by means of (11), and vice versa. 


3.1 Demonstration of (11) 


Consider the multiplier of order N shown in Fig. 4(a) connected to a 
load Ry and driven by a voltage generator having impedance Rf, . 
Let V,(¢) (r = 1, N) be the voltage of the input port (r = 1) and out- 
put port (r = N) of the multiplier in the absence of fluctuations. It 
will be supposed that the input of the multiplier is matched to the 
generator. Therefore, the input impedance of the multiplier is R,. 
The multiplier is assumed to be lossless and without AM <@ PM. 
The characteristic impedance associated with the rth (r = 1, N) port 
is R, , according to the footnote on p. 780. 

Assume now that two small waves v,°(é), vy" (é) are arriving towards 
the multiplier, and that they consist of the frequencies w») + p, Nwo + 2, 
respectively. The problem to be considered is to find the scattered 





v; (t) vn (t) 


GENERATOR 
=, 
vitt) 


Fig. 4— (a) Multiplier of order N. (b) Input and output sidebands of the 
multiplier. 
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waves v,"(t), vy (2), as illustrated in Fig. 4 (in which the carriers V,(d), 
Vy(t) are not shown). 

Let (@;, 65), (a; , a=) be the indexes of v,(#) with respect to the 
carrier V,(é). Then, according to (5), 


v, (t) = 2(Re){(ay + 967) V, exp [jrwo + p)t] 
+ (a,* + 56/")V, exp [ito — p)t]} (12) 


and similarly for v,~(é) (just change the arrows directions). V, is the 
Fourier coefficient of V,(é) (r = 1, N). 

Now let P;.,,., be the power carried by the forward wave of frequency 
rwy + p, and let Pr,,_,, Pr.,-p» Pia.+p have similar meanings. Then 
the total power at the frequency w (w = rw -+ p) is the sum of the two 


components (—), (<-): 
P,=PZ+P°. (13) 


Consider now the fact that the multiplier is lossless, and that there- 
fore Manley—Rowe relations hold.* That is, 


Pye Pan) (Pease Prea-n_| 
otp o-D =, OEP oT P = 0). 14 
( Pee Wo — Pp Noo + p Noo — p ( ) 


From (12) one can calculate PZ , PZ (w = rw) + p). For instance, 


2 
2 LV, | 


OR. = |a, +96, | Po (15) 


Prete = AG gee. | 
in which | | indicates the absolute value. P, is the power delivered by 
the generator to the multiplier at wy) and transferred by the multiplier 
to the load Ry . (The input carrier power is equal to the output carrier 
power because the multiplier is lossless.) 

After calculating PZ , P , one obtains P, (w = rw) -k p) by means 
of (13). Then, by substituting P,, in (14), one obtainsf 


(IM)[N(@tar* — O7ay*) + Oxax* — Gxax*] = 0, (16) 


where (IM) indicates the imaginary part. 
Now if one considers the multiplier PM, AM scattering matrices 


* The minus sign in front of the second addend of (14) occurs because power is 
assumed positive when it is flowing towards the right. This, at the output port, is 
opposite to the usual convention of assuming positive the power absorbed by the 
varactor. See Ref. 4. 

‘ { In deriving (16) the approximation rw) + p & rey has been made so that (14) 
ecomes: 


NPootn — NPoo-p — Procite + Pnwo-» = 9. 
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S, , Sy defined in the previous section 


























Oy a! Oe ay 4 | ay 
= = BS ] 
Ox So Oy ’ ay Sa ay ( 7) 
and substitutes (17) in (16), one obtains 
(IM)[@rar*(V — Npop.* — T6T.*) 
+ Ovay*(1 = pe, NTT .*) (18) 


+ Gjay*(—Nocls* — pr*T%) 
+ Oyay*(—NTGp.* — poT.*)] = 0. 
A relation which has to hold for any choice of (67, 67, etc.) and 
therefore it has to hold for any choice of (67a>*,6yay*,etc.). This is 


possible only if (11) holds. Note that the approximation rwo + p & ror 
has been made in deriving (18). 


IV. SLOW-VARYING FLUCTUATIONS 


Suppose now that the input generator V, is phase modulated and that 
its amplitude is kept constant. If the modulation frequency p is small 
enough, then the output phase deviation will be just N times that of 
the input generator. Furthermore, all amplitudes will remain constant. 
These properties are common to all multipliers and are well known.°”° 
Therefore, they will not be demonstrated here. As a consequence of 
these properties one has* 

lim T; = N lim py = 0. 


p-0 p0 


By combining these two results with (11) one obtains 





TT =N pp = 0 
/ ae | pe = 0 (19) 
Tt ge = — 

NT; NT; 


which apply to the limiting case of slowly varying fluctuations. From 
(19) one can see that 7°; = 1. This has the meaning that if one amplitude 
modulates the input generator, then the output voltage will have an 
identical (percentage-wise) modulation. 

* The multiplier is matched at the carrier frequency wo. Therefore, there is no 
reflection present at the input of the multiplier if modulation is absent. Consequently, 


input reflections will remain absent if the input phase is slowly modulated. Therefore, 
Pg =Oifp20. 
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Note that (19) leaves only two coefficients to be determined: two 
of (p,*,p, ,7;) or two of (p;*,o; ,7.*). Both (11) and (19) represent 
quite general properties of ideal (i.e., lossless, tuned) multipliers. They 
do not depend on the actual circuit configuration of the multiplier, 
on the type of varactor characteristic, on the power level, etc.* Note 
that from (19) one has 


(To*Ta*)(LoT's) = 1. 


This means that either the “PM round-trip transmission” (7'37") 
or the ‘‘AM round-trip transmission” (737%) is larger than unity.f 
Therefore, reflections of the external circuit (the output load and the 
input generator) can produce instabilities.’ 


V. MULTIPLIER OF ORDER N = 2” 
In this section the following type of multiplier will be considered: 


(z) the nonlinear capacitance has a square-law Q-V characteristic, 
(iz) the order of multiplication is N = 2", and 
(iit) it has the least number of idlers’ (Qu), -++ , 2” *w). 


This type of multiplier is most important because it can be exactly 
treated with little difficulty and it is realistic at the same time. In fact, 
abrupt-junction varactors exactly satisfy (2); graded-junctions ap- 
proximately satisfy (2); most practical designs are based on (777) for 
reasons of simplicity. Condition (77) excludes from this treatment two 
important cases: VN = 3, N = 5. On the other hand, the (exact) results 
which will be obtained for N = 2,N = 4,N = 8, N = 16, N = 82, 
etc. may be qualitatively extended to the remaining cases N = 3, 
N = 5, ete. 

In the Appendix, it is shown that, because of (2), (2), and (277), one 
has 7%; = (—1)", pj = [N — (—1)"]/8N. Therefore, by combining 
these results with (19), one has 





8, —_ Pe Ts = 0 (—1) (2a) 
oe ee es bs 
Tf N wate 
a eae aay 
Pa T. OAT!” 
ee | ill. (2b) 
To pa 1 0 


* They are valid both for a single stage and a multiplier consisting of many stages. 
+ Neglecting the possibility that they are both unity. See the next section. 


GENERAL SCATTERING PROPERTIES 787 


Therefore, the multiplier scattering matrices are uniquely determined 
as functions of only the order of multiplication N. 

The physical meaning of the scattering coefficients (2) may be il- 
lustrated by the following two examples. 


5.1 First example. 


Let a small signal Nw, + p be injected into Ry by means of a direc- 
tional coupler connected as shown in Fig. 5(a). Then the voltage across 
Ry will be both amplitude and phase modulated, and the AM, PM 
indexes will have equal amplitudes according to (4). Let a be their 
amplitude. Next, reverse the connection of the directional coupler, 
as shown in Fig. 5(b). Now the following facts will be observed: 


(t) The voltage of Ry contains only PM. That is, the AM component 
of w + p is absorbed without reflection. by the output of the multiplier 
(because p’, = 0). The PM component, on the contrary, is reflected 
back into the load Ry . The PM index across Ry is pa with p> given 
by (2). Note that | >| = 1 for all values of N, and that | p;| > 1 
for N > 4. Therefore, the PM component is reflected with amplifica- 
tion for N > 4. For instance, if N = 8 the gain is 3. All this indicates 
the circuit of Fig. 5(b) as a possible scheme for modulating the output 
phase of a multiplier without producing any AM (and with gain if 
N > 4). 

(ii) The voltage across R, has both AM and PM. The PM index 
is Tja = (—1)"a and the AM index is 7ja= (—1)"/N-a. 






DIRECTIONAL 
77 COUPLER 
7 





Fig. 5 — Injection of a tone at the output port of the multiplier. 
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5.2 Second example. 


Now let the previous experiment be repeated at the input side of 
the multiplier, as shown by Fig. 6. Let a be the amplitude of the indexes 
generated across R, when the generator w + p is connected as shown 
in Fig. 6(a). Then, after connecting w + p as shown in Fig. 6(b), one has: 


(z) The AM index across Ry is Ta = a and the PM index is Tja = 
Na. These are general properties which have already been found and 
discussed in the more general case of the preceeding section [see (19)]. 

(iz) Across R, there is only AM and the AM index is pZa with p; 
given by (2). 


In practical applications it is important to bear in mind the following 
meaning of (2). If only phase modulation is present, then one can cal- 
culate the upper sidebands by simply replacing the multiplier with a 
linear and time-invariant amplifier whose scattering properties are 
given by (2a). Similarly, if only AM is present then one can calculate 
the upper sidebands by using (2b). All this is true provided the upper 
sidebands are normalized with respect to V, , in the AM case, and with 
respect to 7V, , in the PM case. 

Consider now the PM case. From (2a) one obtains that the round- 
trip PM transmission 737'% is 


| T3T; | = N. (20) 





Fig. 6 — Injection of a tone at the input port of the multiplier. 


GENERAL SCATTERING PROPERTIES 789 


Therefore, interactions between the multiplier and the other stages 
can cause phase instabilities.*Furthermore, from (20) one can see that 
the possibility of instabilities increases with the order of multiplication 
N. However, the presence of loss will decrease the round-trip trans- 
mission. More precisely, if y is the multiplier efficiency, f then the round- 
trip power transmission will become approximatelyt 


[Pers P= Ney? (21) 


and the possibility that interactions of the external circuit cause PM 
instability exists when approximately 


1 
1 >" (22) 
Note that amplitude instabilities cannot occur.§ In fact, if the output 
port of the multiplier is connected to an arbitrary passive stage, then 
the input AM reflection of the multiplier is always less than] 


lee | + [Pars |, (23) 


which is never greater than unity, as one may verify from (2b). 


APPENDIX 


Demonstrations for the multiplier of order 2” which has the minimum 
number of idlers and uses a ‘‘square-law”’ varactor. Slowly varying modula- 
tions. 


Let S,(t) be the total elastance of the varactor. It can be separated into 
the sum of a time-varying component S(t) and an average component S, 


S(t) = S@ + So. 


The ‘“‘external circuit’? connected to the varactor can be represented 
by an impedance Z, in series with a voltage generator V, , as it is shown 
in Fig. 7. 

Let Z(w) be the total impedance in series with S(¢). Then Z = Z, + 
So/Jw. 


* This is discussed in detail in the second paper. 

ft Output power divided by input power. 

tf »? is the round-trip transmission through an attenuator which has a forward 
(power) transmission 7. Equation (21) has been obtained by representing the 
multiplier as the cascade connection of such an attenuator and a lossless multiplier. 

§ It is important to note that it is assumed that the “bias circuit” is properly 
designed, so that low-frequency fluctuations of the average varactor capacitance Cy 
are avoided. 

{ This can be obtained by using standard techniques. 
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The hypotheses are: 


Z=0 if wX 2% (s=1,°-- ,n—-J) (24a) 
Z=R, if oLe (24b) 
Z=Ry if w= Noy = Bw (24c) 


Z=o if w isfarfrom 2’w (s=0,---,n). (24d) 


Equation (24a) requires Z = 0 in the neighborhood of the idler fre- 
quencies 2w), -++ , 2" ‘wo. Equations (24b), (24c) require Z to be 
equal to the input impedance R, for w & w., and to be equal to the 
output impedance Ry for w & Nw, . In (24d) it is required that current 
flow be limited to the frequencies of (24a), (24b), and (24c). Because 






iy oe 
ae A 
(v+v)=A(Q+q)? 
(a) (b) (c) 





Fig. 7 — Equivalent circuit of the frequency multiplier. 


of this last requirement, the variable component Q(t) of the charge 
of the varactor (when the sidebands are absent) is of the type 


Qo = E>, : Q, exp (jraol). (25) 


Since the varactor has a square-law Q-V characteristic F(Q), the 
total voltage* V,(é) across the varactor is related to Q(t) through a 
relation of the type 


FQ) = Vi) = AQ’ + QMS. + Vo. (26) 
Therefore, the time-varying part of V; will consist of two components: 
one component V(é) due to AQ(t)” 


_ * The sign convention in this appendix is that power is positive when flowing 
into the varactor. 
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Vit) = A D1 QQ, exp Lr + dwot] = D7 Vi exp (jluot) 
ifr 27) 
(@ = 42°; r= 42°; s,h =0,--- ,n) 


and another component V’ = Q(é)S. which represents the contribu- 
tion of a constant capacitance Cy = 1/S, . The total elastance is obtained 
by taking the derivative of (26) with respect to Q(t): 


S,(t) = S(t) + So = 24Q(t) + a (28) 


and one finds a constant capacitance Cy in series with a time-varying 
elastance component S(t) originating from AQ(t). Therefore, the 
equivalent circuit of the multiplier becomes that of Fig. 7(b) in which 
the nonlinear element has a Q-V characteristic of the type AQ’, and 
Vi), Q(), S(é) represent its voltage, charge, and elastance, respectively, 
when the sidebands are absent. From (28), (25) one has 


S() = 2AQ@ = >> S, exp (jrarol) (29) 
s=0,°%*,n 
r=+28 
with 
S, = 2AQ,. (30) 
Let now 
Qi) = De des exp [ite + ap) al (31) 
r=+28 = 


t=+1 


be the sideband components of the charge of the varactor. Then the 
voltage sidebands across S(t) are 


of) = Sa) = Di rx,+ exp [iro + ép)é]. (32) 
By substituting (29), (81), in (82) one obtains* »,,; 
Uy, 41 = x Seances (r = i; oa eye (33) 
u=+28 
But from (24d) and from the fact that Z(w + p) & Z(w), one has 
Q(r-u),41 = 0 for (r ae u) Fao +], eae +2”, (84) 


In (38), therefore, one has to consider only those values of wu which 
satisfy (put: r = 2”) 
2* —u = +1,--- , +2”. (35) 


* Note that q7,; = @*+,-1, Ur,; = v*,r,-7, (see footnote on p. 780). Therefore, con- 
sideration can be limited to the case r > 0. 


792 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1967 


Then, the only possible values of wu are 
u= or Uu= > Naas ur —2. (36) 


because all other values of wu (remember that wu has to be of the type 
u = +2°) would cause |r — u| to be either zero or odd. Therefore, 
remembering that r = 2", one finds that (33) consists of only three 
terms: 


V;,41 >= Gijarai ris -+ Gay eis + qe <1 Se, (r =1,-:- Ae (37)* 


In conclusion, in order that q(t), v(t) represent sideband components 
produced by sources located only at the input and the output of the 
multiplier, it has to be (necessary and sufficient condition) 


v.41 =O for r= +2,--- , &2”” (38a) 
Gr.41 #0 only for r = +1,--- , +2” (88b) 
Ur 41 = Qr/2, 4194/2 + Gi 4s + Qe e1So- for r = 16.5 ea (38c) 


Equation (38a) follows from (24a); (88b) follows from (24d); (38¢) is 
(87). 

Note that (88a) and (88b), are the constraints given by the “external 
circuit’? (which includes Cy), and (&8c) is the constraint given by the 
elastance S(t) of the varactor. 

The “equilibrium equations” for the ‘‘carriers’” Q, , V, are obtained 
in a similar way. They are 


V,=0 for r= +2,°---,+2"" (39a) 
Q, #0 onlyfor r= +1,--- , +2” (3Yb) 
Vi = jooQ.R, (39¢) 
Vy = —jNoQRn (39d) 
Ve = AlQ,2Q-2 + 2Q2,Q07%] for r=1,--+,2". (89e) 


Equations (39a) and (89b) follow from (24a) and (24d), respectively. 
(39c) follows from (24b) and the hypothesis that the input generator 
is matched to the multiplier. Equation (39d) follows from (24c). Equa- 
tion (89e) can be derived from (27) by using the same procedure used 
to derive (38c) from (32). 


* Note, if r = 1, the first term of the second member is zero; if r = N, only the 
first term of the second member is nonzero. 
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A.1 Demonstration of T; = (—1)"/N 

Theorem 1: If q ,d satisfies (88) and if 
Ge = Gra exp |i (-"| 

(40) 
Dy 41 — —0, 41 exp E (-1'| (r = 2°38 = 0, at ,n) 

then also q, 0 satisfles (38), 7.e., it gives sidebands produced by sources 

at wo & P; Nar =k p. 


Proof: @,% clearly satisfies (88a) and (88b) because g, é does. Therefore, 
it remains to be demonstrated that (38¢c) is satisfied. 
By substituting (40) in (88c) one obtains 


—8,,41 EXP |i ("| = Ger.eSe2 XD liz | 


+ or, 21S-» exp liz ("| ae GF =1Se, exp | -i% |, 


which is satisfied because 6, g satisfies (38), and because 


—exp Hey = exp isco | — exp | -i5-n° |. 


Suppose now that @, 6 is produced by a PM source located at the 
input of the multiplier. Therefore, remembering that p; = O because 
of (19), the input components of ¢, 6 represent a PM forward wave, 
and also the output components constitute a PM forward wave. Con- 
sider now the input (r = 1) components of @, 6. From (40) one has 


n1,41 >= — 1,21 Kiss = JQ1. 41 (41) 


from which one can see that, since 6,41, g1,4: is an FM “forward” 
wave, then 0,,.:, 91,2: is an AM ‘backward’ wave. In a completely 
similar way one finds that also the output components dy ,21, Gv,21 
represent an AM “backward wave’. Finally, remembering (4) and 
the definitions of 7; , 7, , one has 


= On bya Vi 
T; = (2 - —= : A 
@n=0 ’ 











0; Vy Oia 

— ay 011 Vu 
/ (a). Sa ee 
Gn/ a1=0 1 UN, 
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By combining this result with (40) one has 


LU es ee TT, = (—1)". 


1.1 oy ,1 
But 73 = N. Therefore, one obtains the desired result 
-_ 0, 
P= V 
A.2 Demonstration of p, = [N — (—1)"]/N-3 
~ _N — (-1)" 
po = SCY (42) 
First it will be shown that the charge g(é) defined by 
Ort = Oejnt = (1 — p(s))Q, (r = aa s=0,°-:- , n) (43) 
with 
oer es —] n=8s 
p(s) = re 3 (44) 


gives the charge sidebands produced by an AM wave of amplitude 
a; = 1 arriving at the input of the multiplier. 

Let first the fact that g(é) satisfies (38) be demonstrated. Substitute 
(43) in (37) and take into account (30). One obtains v, the voltage 
associated with q, 


Year = 2A{Q.2Q,/2(1 — p(s — 1) 
+ Q2,Q--[(L — ep) + 1 — os + 1)}} 
(r= 23s =1,--- ,n—])). (45) 
From (44) one has 
2 — p(s) — pS +1) = 201 — p(s — I). 
Therefore, by using (39e), (45) gives 
v.21 = 4(1 — p(s —1))V, =0 for r=2,--- ,2"” 


because of (39a). Therefore, g, v satisfies (88). Furthermore, it represents 
AM because the sidebands g,,., are “‘in-phase’”’ with the carriers Q, . 
This is shown by (48). 


Now let the output indexes ay, ay of q, v be calculated. Consider 
the following expressions: (38c) with r = N; (30) with r = N/2; (48) 


GENERAL SCATTERING PROPERTIES 795 


with r = N/2; (89e) with r = N. From them one obtains 
Vie oa Vig . 
From this first relation, by using (39d) and (43) with r = N, one obtains 
Dy.1 > —Ry(jNwogn 1). 
This second relation gives ay = 0, because jNwoqy,. is the current 
at Nw + p. Then the first relation gives at once ay = 1. 
From the fact that ay = 0, ay = 1 and from the fact that (19) gives 
T, = 1, the input indexes must be 
a, = l, a; = pa. 


Therefore, the input voltage sidebands must be 


Ds a1 ag ¢ =F paVi . (46) 
But, by substituting (43) in (37) with r = 1, one has 


D1 = (1 = p(1))Q2S-1 ae (1 or: (0)) QiS2 


from which, remembering that Q.S._, = Q*S,; = V, because of (39e) 
and (30), one has 


Yi = (2 — p(1) — 0(0))Vi = 1 + 0) Vi (47) 


in which use has been made of (44). Finally, by combining (47) with 
(46), one obtains 


NS te) 


pa = p(0) = 3N 


which is the desired result. 


A.3 Demonstration of the fact that the multiplier does not produce AM ss 
PM conversion, if p 1s small enough. 


An input PM source does not cause AM if p is small enough. In 
fact, one may verify that the sidebands produced by an input PM 
forward wave of amplitude 67 = @ are 


Qr,1 = rQ,8, Qr,-1 = jrQ,0* (48) 


which are in ‘‘quadrature’’ with respect to the carriers Q,. To dem- 
onstrate that there is no PM — AM conversion it therefore remains 
to be shown that an output PM source does not produce AM. This 
can be done by applying transformation (40) to the AM sidebands 
given by (43). In fact, in this way one obtains PM sidebands. They are 
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produced by PM sources located both at the input and the output 
of the multiplier. Therefore, one concludes that there is no PM — AM, 
ifp=o0. 


By using the same procedure one can show that there is no AM — 


PM conversion. In fact, the discussion following (48) shows that an 
input AM source does not produce PM. Then, by applying (40) to (48), 
one finds that an output AM source does not produce PM. This con- 
cludes the demonstration. 
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Phase and Amplitude Modulation in High- 
Efficiency Varactor Frequency Multipliers of 
Order N = 2’-Stability and Noise 


By C. DRAGONE 
(Manuscript received April 20, 1966) 


A general analysis of the stability conditions of frequency multipliers 
of order N = 2” ts presented. The frequency multipliers to be considered 
are harmonic generators which use varactor diodes as nonlinear elements. 
The type of instability investigated 1s that which causes spurious tones 
to appear at the output of a multiplier in the vicinity of the desired harmonic. 
It arises because an efficient multiplier 1s potentially unstable with respect 
to the quadrature components of its sidebands. 

This paper shows how instability arises and how it can be avoided. 
One of the main results is that, to obtain stability in practical cases, it 
ts sufficient that the bandwidths of the various resonant circuits satisfy 
some simple conditions. 


I. SUMMARY OF RESULTS 


A general analysis of the stability conditions of frequency multipliers 
of order N = 2” is presented. The frequency multipliers to be considered 
are harmonic generators which use varactor diodes as nonlinear elements. 

That stability is one of the most serious problems in high-efficiency 
multipliers is a widely known experimental fact.'’” 

At the present time little is known of the restrictions placed by the 
condition of stability on the available circuit configurations. Con- 
sequently, present design procedures leave the problem of stability 
to be solved experimentally, and this is often done at the expense of 
efficiency. Furthermore, multipliers which are individually stable may 
become unstable when connected together to form a chain. As a result, 
isolators are often needed. The isolators will guarantee stability but 
will lower the overall efficiency. 

The practical importance of the problem to be analyzed is illustrated 
by the following example. 
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Suppose that one has designed a stable and very efficient octupler, 
and that one wants to reduce the output noise by using a bandpass 
filter consisting of a high Q resonant circuit connected in series with the 
load. In general, one will not obtain the desired result. In fact, the out- 
put noise will, in general, increase rather than decrease. Furthermore, 
if the bandwidth of the filter is too narrow, then the multiplier will 
become unstable.* 

This example illustrates the important fact that, as the circuit 
approaches an unstable condition, the output noise level increases 
indefinitely. Then, when the multiplier is on the point of becoming un- 
stable, the output noise becomes very large at some frequencies. Further- 
more, when instability arises, spurious tones appear in the vicinity 
of the carriers. Therefore, not only is it important that a multiplier be 
stable, but it is also important that it be far from instability, if one 
wants a low-noise multiplier. 

This paper shows how instability arises, and how to avoid it. It also 
shows how to derive the output noise from a knowledge of the various 
noise sources. Some of the results which have been obtained are sum- 
marized by the following statements. 


1.1 Doubler 
Consider the doubler first. 


(z) The stability of the doubler will not, in general, depend on the 
impedance presented by the input circuit (to the varactor) in the 
“vicinity” of the input carrier w) . More precisely, in the design of the 
input circuit, consideration can be limited to those frequencies whose 
distance from the carrier w. is larger than half the bandwidth of the 
output circuit [roughly; see (44), (45), (46)]. At those input frequencies 
which are not in the ‘vicinity’ of the input carrier, 

(77) the impedance presented by the input circuit to the varactor 
should be large rather than small, as compared to RF, [see (13), Theorem 2, 
Theorem 3], where R, is the impedance presented by the input circuit 
at wo . Furthermore: 

(zat) The output bandwidth should be large compared to the input 
bandwidth [see (35), ete.]. Furthermore: 

(iv) An efficient chain of more than two doublers which are in- 
dividually stable will, in general, be stable if each doubler has been 
designed according to (77) and if, in addition, each doubler is sufficiently 
broadband with respect to the preceding one [see (58)]. Furthermore: 


* This example is derived at the end of this paper. See (51), (52), ete. 
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(v) Those circuit configurations should be preferred which produce 
PM = AM conversion [see (55)]. Therefore, low-pass circuit configura- 
tions at the input, and high-pass circuit configurations at the output, 
are, in general, preferable to bandpass circuit configurations. 


1.2 Multiplier of Order N = 2" 


Consider a multiplier of order N = 2” > 2 which has the minimum 
number of idlers. This memorandum shows that such a multiplier is 
equivalent to a chain of doublers.* Therefore, results obtained for the 
doubler can be extended to this type of multiplier. For instance (2) 
gives: 


An efficient multiplier of order 2” which has the minimum number 
of idlers should have 


B,KBoK+++ KBr KB, 


where B, is the bandwidth of the input circuit, B. is the ‘‘equivalent”’ 
bandwidth of the first idler, etc. 

Note that the above results} apply to the case of a lossless multiplier. 
The presence of losses reduces the restrictions placed by the condition 
of stability. If 7 < 1/N, then the question of stability does not arise 
any more, in general. This has been shown by Ref. 3. 


II. PRELIMINARY CONSIDERATIONS 


2.1 Method of approach of the mathematical description of the multiplier 
and its signals 


This paper is concerned with the presence of amplitude and phase 
fluctuations in the multiplier signals. Suppose, for the moment, that 
these two types of fluctuations are independent from each other. That 
is, suppose that the multiplier does not produce AM-to-PM conversion, 
and vice-versa. Suppose, furthermore, that only PM is present. Then 
each signal will consist of a carrier and of a pair of sidebands in quad- 
rature with respect to the carrier. Since either sideband can be obtained 
from a knowledge of the other, then one may consider only one of the 
two sidebands and ignore the other one. Let then, for instance, the 

* See Theorem 1 of the second section. Note that the equivalence is exact only 
if the varactor has a square-law Q-V characteristic. 

+ These results implied that the impedance presented by the circuit (to the 
variable capacitance) is “large’’ at the frequencies which are far from the various 


carriers. Therefore, subharmonic oscillations and bias instabilities are not included 
in this analysis. 
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upper sideband be chosen as the variable, and let it be described in 
terms of propagating waves. 

Then the input variables of the multiplier will be the two waves* 
U5.+p) Yae+p Which constitute the upper sideband of the input carrier 
wo. Similarly, the output waves vy.,,1,) Yvo.+p represent the output 
variables. wy is the input ‘‘carrier’’ frequency of the multiplier, and p 
is the frequency of the fluctuations. At this point the scattering formalism 
furnishes a convenient way of describing the properties of the multiplier. 
More precisely, one may define the PM scattering parameters of the 
multiplier as the reflection and transmission coefficients which relate 
the “scattered” waves v,,.,, Uve.+p to the “incident” ones v;,,, , 
VNwe+p- One obtains in this way the PM scattering matrix S, of the 
multiplier. 

Through the analysis v2,4,, User» Unoocep) Ynootp» Have been nor- 
malized with respect to the carriers. In this way the variables become 
the dimensionless coefficients 


~ - ~ — 
Veotp 50, eo Vootp UNwotr iOS nan UNwo+p 
: 1 = N ~~ 


Vo, y] ae b jOn = Vinee ? Vuwe d 
which represent ‘‘modulation indexes’. Note that V., is the Fourier 
coefficient of the input carrier w) , and Vy,, is the Fourier coefficient 
of the output carrier Na . 

In a completely similar way the AM case is treated. By normalizing 
the AM upper sidebands with respect to V.,, Vy, , one obtains four 
modulation indexes aj, a], @y , @y which represent the AM variables. 
Furthermore, by considering how these waves are scattered by the 
multiplier, one finds the AM scattering matrix S, which describes the 
multiplier AM properties. 


in = 


2.2 General Scattering Properties of an “Ideal” Multiplier 


Consider a multiplier of order N = 2” which has the following prop- 
erties: it is lossless, it has the minimum number of idlers (n — 1), and 
it uses a varactor having a square-law Q-V characteristic. The first 
paper (See Ref. 3) has shown that, if the fluctuations are slow enough 
(p = 0), then such a multiplier does not produce AM $3 PM conversiont+ 
and it has the following scattering matrices: 


*v,, vs designate the Fourier coefficients of the voltage components of fre- 
quency w. They propagate in the directions indicated by the arrows. For more 
details on the mathematical description of the multiplier and its signals see Ref. 3. 

{ It is important to note that it has been assumed that the ‘‘bias circuit’’ of the 
multiplier is properly designed, so that low-frequency fluctuations of the average 
capacitance Cy are avoided. See Ref. 4. 
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Va ee. 








pa Vy. aa 

8, = S| ae (1) 
PP op. 1 0 

“ a ee 0 —j)" 

S, = Po 6 =. ( ) (2) 
Ts pe Nn i = 1)"N 


A multiplier for which (1), (2) apply will be called ‘‘ideal’’. According 
to the preceding discussion, one has that, if only PM is present and if 
a multiplier is ‘‘ideal’’, then 





Veotp = pe Vaotp + Tl UNwotp 
—_ 6 0 
Veg Ve Vue. (3) 
Ve aes: fA Vendy -_ pe 
ee = 15-55 4 Po a= 
Vues Von Vuwe 


If only AM is present, the relations existing among the upper sidebands 
are similar to (3) (Replace pj, 7%, etc. with p; , T’, , ete.). 

In conclusion, one can calculate the sidebands present in an ‘‘ideal’’ 
multiplier in the following way. First, separate the PM components 
from the AM components. Next, solve separately the two cases AM, 
PM. The PM case is solved by means of (3). That is, one replaces the 
multiplier with a linear amplifier which has 


(c) a forward transmission coefficient 7"; , 
(iz) a reverse transmission coefficient 7'5 , 
(272) an input reflection coefficient p, , and 
(zv) an output reflection coefficient p%, . 


Then one supposes that the variables consist of the PM upper side- 
bands only, as it is illustrated by (3), and one readily calculates them. 
In a completely similar way one calculates the AM upper sidebands. 


2.3 General Considerations on Stability 


The frequency response of an unstable circuit which is on the point 
of becoming stable is infinite.* Hence, if the parametersf of a circuit 
are continuously varied, a circuit can go from a stable to an unstable 
situation if and only if its frequency response becomes infinite. There- 
fore, a certain situation A is stable iff 

* See Ref. 5, p. 316, or Ref. 6, p. 112. 


+ These parameters are discussed in the next sections. 
t Note that there are stable circuits for which (4) is not satisfied (conditionally 
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one can cause the circuit to pass from situation A to a stable 
situation B by continuously varying its parameters, without (4) 
causing its frequency response to become infinite. 


For instance, consider a multiplier which is individually stable and 
which is connected to a generator and a load which have reflections. 
Suppose, furthermore, that one is concerned about the possibility that 
these reflections cause instability. Then, according to (4) one may find 
out whether or not the circuit is unstable (or conditionally stable) by 
examining whether or not the circuit response becomes infinite as the 
reflections are decreased in amplitude. This can be done by inserting 
attenuators between the multiplier and the other stages (the load and 
the generator). In this example, situation B occurs when all reflections 
become zero. The importance of (4) and of this example will be better 
understood in the next sections. 

Consider now the case of a multiplier which does not produce AM @ 
PM conversion. Then only phase instabilities can occur, as has been 
shown by Ref. 3. Therefore, one is concerned about its behavior with 
respect to PM only. 

In general, when the multiplier is on the point of becoming unstable, 
all four PM scattering coefficients (7~,T",p’,p~) become infinite si- 
multaneously.* Therefore, one may choose any one of them to study 
the stability of the multiplier. However, since the forward transmission 
T~ is of special interest because it gives the output response to input 
PM signals, it will be selected as the multiplier frequency response to 
be analyzed. 

The stability of a multiplier depends on the values of certain param- 
eters which will be discussed in the following sections. Those values of 
these parameters for which the circuit is stable constitute the so called 
“region of stability”. At the boundary of the region of stability the PM 
forward transmission becomes infinite at one or more frequencies, as 
has already been pointed out. Therefore, in order to find the stability 
region, one must find those “critical”? values of the parameters for 
which 77” = o. Then, once the boundary is found, one must identify 
which one of the two regions separated by the boundary is stable. This 
can be done by applying the stability test (4) or any of the usual stability 


stable circuits. See Ref, 6, p. 162). However, for practical reasons, one is interested 
in designing a multiplier which does not become unstable if the losses are increased. 
(i. ey which is unconditionally stable.) For such a multiplier (4) is necessary and 
sufficient. 

* See Ref. 6, p. 164. Since the discussion will be mostly confined to the PM case, 
the subscripts a, 6 used in (1), (2), and (3) to distinguish the AM case from the 
PM case will be omitted unless necessary. 
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tests (such as the ‘Nyquist plot’, etc.) to one point of either one of 
the two regions. 

As will be seen, a multiplier can be represented by means of a chain 
of stages which are individually stable. The analysis of such a chain 
may be carried on as follows. Consider the cascade connection of the 
two stages (b), (c), shown in Fig. 1. If the two stages are individually 
stable the overall forward transmission becomes infinite if and only if 


PoPe = 1 (5) 


where p*, is the output reflection of the first stage* and p; is the input 
reflection of the second stage. Equation (5) may be obtained in the 





Fig. 1— “Loop transmission” at the interconnection between two stages. 


following way. Consider the interconnection between the two stages 
(b), (c), and suppose that the reverse path is separated from the forward 
path, as illustrated in Fig. 1(a). Next, break the reverse-path connec- 
tion as it is shown in Fig. 1(b). In addition, terminate P, in the char- 
acteristic impedance Z, of the interconnection, and apply a unit voltage 
to P,, as illustrated in Fig. 1(b). Then the voltage which appears 
at P. is p;p.. Therefore, pp; represents the ‘‘return voltage” (i.e., 
the familiar loop transmission p@)} of the loop indicated in Fig. 1(a), 
and (5) is demonstrated. 

In some cases, one will be interested in knowing the reflection p~ 
presented by the output of a stage when its input port is connected 
to a generator which has an output reflection p> + 0 (see Fig. 2). 
The output reflection p~’ is given by | 


/ 


* Which occurs when the other port of the stage is terminated in its characteristic 
impedance. 
1 See Ref. 6, p. 44. 
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Fig. 2 — Effect of p> on the output reflection of a stage. 


pS py. ee 5 (6) 
Lap Os 
where p , p', 7’, T~ designate the scattering coefficients of the stage 
when its ports are terminated in the respective characteristic impedances. 
Equation (6) may be obtained by using standard techniques.’ 


III. EFFECT OF INTERACTIONS BETWEEN A MULTIPLIER AND TWO 
PASSIVE STAGES 


Consider a multiplier connected between two stages as shown in 
Fig. 3(a). Such a problem is encountered in two practical cases. One 
is the design of the input and output networks of a multiplier, as will 
be seen in the next section. The other case occurs when the load and the 
generator have reflections and one is concerned about the possibility 
that these reflections may cause instability. This is the case considered 
here and it is illustrated in Fig. 3(b). It should be clear, however, that 
results obtained for either case can be extended to the other one. 

In this section it is supposed that the multiplier is lossless. Further- 
more, it is assumed that interaction only occurs in the vicinity of the 
carriers* so that one can approximate the multiplier properties by means 
of (2) (which is valid if losses are absent and if p = 0). It is assumed, 
in addition, that the generator and the load do not produce AM = PM 
conversion. This is the worst case, as is explained in the next section. 

Since there is no AM = PM conversion, one is only concerned about 
phase instabilities and, therefore, it will be assumed that there is only 

* This will be justified by the results obtained in the next section. In fact, the 
next section shows that, if the multiplier is ‘“‘properly”’ designed, then its PM trans- 
missions (7'~, 7") are maximum for p & O and decrease monotonically with p. 


Therefore, interaction is in general more likely to occur for p & 0 (i.e., for those 
values of p for which (2) holds), in practical applications. 
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PM. Therefore, consideration can be limited to the upper sidebands only 
as it is indicated in Fig. 3 where Z,(w.) + p), Z~,(Nwo + p) are the im- 
pedances of the generator and the load at the upper-sidebands, respec- 
tively. Let p, , p,* be the reflection coefficients of the generator and the 
load, respectively. Then 


= Z (wo + p) — R, = Z(Nowo + p) — Ry 
‘ Z,(wo + p) + R, we Z1(Nowy + p) + Ry’ 


where R, , Ry are the values of Z, , Z, for p = 0.¢ Now let p’ indicate* 
the reflection presented by the multiplier to the load when the input 


(7) 





Z_(N@o +p) 





~1e 


PL[™ 


- | 
xN 
~1Pg 
(b) 


Fig. 3— Multiplier connected between two stages. 
port is connected to the generator. From (6), taking into account the 
fact that now p’ = 0 [see (2)], one has 
p =p +p," T, (8) 


where p", p , 7’, T” designate the multiplier PM scattering coefficients. 
From (2) one has 


p= (-1)" 


* The arrows will be omitted, unless necessary. 
+ R, is the input impedance of the multiplier at w) (the multiplier is assumed to 
be matched: Z,(wo) = 1). 


Sa TT" = (—-1)N. (9) 
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According to (5), infinite transmission occurs if (and only if) 


orp’ = 1. (10) 


Equation (10) is obtained from (5) by considering the cascade con- 
nection of the generator and the multiplier as stage (b), and by con- 
sidering the load as stage (c). By substituting (8) in (10), and by taking 
into account (9), one obtains 


sieht (—1)'po( Wp, = N= (Cw) og: (11) 


According to the stability test given by (4), one has that the circuit 
is unstable depending on whether or not at some frequencies* 


(Re)(pip’ — 1) 2 0 
(IM){ez0’] = 0. 


This can be shown by applying (4) as illustrated by the example of 
the preceding section. In fact, consider a circuit for which (12) occurs 
at some frequencies. Next, insert an attenuator} between the load and 
the multiplier. Then, as the attenuator is varied from zero to infinite 
attenuation, the forward transmission becomes infinite, because of (12), 
(11). Therefore, (12) guarantees instability. Next, consider a circuit 
for which (12) never occurs. In this case, as the attenuator is increased, 
the forward transmission never becomes infinite. Therefore, one con- 
cludes that the circuit is stable if and only if (12) never occurs. 


(12) 


3.1 Discussion of (11) 


Let the effect of p, on p’ be considered first. Since both N and 
[NV — (—1)"]/3 are always positive, one has from (11) that the mag- 


nitude of p’ is maximum when p, = —1. Therefore, 
The output reflection (p’) is maximum when 
|Z,|<«R,. (13) 


Since | p, | S 1, the first relation of (12) gives that instability can occur 
only if 


3 
| eo | = ap (aap = Po (14) 


* (RE) means “Real part of”. (IM) means ‘Imaginary part of”. 

{| The attenuator is supposed to be ideal, i.e., without phase delay. Note that 
the circuit is stable when the attenuator provides infinite attenuation. In fact, in 
this case | pz | = O and therefore (12) is never satisfied. 
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Therefore, 


In order to guarantee stability it is sufficient to require 

loz| < |pzo|. Note that, in a practical case, it will suffice to 
require (14) in the vicinity of the output carrier only.* That (15) 
is, consideration can be limited to those frequencies which 

fall within the pass-band of the output circuit. 


Now put (11) into the form 


(y(n, - P= EB) _ 2, = 0. (16) 


By considering the imaginary part of (16) one finds that (11) [and 
therefore (12), also] is satisfied only if either 


Pri = Por = 0 (17) 
or 
(—1)"" purpor > 0, (18) 


where pz; , p,; are the imaginary components of pz, p,. According to 
(17), (18) one can say 


Instability can occur only if at some frequencies p either 
one of the following two situations occur: (19) 


(2) both Z, , Zz are real. 
(71) X,X,(—1)""* > 0(N = 2"). 


Consider now the imaginary part of (11). After some manipulations 
one obtains 


N—(-"] _ —Np? 
pua| Nove x Na (ou! _ Pe (pyrpxt)- (20) 
PLI 


Therefore, if p,; , prr ¥ O, then (20), (18) give 
Ye N — (-1)” 
(—1) ozn| No —_ Y— Cu! > 0. (21) 


Note that, if p,r = pir = 0, then (21) follows directly from the first 
relation of (12). According to (21) one can say 

Instability can occur only if, at some frequencies, 

(21) is satisfied. (22) 


* See footnote on p. 806. 
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It is recalled that p,r > 0 or prr < 0 depending on whether or not 
|Z, | > Ry (and similarly for p,r , | Z,|, 21). The importance of (19), 
(22) follows from the fact that they allow the question of stability 
to be answered in many cases simply by looking at the signs of X, , 
X, and at the magnitudes of | Z, |, | Zz | (more precisely, of | Z, |/R: , 
|Z, |/Ry). This will be illustrated by the examples of the next section. 
Note that the presence of (—1)” in (11), (18), (21), follows from the 
fact that there is a reversal in the sign of both p", T’7 each time NV 
is increased by a factor of two. This is explained by Theorem 1 of the 
next section which shows that a multiplier of order 2” is equivalent 
to a chain of doublers. Since the round-trip transmission of a doubler 
is negative, each doubler gives a contribution of 180° to the phase of the 
overall round-trip transmission (from which (—1)” follows). 

Consider now some special cases. If N = 2, 4, then (11) becomes 


pip’ = (—1)"p,(Np, — 1) = 1. 
Thercfore, 


if N = 2 or4, instabilities can occur only if both the generator 
and the load are interacting (i.e., both pz , p, ¥ 0). (23) 


If N > 4, on the other hand, then | p’ | > 1 evenif | p, | = 0. Therefore, 


if N > 4, then instability can arise even if the generator 
does not interact (i.e., p, = 0). (24) 


This discussion of interactions will be concluded by emphasizing that 


if pp = O there is no interaction at all. Therefore, one may 
say that the load reflection p, is the primary cause of in- (25) 
stabilities. 


IV. ANALYSIS OF THE FREQUENCY DEPENDENCE OF THE PROPERTIES OF 
A MULTIPLIER OF THE TYPE (26) AND DISCUSSION OF THE STABILITY 
CONDITIONS 


Consider a multiplier which has the following characteristics: 


The order of multiplication is N = 2”. It has the minimum 
number (n — 1) of idlers. It is tuned. (26) 


The analysis of the frequency dependence of the properties of such 
a multiplier is based on the following theorem which is demonstrated 
in the Appendix. 


STABILITY AND NOISE 809 


Theorem 1:* <A multiplier of the type (26) is equivalent to the chain of 
stages shown in Fig. 4. The doublers of Fig. 4 are ‘‘ideal’’. More precisely, 
they are frequency independent and lossless, and they do not produce 
AM @ PM conversion. Their scattering matrices are given by (2) (with 
N = 2). That is, 








§ = a (Uf only PM 1s present) 

2 1 (27) 
2s 1 
ss * (If only AM is present). 

1 0 








LTH DOUBLER (L+1)TH DOUBLER 
(r=2') 





Fig. 4— Chain of doublers equivalent to the multiplier of order N = 2” shown 
in Fig. 5. 


Note that this theorem applies to the general case of a multiplier which 
is lossy and produces AM @ PM conversion. Fig. 5 shows the equivalent 
circuit of the actual multiplier from which the chain is derived. Z, 
is the impedance presented by the external circuit to the varactor. 
R, is the series resistance of the varactor. Cy is the varactor average 
capacitance. The generator V, consists of a carrier voltage V, of 
frequency w), and of noise terms 2,;, U,2., °**, Uy Which 
correspond to the various sidebands wy - p, 2wo + 7p, ete. 

The impedance Z, resonates with Cy at the ‘‘desired”’ carriers wo , 
209, °°: , Nwo. Furthermore, Z is so large for w far from wy , 2w9, --- , 
No that current flows through the varactor only at the desired carriers 
and their sidebands. 

Let now the chain of Fig. 4 be examined. Consider the impedance 

* In deriving (27) the approximation w) + p & w has been made. For more 


details, see the considerations made in the Appendix on the frequency dependence 
of the properties of the ‘‘ideal’’ doublers of Fig. 4. 
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Fig. 5 — Multiplier equivalent circuit. 


connecting the output of the zth doubler to the input of the (¢ + 1)th 
doubler (¢ < n). Current is flowing only at the frequencies rw , rwo - p 
(r = 2°), through this impedance. Furthermore, at w & rw (i.¢., wo = 
Two , TWH) + p) Z, gives the impedance presented by the rw,-idler circuit 
to the varactor. Therefore, the impedance connected to the output of 
the ith doubler corresponds to the 7th idler and can be represented by 
means of a series resonant circuit (resonant at rw,) as illustrated in Fig. 6. 
The resistance R, shown in Fig. 6 represents the losses of the idler circuit 
and it is equal to R, , the series resistance of the varactor, if the varactor 
is the only lossy element. The chain of Fig. 4 can be represented by 
the diagrammatic circuit of Fig. 7 in which, according to the preceding 
considerations, S;,, represents the zth idler circuit and can be approxi- 
mated as illustrated in Fig. 6. The impedance Z{‘*” is the characteristic 
impedance of the input port of the (¢ + 1)th doubler. It is equal to the 
impedance presented by the input port of the (z -+ 1)th doubler at the 
carrier frequency rw) . Z‘°) is the characteristic impedance of the output 
port of the 7th doubler and it is related to Zj;’ through the formula 
Zo) = ZO 4 R, (with r = 2°). Z{® is calculated in the Appendix. 

At w & No, , the impedance Z, is the impedance presented to the 


ORR oop 





le ee 


(i +14) TH INTERSTAGE 


TH 
DOUBLER 
(r=2l) NETWORK Si. 


Fig. 6— (¢ + 1)th interstage network Sts. 
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varactor by the output circuit of the multiplier. Accordingly, the 
impedance connected to the output of the last doubler of Fig. 4 can be 
represented as shown in Fig. 7, where S,,, corresponds to the output 
circuit of the multiplier. The input circuit of the multiplier is represented 
in Fig. 7 by means of a network S, connected between the generator 
and the first doubler. Note that S,, S,: include R, and Cy. Further- 
more, they are tuned at the carrier frequencies and therefore they 
provide unity transmission at the carriers if losses are absent. 





O 
x2 
(IDEAL) 
Oo 


UTH UTH 
INTERSTAGE DOUBLER 
NETWORK 


Fig. 7—Schematic representation of the chain of doublers equivalent to the 
multiplier of order N = 2”. 


4.1 Analysis of the Worst Case (No Losses, No AM = PM) 


Assume that the multiplier does not produce AM = PM conversion, 
and that it is lossless. (That is, S, , S2,--- , S,+; are lossless and do not 
produce AM = PM conversion.) 

Let the doubler be examined in detail. It is the most important 
multiplier because it represents the elementary constituent of any 
multiplier or chain of multipliers. 


4.2 Discussion of the Doubler 


According to the preceding discussion, the equivalent circuit of a 
doubler consists of three stages: an input stage S, , an “‘ideal’’ doubler, 
and an output stage S,. Therefore, the general considerations of the 
preceding section on the interactions between a multiplier and two 
other stages are applicable to the analysis of a doubler. Suppose that 
the scattering coefficients of the input and output stages are labeled 
1 and 2, respectively. Then the overall forward PM transmission is 
given by* 


* (28) can be obtained by means of standard techniques. See for instance Ref. 7. 
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5 2Ea s 

T;, = = — 

1 — p2 + 2pip2 
Since in this section it is assumed that S; does not produce AM = PM 
conversion and is lossless, its PM scattering parameters are simply 
given by its reflection and transmission coefficients at the frequency 
w + p(w = wp or 2w) depending on whether z = 1 or z = 2). The char- 
acteristic impedances of the two ports of S; are equal because S; is 
lossless; they will be denoted Z‘”. If, for instance, S; is either one of the 
two simple circuits of Fig. 8, then 


(28) 











(a) 


Fig. 8— Single-tuned resonant circuits. 





witht 
B. 
PS ees 
B; + yp 
pi = ps = ps = 41 — T,) (29) 
ZS? 
B= T°? 


where the negative sign applies to the case of Fig. 8(b). Note that B; 
is the 3-dB bandwidth of T; . Note, furthermore, that the only difference 
between the two circuits is that they have opposite reflections. If p ¥ 0, 
o , the signs of the real and imaginary components of the reflections are 


circuit (a) Pir > 0, pi > O 
circuit (b) Pir < O, pir < O, 


(30) 


+ The approximation w 2 w + p has been made in deriving (29). That is, the 
resonant circuit has been supposed to be very narrow band. 
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where p;r is the real component of p; , and p;; is the imaginary com- 
ponent. 

Suppose now that both S,, S, can be approximated by either one 
of the two circuits (a), (b) of Fig. 8. From (18), (21) one has that 
instability may occur only if 


PirP2r > 0, por(1 — 2 pir) > 0. (31) 


From the first of (81) and from (80) it follows that the circuit is unstable 
only if pizper > O. Therefore, by combining this result with the second 
inequality of (31), one has that instability requires 


Por > 0, Pir > 0. (32) 


Therefore, the only case in which instability may occur, is that in 
which both S,, S, are of the type (a). Let now this case be examined 
in detail (see Fig. 9). 


4.3 Discussion of the Case in Which Both S, and S, can be Approximated 
by Means of a Single-Tuned Series Resonant Circuit 


This case is most important for two main reasons. A first reason is 
that the equivalent circuit of a varactor includes the series connection 
of Cy and the inductance of the diode mount. Therefore, in most practical 
cases it will be possible to represent both the input and output circuit 
of a multiplier by means of a series resonant circuit, by first approxima- 
tion. The second reason is that, in the case of an idler circuit, it has 
already been pointed out that the equivalent circuit is a series resonant 
circuit, to a first approximation. 

Let the frequency p be normalized with respect to the output band- 
width (that is, B, = 2) and lety = B,/B, . From (28) and (29) one has 


en 27 
| gee peeer an eur aa ae 
me (y — 2p") + yp oe 


x2 
(IDEAL) 





Fig. 9— Case in which both the input and output circuit of the doubler can 
be opr by means of single-tuned series resonant circuits. 
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By analyzing (83) one finds that the behavior of T7 is different de- 
pending on whether y < } ory > 4. More precisely, 


(1) if y < 4, then maximum transmission (7'y;) occurs at p = 0; 
therefore, Ty, = 2; and 

(iz) ify > 4, then maximum transmission occurs at p = ~/(4y — 1)/8 
and it is given by 


- 27 
rns 
1/4 + jv (4y — 1)/8 


All this is illustrated in Fig. 10, where 20 log, |77 | is plotted as a 
function of p for different values of y(= B,/B,). One can see from Fig. 10 
that if the output bandwidth is enough greater than the input bandwidth 
(i.e. y < 4), then the transmission curve decreases monotonically with 
p. On the other hand, if y > 4, then 


a peak appears in the transmission curve, and the maximum 
transmission increases indefinitely with y. oe, (34) 


eo deed 
pt | AAA TL 
a 
ea eee 
Lhe ALN! 


Shee es 
BEAR 








20 LOG yo | Te| 
bh 






ae 10 — Effect of y (=B:/Bz) on the forward PM transmission of the doubler 
of Fig. 9. 
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a 11— Effect of y (=B:/Bz) on the input PM reflection of the doubler of 
ig. 9. 


This is in accordance with (2). One concludes from (84) that small 
values of B,/B, are desirable if large values of | 77 | are to be avoided. 
Small values of B,/B, are desirable also for the following stability con- 
siderations. If B,/B, is large, then the input PM reflection p; becomes 
larger than unity, as shown in Fig. 11. In addition, large values of B,/B, 
enhance the reverse overall transmission 7*;, as shown in Fig. 12. 
Consequently, one may conclude 


Large values of y(B,/B.) deteriorate the stability of the 
doubler and enhance the PM forward transmission (i.e., the 

output noise) of the multiplier in the vicinity of the carriers, (385) 
if both S, , S. are of the type shown in Fig. 8(a). 


It is important to note that the circuit becomes unstable only in the 
limiting case ~y = «©. Remember that it has been pointed out in the 
preliminary considerations of the first section that all scattering coeffi- 
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Fig. 12 — Effect of y (=B:/B:z) on the reverse PM transmission of the doubler 
of Fig. 9. 


cients become infinite when the multiplier is on the point of becoming 
unstable. 


4.3.1 Remaining cases 


In the remaining cases, [in which both S, , S, are of the type (a), (b) 
and at least one of S,, S» is of the type (b)], the overall transmission 
7; decreases monotonically with p for all values of B,/B,. This is in 
accordance with (32) and it can be directly verified by using (28), (29). 


4.4 Discussion of the General Case 


In most practical cases the output and input circuits of a multiplier 
are far more complicated than the simple cases considered in the preced- 
ing discussion. For instance, the impedance of the external circuit 
Z,. + 1/jwCo may have spurious resonances at the sidebands, etc. 
It is therefore important to consider the case in which the input and 
output circuits of a multiplier have more complicated configurations. 
Let first the input circuit be considered in detail. 
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4.5 Discussion of the Case in which only the Output Circuit can be 
Approximated by Means of a Single-Tuned Series Resonant Circurt 


According to the terminology of the first section, let Z, be the im- 
pedance presented by the input circuit (i.e., by the cascade connection 
of S, and R,) to the “‘ideal’”’ doubler of Fig. 18. Similarly, let Z, be the 
impedance presented by the output circuit to the output port of the 
‘deal’? doubler. 

Consider now the two cases illustrated in Fig. 14. The parameters 
B,, , By , B. represent the bandwidths of the various resonant elements 
(when they are considered individually). In order to investigate the 





Fig. 13 —Impedances Z, and Zz presented by the input and output circuits to 
the “ideal” doubler. 


possibility of instabilities, consider the reflections p,, p,, defined in 
the previous section (p, = p{, px = p32). Then 


ip 
= — iP _. 36 
PL B, a jp ( ) 
p, is different in the two cases (a), (b). In case (a), p, is given by 
==> —p' (Bra gine it) Nome jp(Bia Day eae 2p") . 
°  (BLBy — p (8B, + By) + ip2Bi. + Bis — 2p’)] 


Case (b) can be obtained simply by changing the sign of (37) and by 
interchanging B,,, By. 
Consider now the infinite gain condition* 


1 — pz + 2pzp, = 0 (38) 


which is obtained by letting the denominator of (28) equal zero. Those 
values of B,,, By, , B. for which (38) is satisfied at one or more fre- 


* Note that (38) is equal to (11) (n = 1, N = 2). 


(37) 
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quencies furnish the “‘boundary”’ of the region of stability (see the 
general considerations on stability of the first section). 


Consider now the circuit of Fig. 14(a). By substituting (86), (37), in 
(88) one finds that the circuit becomes unstable under certain conditions. 
These conditions are illustrated in Fig. 15. One can see from Fig. 15 
that, if both B,,/B,, , B./B,, are too small, then the circuit is unstable. 

If B,, > B,,., then the input circuit does not have spurious resonances, 





Fig. 14— Two typical examples in which only the output circuit can be ap- 
proximated by means of a single-tuned series resonant circuit. 


as one may verify from (37).* But Fig. 15 shows that the circuit may 
still be unstable. Therefore, one can say 


instability may occur even if the circuit does not have 
spurious resonances. (That is, even if the circuit impedance 
is real only at the desired carriers wo , 2wp .) (39) 


Fig. 15 shows that if B,, >> B,,, then stability is secured. This can 
be explained by noticing that for B,, >> B,, the circuit reduces to that 
of Fig. 9, which has already been found to be stable. 

Fig. 15 shows also that stability is secured if B, > 1.2B,, (approxi- 
mately). This can be explained as follows. Consider first the following 
theorem. 


*If By > Bi, then p,zy ~ Ofor0 <p < o. 
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Theorem 2: If at all frequencies p either one of the following two conditions 
lex | < | pro | (40) 
|Z, |>R, (41) 


is satisfied, then the doubler is stable. | pro | is 4 af losses are absent, and 
it is greater than 3 if losses are present. 


Proof: The preceding section has shown* that the first condition 
guarantees 7’; ~ o. Consider therefore, the second condition. Note 


2.5 





2.0 


STABLE 


0.5 1.0 1.5 2.0 
B2/Bia 


Fig. 15 — Stability region of the doubler of Fig. 14(a). 


that |Z, | >> R, is equivalent to p, & 1. Therefore, by substituting 
p, = 1 in (88) and neglecting the exceptional case p, = —1, one finds 
that (41) also guarantees 777 + o. One concludes therefore, that either 
one of (40), (41) guarantees T7 ~ o. Consider now a circuit which 
satisfies (40), (41). Next, decrease the amplitude of | pz, |. Since, as 
| pz, | is decreased, the forward transmission never becomes infinite 
one concludes that the circuit is stable, according to (4). 
Consider now the circuit of Fig. 14(a) and suppose B,, « B,. Note 
that condition (40) is satisfied for 
1 
p< 94/2 B,, (42) 


* See (14). Note that | pz | < 4 corresponds to the worst case (no losses; py = 1). 
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as one can verify from (86). Furthermore, since B,, « B, , one has that 
condition (41) is satisfied at the frequencies p 


p> a By > By. (43) 


Therefore, the circuit is stable, as shown by Fig. 15. 

The preceding demonstration used the fact that the condition of 
stability places restrictions on only those values of p, (i.e., of Z,) which 
occur at those frequencies p for which | p, | < | pzo |. This property 
is most important and will be further emphasized by stating 


the condition of stability places restrictions only at those 
frequencies p which are rejected by the output filter. There- 

fore, at those frequencies which are passed by the output (44) 
filter, the impedance Z, produced by the output filter can be 

quite arbitrary (all this applies to the case of a doubler). 


In this statement, the frequencies which are ‘‘rejected” by the output 
filter are given by the complement of (42) 


1 

> eB (45 
p 9 4/2 2 ) 
if losses are absent. The frequencies which ‘“‘pass through” the output 
filter are given by (42) 

1 
Bp (46 
22° 

Consider now the circuit of Fig. 14(b). Consider first the following 
theorem. 


p< 


Theorem 3: If R, > R, , then infinite transmission never occurs. 
Proof: Let Z, = R, + 7X, , and assume 

R,> h,. (47) 
Then, after some manipulations one finds 


(R, ag R,) i 2(k, at R,) = jX, 
(Rk, + R,) + 5X, 
Therefore, (38) is never satisfied for all values of p, . This demonstrates 


Theorem 3. 
Consider now the circuit of Fig. 14(b). Since X, > 0, from (19) 


|1 — 2p, | = = 1. 
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one has that infinite gain occurs only if 
xX, > 0. (48) 
But one can directly verify that the circuit of Fig. 14(b) has the fol- 
lowing properties 
R, > R, when X, > 0. 


Therefore, infinite transmission never occurs because of Theorem 3. 
This demonstrates that the circuit is always stable. Theorem 3 has a 
very important consequence. 


Spurious resonances which occur at the input sidebands 

do not produce instability, provided they produce high im- (49) 
pedances (R, > R,) at the varactor terminals (this is valid 

for a doubler). 


Note that (49) is valid for any arbitrary output circuit. 
4.6 Case in which the Input Circuit can be Approximated by Means 
of a Single-Tuned Series Resonant Circutt 


Suppose now that the input network consists of a series resonant 
circuit and that the output circuit is quite arbitrary (see Fig. 16). 
Then one has 


By. — Ip 
, = 1 _ 2 9 = als, 
Po By + ip 
Therefore, | p’ | = 1. Consequently instability may arise in the limiting 
case | p, | = 1. If one neglects this limiting case, one can say 


If the input circuit of the doubler can be approximated by a 
single-tuned resonant circuit in series with the generator, (50) 
then stability is secured for any arbitrary output circuit 
configuration (in the case of a doubler). 





Fig. 16 — Case in which only the input circuit can be approximated by means 
of a single-tuned series resonant circuit. 
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Note that (50) is in agreement with (13), (41) which show that it is in 
general desirable that | Z,| > R,. 


4.7 Discussion of a Frequency Multiplier of Order N > 2 (both a Single 
Stage and a Chain of Stages) 


The following discussion is concerned with the problem of designing 
a stable chain of doublers. This problem arises both in the design of a 
single stage multiplier of order 2" > 2, and in the design of a multiplier 
of order 2” consisting of many stages. That the chain may be unstable 
even if each doubler is individually stable, is a widely-known experi- 
mental fact. It will be shown here that if N > 4, then unstability may 
arise even if: each doubler is individually stable and it is, in addition, 
of the type shown in Fig. 9 (which is the simplest case which may occur). 

Consider an octupler and suppose that both its input and output 
circuits consist of single tuned resonant circuits connected in series 
to R, and R, . Then, according to Fig. 6, its equivalent circuit is that 
shown in Fig. 17 and the bandwidths B, , B. , B; , B, are the parameters 
on which the frequency dependence of 77; (and therefore, also the 
stability of the circuit) depends. To simplify the analysis assume the 
following conditions: 


2a SP Shs (61) 


and let the frequency p be normalized with respect to B,/2 (i.e., By = 2). 
In this way 7’; depends on the parameter y, only. T’7 has been cal- 
culated for different values of y and the results are shown in Fig. 18. 
One can see from Fig. 18 that the multiplier is unstable if 


y = 1.425 (approximately). (52) 


This result is in agreement with (35), which pointed out that small values 
of y are desirable. The conclusion from (52) is that 


x2 x2 x2 
(IDEAL) (IDEAL) (IDEAL) 





_Fig. 17 — Equivalent circuit of an octupler in which both the input and output 
circuits consists of single-tuned series resonant circuits. 
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if three or more doublers of the type shown in Fig. 9 are 
connected in cascade, the chain will be unstable unless each (53) 
stage is broadband enough with respect to the preceding one. 


Consider now a doubler which is individually stable, and suppose 
that the input circuit is made more narrowband than the output circuit. 
Then, if a stage is connected to the input port of the doubler, instability 
will not occur. In fact, according to (44), instability may arise only 
from interactions occurring at the frequencies which are rejected by 





20 LOG ara 





Fig. 18— Effect of y (=B:/Bis:) on the forward PM transmission of the octu- 
pler of Fig. 17. 


the output circuit. But these frequencies are also rejected by the input 
circuit. Therefore, the input filter prevents these ‘‘dangerous”’ interac- 
tions. Note that if a broadband* stage is connected to the output of 
the doubler, then instability does not occur because the output circuit 
of the doubler will reject those frequencies at which the broadband 
stage presents reflections. The conclusion can be stated as follows. 

Consider a chain consisting of doublers which are individually 

stable and which have B, «< B, (where B, , Bz are the input 

and output bandwidths of the general doubler, respectively). (54) 

The chain will be stable, provided the input circuit of each 

* Broadband with respect to the output circuit of the doubler. 


824 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1967 


doubler is broadband with respect to the output circuit of the 
preceding doubler. 


This is also in accordance with (52) and with the conclusions derived 
from it. 


4.8 Discussion of the Effect of AM <= PM conversion 

Consider a doubler and suppose, for instance that the output circuit 
S, produces AM = PM conversion. Then, if a unit PM wave is incident 
to the input port of S. , two waves will be reflected back, 


P(PM,PM) » P(PM,AM) 


where the subscript ( , ) indicates the type of incident wave and the 
type of reflected wave (in this order). Conservation of energy requires 


| P(PM,PM) | s1= | PL(PM,AM) | 

which shows that the amplitude of the reflected PM wave will in 
general be decreased by the presence of PM — AM conversion. This 
1s a desirable effect. Note that the AM wave generated by S, may be 
reflected back and converted again to PM by S,. However, for this 
to happen, the AM wave must travel through the “‘ideal’’ multiplier 
in the reverse direction, be reflected back by S,, and travel through 
the multiplier in the forward direction. Through this path the AM wave 
is attenuated because the AM round-trip (power) transmission of the 
‘ideal’? doubler is +, (see (1)). The conclusion is that, when the AM 
wave is converted back to PM by S; , it has small effect in most practical 
cases. Similar arguments can be applyed to the case in which S, has 
AM 2 PM conversion. The conclusion is 


input and output circuits which produce AM — PM con- 
version are in general preferable to those which do not (55) 
produce AM = PM conversion. 


Because of the qualitative nature of the discussion leading to (55) 
a more convincing argumentation may be desirable. A rigorous demon- 
stration requires considerably more space and is beyond the scope of 
this paper. 


V. CONCLUSIONS ON THE DESIGN OF A STABLE MULTIPLIER 
Consider first a doubler. From (35), (44), and (54), one has 


the output bandwidth should be broadband with respect 
to the input bandwidth. That is, (56) 


B,>B, . 


STABILITY AND NOISE 825 


The input circuit can have any arbitrary frequency dependence at 
those frequencies » which correspond to the passband of the output 
filter [see (44)]. In fact, this will not cause any instability. 

Instability may arise only 7f both the input circuit and the output 
circuit satisfy simultaneously certain conditions. For instance it must 


be 
XX 1, 2 0, R, < Ry ’ | PL | > 3 (57) 


for some frequencies in order that instability may occur. From (57) 
one can see that spurious resonances of the input circuit do not produce 
instability if they cause 


k, > kh, . 


Furthermore, in the design of the input and output circuits, those 
circuit configurations should be preferred which produce AM @ PM 
conversion. Therefore, low-pass circuit configurations at the input, 
and high-pass circuit configurations at the output, are preferable to 
bandpass circuit configurations. 

Finally, under certain circuit conditions, instability may arise even 
in the absence of spurious resonances. However, instability can be 
avoided by making the output circuit broadband enough with respect 
to the input circuit. This last part of the statement follows from (44). 

The need for condition (56) is further emphasized by considering 
the stability conditions for a chain of doublers (or for a multiplier of 
order N = 2”, which has been shown to be equivalent to a chain of 
doublers). 


A chain of doublers will be stable if the following conditions 

are satisfied: each doubler is individually stable and satisfies (58) 
(56); furthermore, each doubler is broadband enough with 

respect to the preceding one [see (54)]. 


If (58) is not satisfied, (53) shows that instability may arise even in 
the simplest case in which the input and output circuits of each doubler 
consist of (single-tuned) series resonant circuits. 

Note that (56), (57) etc. apply also to the design of a multiplier of 
order N = 2”. In this case, one should design the input circuit by con- 
sidering the doubler which consists of the cascade connection of: the 
input circuit, an ‘‘ideal’’ doubler, the first idler (connected in series to 
the output of the ‘“‘ideal’”’ doubler). The output circuit should satisfy 
the stability conditions of the doubler which consists of the last idler 
connected in series to the input of the ‘‘ideal doubler’’, and of the output 
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circuit connected to the output of the “ideal” doubler. Note that 
requirement (56) gives 


B,KB,K +++ < Bi <K Bas , (59) 


where B, is the input circuit bandwidth, B, is the first idler bandwidth 
[see (29) and the Appendix], B, is the last idler bandwidth, and B,,, 
is the output circuit bandwidth. 

It is important to emphasize that all these requirements have been 
derived in the case of no losses. The presence of losses reduces the 
limitations [such as (56) and (59)] placed by the condition of stability. 
If the losses are so large that the multiplier efficiency is less than 1/N 
[see (24) of Ref. 3] then the multiplier will be in general stable. 

Note that it has always been implied that the impedance presented 
to the varactor terminals by the external circuit* is “‘large’”’ at frequencies 
very far from the carriers. If this is not true, other types of instabilities 
may arise, such as subharmonic generation, etc. 


APPENDIX 


Equivalence of a Multiplier of Order N = 2” to a Chain of Doublers 


Consider a multiplier of the type defined by (26). Its equivalent 
circuit is shown in Fig. 19. V, is the input voltage generator of frequency 


Z=Zet1/jwCotRs 


Vg (~) 


Vg=Vgi t Vga t eee +Vgn (~) 


V=AQ*= F(Q) 


Fig. 19 — Multiplier equivalent circuit. 


Wo» Ugi» Vga y °** y Ven are the noise terms present at the various side- 
bands wo + p, 2m) + p, --: , Nwo & p. The series resistance R, and 
the average capacitance C, of the varactor are included in the impedance 
Z(w). Therefore, the nonlinear capacitance shown in Fig. 19 represents 
the nonlinear part of the capacitance of the varactor and it has a Q-V 


_ * The external circuit includes Co, the average capacitance of the varactor, the 
inductance of the varactor, etc. 
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characteristic of the type: V = F(Q) = AQ? (see the Appendix of 
Ref. 3 for more details). 
The hypothesis are 


Doi = Be oe Sot SAG oN) (60a) 
Z(w) = © forwfarfromre, (r = 2° = 1,2,---,N). (60b) 


In (60a) it is assumed that the circuit is resonant at the idler fre- 
quencies and at the input and output frequency. R, is the impedance 
of the input generator; R., --- , Ry. represent the losses of the idlers; 
Ry is the impedance of the load at Nw). In (60b) it is required that cur- 
rent flow be limited to the frequencies wy), 2w), --- , Nw and their 
sidebands. In this appendix both positive and negative frequencies 
will be considered, as illustrated by the Fourier Series (61), (62), etc. 
In (60) consideration is confined to w > 0 because the case w < 0 is 
given by Z(w) = Z(—w)*, where ( )* indicates the complex conjugate. 


When the sidebands are absent (i.e., when v,; = 0,2 = +++ = Uw = 0) 
the charge of the nonlinear capacitance is of the type 
Q(t) = 2 Q, exp (jrwot). F (61) 


8=0,°°°,n 


Notice that Q, = Q*, because Q(t) is real. The elastance S(¢) is given by 
S(t) = 2AQ() = >> S, exp (jraot) (62) 
r=+24 


s=0,°°¢.n 
with 
S, = 2AQ, . 63) 


Consider now the voltage V(¢) across the nonlinear eae 
From V = AQ’ and (61) one has 


Vi) =A 2 Q.Q, exp [jr + wot] = DO Vi exp [jlwot] 


(¢ = 42°; r= 42°; s,h =0,-:- ,n) 


from which one obtains the Fourier coefficient of V(t) relative to 
WO = TWo 


Vi = AlQ.2Q:2 + 2Q2Q--] (= £2°38 = 0,---,n). (64) 


The linear circuit connected to the nonlinear capacitance gives an 
additional expression for V, 


V, = —jroQ,Z (rw) + Vor (r = 42°;s =0,--- ,n), (65) 
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where V,, is the Fourier coefficient of V, relative to w = rw. There- 
fore, V,, = 0 for |r| # 1. From (65), (64) one obtains 


Vie jraoQ,Z (Tw) = A(Q,/2Q+/2 ae 2Q2,Q-,) 
(ry = +2°;58 =0,-:-- ,n). (66) 


Equation (66) gives the ‘‘equilibrium equations” of the circuit of Fig. 19 
at the carriers -kwy , -k2w9, °-- , Nar. 
Consider now the sidebands. Let 


QO) = 2 ars exp [ilreo + ip)d) (67) 


r=+2° 
t=41 


be the sidebands of the charge of the nonlinear capacitance. Then the 
voltage sidebands are given by 


v(t) = S(Hq(d). (68) 
By substituting (67), (62) in (68) one obtains », ; , the Fourier coefficient 
of v(t) relative to w = rwy + ip 
Ves = Qr/2,i8r72 + Qer,iS-r + Q-7, Sar (@ = +1). (69) 
For more details on the derivation of (69) see the Appendix of Ref. 3. 
The linear circuit gives a second expression for 2, ,; 
Vrs = Uy)r.s — Jo + tp)qr,:Zrwo + tp), (70) 


where (v,),,; indicates the Fourier coefficient of v, = v,, + ++: + v,y 
relative to w = rw -+ tp. By combining (69) with (70) one obtains 
the “equilibrium equation” of the circuit of Fig. 19 at the sidefrequency 
To + tp 


(Vy). = j(TwWo + tp) dr, :Z (ra + ip) 
+ Qr/2,iSr/2 + | PT ae -+- G6. 4927 : 


Consider now the chain illustrated in Fig. 20. Let the Q-V char- 
acteristic of the nonlinear capacitances of Fig. 20 be equal to that of 
Fig. 19. Furthermore, let Z,(w) be equal to Z(w) for w = rw, and be 
infinite for w far from rw, . That is, 


(71) 


Z,(w) = Zw) for w = Two , Tw + p 
Zw) = © for w farfrom rw (r = 2’ = 1,2,--- ,N), 


(72) 


where only positive frequencies are considered because Z,(w) = Z*(—w). 
Note that because of (72), the spectrum of the charge of the (s + 1)th 
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Zi Voz, Mae: 7, VgN 


F(Q) F(Q) ZN 







TW ,TWo tp 2PWo,2TWotP 
_—_—> —_ 


Zlerl| 


(S+1)TH VARACTOR (T= +25) 


Fig. 20 — Chain of doublers equivalent to the multiplier of Fig. 19. 


nonlinear capacitance of Fig. 20 is restricted to the carrier frequencies 
Two, 27wo (r = +2°) and to their sidebands rwy + p, 2rw) + p. Let the 
symbol ( )’ be used to distinguish the variables of the circuit of Fig. 19 
from those of the circuit of Fig. 20. Then the charge components of 
the (s + 1)th capacitance of Fig. 20 are Q/, Q3., q6.;, 3.2, with 
r= +2’ andi = 1. 

It will be shown that, under these hypothesis, the chain of Fig. 20 
is equivalent to the multiplier of Fig. 19. More precisely, it will be 
shown that the charges and voltages of the two circuits are equal. 


Demonstration: Let first the fact that Q,, --- , Qw are equal to Q/, 
-, Q4 be demonstrated. Consider therefore, the ‘‘equilibrium equa- 
tion” of the chain of Fig. 20 for w = rw (r = +2*). With reference 
to Fig. 21, it is obtained by applying Kirchoff’s law to the (s + 1)th 
loop. 
Across the first capacitance of Fig. 21, the voltage component of 
frequency 7w is 


AQT Qi 2 (ry = +2°, |r| 41, N). (73) 


The voltage component of frequency rw) produced by the second capac- 
itance is 


2AQ.,Q; (r= +2°, |r| #1,N). (74) 


Since only the current jrw,Q/ is flowing through Z,,,, the voltage 
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Ir=jrwoQi 
— 





Fig. 21 — Sth loop (r = 2°) of the chain of Fig. 20. 


across Z|, 1S 
FQ | - 1 (To) = Fro: Z (rao) = (7 = 2’, |r| #1,N). (75) 


Equilibrium requires that the sum of these three voltages be zero. 
Therefore, one obtains 

0 = Fro: Z(rwo) + AlQt Qt. + 2Q3,Q/,] (76) 
which is identical to the “equilibrium equation” of the circuit of Fig. 19 
for w = Two (with r = +2’, |r| #1, N). 

In a completely similar way, one finds that also for |r| = 1 and 
|r| = N the “equilibrium equations” of the two circuits are identical. 
If r = 1, for instance, the loop consists of the voltage generator V, 
in series with Z, and the first nonlinear capacitance of Fig. 20. One 
obtains 


Var oF Je0Qi Z (wo) = 2AQ20'; 
which is identical to (66) for r = 1. The case r = N is obtained by 
considering the last loop and one finds 

—jNwQnZ(Nowo) = AQ 2Qh2 


which is identical to (66) for r = N. All this demonstrates that Q,, --- , 
Qy are identical to Q/, --- , Q{. Therefore, the two circuits are equiv- 
alent at the carrier frequencies. 

Consider now the sidebands. Notice that the elastance components 
of the (s + 1)th nonlinear capacitance of Fig. 20 are S, , S:, (r = +2’). 
Therefore, the equilibrium of the (s + 1)th loop for w = rw) + tp 
requires 


(CA ae = j(Tw9 + tp) qr, :Z (Ta aie tp) (77) 
Q7/2,iSr/2 + SerQ.r,i a S_-Q3r,i (r= +2’, i= +1). 


STABILITY AND NOISE 831 


The first term of the second member represents the voltage component 
across Z,; , the second term corresponds to the voltage across the first 
capacitance, and the last two terms correspond to the voltage across 
the second capacitance. Equation (77) is equivalent to (71). Therefore, 
the two sets of charge components (q,,;), (@.;) (|r| = 1, 2, --- , N; 
t = +41) are identical. This concludes the demonstration of the equiv- 
alence of the two circuits of Figs. 19 and 20. Note that the voltage 
across Z;,; 18 equal to the voltage components of Z due to the frequencies 
--rw, and their sidebands. 


A.1 Discussion of the Circuit of Fig. 20 


The circuit of Fig. 20 can be represented by the chain of doublers 
illustrated in Fig. 22. The equivalent circuit of the (s + 1)th doubler 
of Fig. 22 is shown in Fig. 23. It consists of an “‘ideal’”’ input filter F{¢*” 
connected in series to the input, an “ideal” output filter PF“&*” con- 
nected in series to the output, and an “‘ideal’’ lossless varactor which 
has a Q-V characteristic of the type V = AQ’. The filter F{%*” has 
zero impedance at the input frequencies rw) , rw) + p and it has infinite 
impedance for w far from rwo (r = +2°). Similarly, F&%*” limits the 
output current to the frequencies 2rw) , 2rw) + p. Consequently, the 
series connection of FS. , Z, F{*” is equivalent to Z,,; , according to 
(72). Therefore, the two circuits of Figs. 22 and 20 are equivalent. 
By first approximation, the properties of the “‘ideal’’ doublers of Fig. 
22 are independent of the modulation frequency p. More precisely, the 
frequency dependence of the ‘“‘ideal’’ doublers can be neglected with 
respect to the frequency dependence of the impedances Z,, Z., etc. 
This is more precisely explained by the following considerations. 

Consider (77). Let 2,,; indicate the Fourier coefficient of the current 
flowing through Z,,,; for w = rw, + tp. Then 





Fig. 22 — Chain of doublers equivalent to the circuit of Fig. 20. 
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1,2 = Jao + tp)q,,: o> ceM). (78) 


If one approximates (78) by means of 7,,; = jrwog,,; , one obtains from 


(77) 


Z (rar Up)t,, i — S,/2br/2, i oe J Beste = iL S_+lor,s 
r r 2r 
(79) 
= (Ug)r,3 ° 


Equation (79) gives the approximate version of (77) which is obtained 
by neglecting the frequency (p) dependence of the behavior of the 
nonlinear capacitances of Fig. 20. In fact, one can see that the only term 
in (79) which depends on p is the first term, and this term is caused by 
the impedance Z,,, of Fig. 20. In this analysis the approximate expression 
(79) is valid because it is assumed Z(w) & © for w far from rw, . There- 
fore, at those frequencies p for which the approximation 7,,; = rwod,,; 
becomes invalid, the only important term in (79) is the first, which does 
not contain any approximation. However, it is important to notice 
that (79) is exact only in the limiting case of a narrow-band multiplier. 
That is, when Z,,,;(w) varies so rapidly with p that it is infinite outside 
a very narrow band around w = rw. 

According to the preceding considerations, the properties of the 
“ideal” doublers of Fig. 22 can be assumed independent of the modula- 
tion frequency :, by first approximation. 

In the following part of this Appendix, it will be shown that the 
impedance presented to the output of each “‘ideal’’ doubler by the 
following part of the chain is real at the carrier frequency. Furthermore, 
cach doubler is lossless. Therefore, according to the results derived in 
Ref. 3, one concludes that each doubler does not produce AM = PM 


TW 2rwo Be? (S+1) 


P : Fo UT 





V=F(Q)= AQ? 


(S+1)TH DOUBLER 
(r =+28) 


Z 


Fs 





Fig. 23 — Constituents of the chain of Fig. 22. 
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and that it is characterized by the two scattering matrices given by 
(27). The characteristic impedances of the input and output ports of 
the various doublers are derived in the following part. 


A.2 Calculation of the Output and Input Impedances of the Doublers 
of Fig. 22 
At the carrier frequency 7) , the input impedance of the (s + 1)th 
doubler (r = 2’) is given by 
Vie = 2AQ_,Qo, 
jroo, jroo, 
where V{&:t?, the input voltage of the (s + 1)th doubler, has been 


obtained from (74). 
The output impedance of the (s + 1)th doubler is given by 


Zoue = Z(2re) + Zn? (ry = +2"). (81) 





(s+1) _ 
Zin = 


(ry = +2°), (80) 


But one also has 


=e. Vi —AQ Q 
(stl) __ out = rr = 8 
Lieut a j2raQo, 2reQeo, (r +2 ys (82) 


where V ‘3? is the output voltage of the (s + 1)th doubler and it is 


out 


given by (73) with r/2 replaced with r. Therefore, from (80) and (82) 


one obtains 
2 2 
gage - AIOE 4 8 yon, (83) 


2 
r Wo 


By combining (81) and (88), one has 


ZEA 4 Zr) = ASL - q. (7 = +2'). (84) 


in as 
The output impedance of the last doubler (s = n) is 
Zone = Z(Noo) = Ry . (85) 


Equation (85) and the recurrent formulas (81) and (84) allow the output 
and input impedances of the various doublers to be readily calculated 
(Ry , A, | Q, | are known). 

Note that Z is real at the frequencies wy , 20), --- , Nw, because 
the circuit is resonant at these frequencies. Z(2w.), -++ , Z(N/2w) 
represent the idler losses. If, for instance, the varactor is the only lossy 
element of the circuit, then . 


VACTR EAC ee ey | (86) 
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where R, is the series resistance of the varactor. Therefore, since (85), 
(84), and (81) are real, one concludes that Z‘%, 23) are real. 

Finally, according to the convention of the first paper, the input 
and output characteristic impedances of the sth doubler are given by 
Zi, Z2., respectively. 


in ) out ) 
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